Background

Methods for inferring the three-dimensional (3D) configuration of chromatin from conformation capture assays that provide strictly pairwise interactions, notably Hi-C, utilize the attendant contact matrix as input. More recent assays, in particular split-pool recognition of interactions by tag extension (SPRITE), capture multi-way interactions instead of solely pairwise contacts. These assays yield contacts that straddle appreciably greater genomic distances than Hi-C, in addition to instances of exceptionally high-order chromatin interaction. Such attributes are anticipated to be consequential with respect to 3D genome reconstruction, a task yet to be undertaken with multi-way contact data. However, performing such 3D reconstruction using distance-based reconstruction techniques requires framing multi-way contacts as (pairwise) distances. Comparing approaches for so doing, and assessing the resultant impact of long-range and multi-way contacts, are the objectives of this study.

Results

We obtained 3D reconstructions via multi-dimensional scaling under a variety of weighting schemes for mapping SPRITE multi-way contacts to pairwise distances. Resultant configurations were compared following Procrustes alignment and relationships were assessed between associated Procrustes root mean square errors and key features such as the extent of multi-way and/or long-range contacts. We found that these features had surprisingly limited influence on 3D reconstruction, a finding we attribute to their influence being diminished by the preponderance of pairwise contacts.

Conclusion

Distance-based 3D genome reconstruction using SPRITE multi-way contact data is not appreciably affected by the weighting scheme used to convert multi-way interactions to pairwise distances.

Background

Comprehensive molecular profiling has revealed somatic variations in cancer at genomic, epigenomic, transcriptomic, and proteomic levels. The accumulating data has shown clearly that molecular phenotypes of cancer are complex and influenced by a multitude of factors. Conventional unsupervised clustering applied to a large patient population is inevitably driven by the dominant variation from major factors such as cell-of-origin or histology. Translation of these data into clinical relevance requires more effective extraction of information directly associated with patient outcome.

Methods

Drawing from ideas in supervised text classification, we developed survClust, an outcome-weighted clustering algorithm for integrative molecular stratification focusing on patient survival. survClust was performed on 18 cancer types across multiple data modalities including somatic mutation, DNA copy number, DNA methylation, and mRNA, miRNA, and protein expression from the Cancer Genome Atlas study to identify novel prognostic subtypes.

Results

Our analysis identified the prognostic role of high tumor mutation burden with concurrently high CD8 T cell immune marker expression and the aggressive clinical behavior associated with CDKN2A deletion across cancer types. Visualization of somatic alterations, at a genome-wide scale (total mutation burden, mutational signature, fraction genome altered) and at the individual gene level, using circomap further revealed indolent versus aggressive subgroups in a pan-cancer setting.

Conclusions

Our analysis has revealed prognostic molecular subtypes not previously identified by unsupervised clustering. The algorithm and tools we developed have direct utility toward patient stratification based on tumor genomics to inform clinical decision-making. The survClust software tool is available at https://github.com/arorarshi/survClust .

Gene regulation during cell-cycle progression is an intricately choreographed process, ensuring accurate DNA replication and division. However, the translational landscape of gene expression underlying cell-cycle progression remains largely unknown. Employing genome-wide ribosome profiling, we uncover widespread translational regulation of hundreds of mRNAs serving as an unexpected mechanism for gene regulation underlying cell-cycle progression. A striking example is the S phase translational regulation of RICTOR, which is associated with cell cycle-dependent activation of mammalian target of rapamycin complex 2 (mTORC2) signaling and accurate cell-cycle progression. We further identified unappreciated coordination in translational control of mRNAs within molecular complexes dedicated to cell-cycle progression, lipid metabolism, and genome integrity. This includes the majority of mRNAs comprising the cohesin and condensin complexes responsible for maintaining genome organization, which are coordinately translated during specific cell cycle phases via their 5' UTRs. Our findings illuminate the prevalence and dynamic nature of translational regulation underlying the mammalian cell cycle.

Motivation

The translational landscape of diverse cellular systems remains largely uncharacterized. A detailed understanding of the control of gene expression at the level of messenger RNA translation is vital to elucidating a systems-level view of complex molecular programs in the cell. Establishing the degree to which such post-transcriptional regulation can mediate specific phenotypes is similarly critical to elucidating the molecular pathogenesis of diseases such as cancer. Recently, methods for massively parallel sequencing of ribosome-bound fragments of messenger RNA have begun to uncover genome-wide translational control at codon resolution. Despite its promise for deeply characterizing mammalian proteomes, few analytical methods exist for the comprehensive analysis of this paired RNA and ribosome data.

Results

We describe the Babel framework, an analytical methodology for assessing the significance of changes in translational regulation within cells and between conditions. This approach facilitates the analysis of translation genome-wide while allowing statistically principled gene-level inference. Babel is based on an errors-in-variables regression model that uses the negative binomial distribution and draws inference using a parametric bootstrap approach. We demonstrate the operating characteristics of Babel on simulated data and use its gene-level inference to extend prior analyses significantly, discovering new translationally regulated modules under mammalian target of rapamycin (mTOR) pathway signaling control.

Common environmental contaminants such as bisphenols and phthalates and persistent contaminants such as polychlorinated biphenyls are thought to influence tissue homeostasis and carcinogenesis by acting as disrupters of endocrine function. In this study we investigated the direct effects of exposure to bisphenol A (BPA), mono-n-butyl phthalate (Pht), and polychlorinated biphenyl 153 (PCB153) on the proteome of primary organotypic cultures of the mouse mammary gland. At low-nanomolar doses each of these agents induced distinct effects on the proteomes of these cultures. Although BPA treatment produced effects that were similar to those induced by estradiol, there were some notable differences, including a reduction in the abundance of retinoblastoma-associated protein and increases in the Rho GTPases Ras-related C3 botulinum toxin substrate 1 (Rac1) and cell division cycle protein CDC42. Both Pht and PCB153 induced changes that were distinct from those induced by estrogen, including decreased levels of the transcriptional corepressor C-terminal binding protein 1. Interestingly, the three chemicals appeared to alter the abundance of distinct splice forms of many proteins as well as the abundance of several proteins that regulate RNA splicing. Our combined results indicate that the three classes of chemical have distinct effects on the proteome of normal mouse mammary cultures, some estrogen-like but most estrogen independent, that influence diverse biological processes including apoptosis, cell adhesion, and proliferation.

Large-scale cancer genome projects, such as the Cancer Genome Atlas (TCGA) project, are comprehensive molecular characterization efforts to accelerate our understanding of cancer biology and the discovery of new therapeutic targets. The accumulating wealth of multidimensional data provides a new paradigm for important research problems including cancer subtype discovery. The current standard approach relies on separate clustering analyses followed by manual integration. Results can be highly data type dependent, restricting the ability to discover new insights from multidimensional data. In this study, we present an integrative subtype analysis of the TCGA glioblastoma (GBM) data set. Our analysis revealed new insights through integrated subtype characterization. We found three distinct integrated tumor subtypes. Subtype 1 lacks the classical GBM events of chr 7 gain and chr 10 loss. This subclass is enriched for the G-CIMP phenotype and shows hypermethylation of genes involved in brain development and neuronal differentiation. The tumors in this subclass display a Proneural expression profile. Subtype 2 is characterized by a near complete association with EGFR amplification, overrepresentation of promoter methylation of homeobox and G-protein signaling genes, and a Classical expression profile. Subtype 3 is characterized by NF1 and PTEN alterations and exhibits a Mesenchymal-like expression profile. The data analysis workflow we propose provides a unified and computationally scalable framework to harness the full potential of large-scale integrated cancer genomic data for integrative subtype discovery.

Germ Cell Tumors (GCT) have a high cure rate, but we currently lack the ability to accurately identify the small subset of patients who will die from their disease. We used a combined genomic and expression profiling approach to identify genomic regions and underlying genes that are predictive of outcome in GCT patients. We performed array-based comparative genomic hybridization (CGH) on 53 non-seminomatous GCTs (NSGCTs) treated with cisplatin based chemotherapy and defined altered genomic regions using Circular Binary Segmentation. We identified 14 regions associated with two year disease-free survival (2yDFS) and 16 regions associated with five year disease-specific survival (5yDSS). From corresponding expression data, we identified 101 probe sets that showed significant changes in expression. We built several models based on these differentially expressed genes, then tested them in an independent validation set of 54 NSGCTs. These predictive models correctly classified outcome in 64-79.6% of patients in the validation set, depending on the endpoint utilized. Survival analysis demonstrated a significant separation of patients with good versus poor predicted outcome when using a combined gene set model. Multivariate analysis using clinical risk classification with the combined gene model indicated that they were independent prognostic markers. This novel set of predictive genes from altered genomic regions is almost entirely independent of our previously identified set of predictive genes for patients with NSGCTs. These genes may aid in the identification of the small subset of patients who are at high risk of poor outcome.

Recent prostate-specific antigen-based screening trials indicate an urgent need for novel and noninvasive biomarker identification strategies to improve the prediction of prostate cancer behavior. Noncoding microRNAs (miRNA) in the serum and plasma have been shown to have potential as noninvasive markers for physiologic and pathologic conditions. To identify serum miRNAs that diagnose and correlate with the prognosis of prostate cancer, we developed a multiplex quantitative reverse transcription PCR method involving the purification of multiplex PCR products followed by uniplex analysis on a microfluidics chip to evaluate 384 human miRNAs. Using Dgcr8 and Dicer knockout (small RNA-deficient) mouse ES cells as the benchmark, we confirmed the validity of our technique and uncovered a considerable lack of accuracy in previously published methods. Profiling 48 sera from healthy men and untreated prostate cancer patients with differing CAPRA scores, we identified miRNA signatures that allow us to diagnose cancer patients and correlate with a prognosis. These serum signatures include oncogenic and tumor-suppressive miRNAs, suggesting functional roles in prostate cancer progression.