Mouse and rat models are mainstays in pharmacology, toxicology and drug development -- but differences between strains and between species complicate data interpretation and application to human health. Dioxin-like polyhalogenated aromatic hydrocarbons represent a major class of environmentally and economically relevant toxicants. In mammals dioxin exposure leads to a broad spectrum of adverse affects, including hepatotoxicity of varying severity. Several studies have shown that dioxins extensively alter hepatic mRNA levels. Surprisingly, though, analysis of a limited portion of the transcriptome revealed that rat and mouse responses diverge greatly (Boverhof et al. Toxicol Sci 94:398-416, 2006).

Teneurin C-terminal associated peptide (TCAP)-1 is a member of a novel family of neuropeptides that has been highly conserved throughout evolution. TCAP-1 is expressed in the limbic system in areas such as the hippocampus and amygdala. In vitro, TCAP-1 increases cytoskeletal proteins in immortalized neurons and modulates neurite outgrowth in cultured primary hippocampal neurons. In vivo, TCAP-1 blocks stress-induced c-Fos in the hippocampus and amygdala, and modulates stress-induced anxiety-like behaviors. This suggests that TCAP-1 plays a role in the remodeling of limbic system networks to alter stress behaviors. Dendritic spines on the apical and basilar shafts of hippocampal neurons are sensitive to stress and many receive incoming excitatory synaptic connections. In this study, repeated daily injection of TCAP-1 for 10 days increased spine density in the CA1 and CA3 regions of the hippocampus without affecting spine density in the amygdala. Further investigation of the CA3 region indicated that TCAP-1 did not affect the morphology of apical dendrites, but decreased branching in the basilar dendrites 90-130 μm away from the soma. Moreover, TCAP-1 treatment increased open arm time and decreased closed arm entries on the elevated plus maze, a test of anxiety-like behavior. These results suggest that TCAP-1 may be associated with anxiety-like behavior via regulation of dendritic morphology in the hippocampus, independent of amygdalar modification.

MicroRNAs are small regulatory RNAs that post transcriptionally control gene expression. Reduced expression of DICER, the enzyme involved in microRNA processing, is frequently observed in cancer and is associated with poor clinical outcome in various malignancies. Yet, the underlying mechanisms are not well understood. Here we identify tumour hypoxia as a regulator of DICER expression in large cohorts of breast cancer patients. We show that DICER expression is suppressed by hypoxia through an epigenetic mechanism that involves inhibition of oxygen-dependent H3K27me3 demethylases KDM6A/B and results in silencing of the DICER promoter. Subsequently, reduced miRNA processing leads to derepression of the miR-200 target ZEB1, stimulates the epithelial to mesenchymal transition and ultimately results in the acquisition of stem cell phenotypes in human mammary epithelial cells. Our study uncovers a previously unknown relationship between oxygen-sensitive epigenetic regulators, miRNA biogenesis and tumour stem cell phenotypes that may underlie poor outcome in breast cancer.

The complexities of tumor genomes are rapidly being uncovered, but how they are regulated into functional proteomes remains poorly understood. Standard proteomics workflows use databases of known proteins, but these databases do not capture the uniqueness of the cancer transcriptome, with its point mutations, unusual splice variants and gene fusions. Onco-proteogenomics integrates mass spectrometry-generated data with genomic information to identify tumor-specific peptides. Linking tumor-derived DNA, RNA and protein measurements into a central-dogma perspective has the potential to improve our understanding of cancer biology.

The reproducibility of transcriptomic biomarkers across datasets remains poor, limiting clinical application. We and others have suggested that this is in-part caused by differential error-structure between datasets, and their incomplete removal by pre-processing algorithms.

Numerous recent studies have suggested that phenotypic effects of DNA sequence variants can be mediated or modulated by their epigenetic marks, such as allele-skewed DNA modification (ASM). Using Affymetrix SNP microarrays, we performed a comprehensive search of ASM effects in human post-mortem brain and sperm samples (total n = 256) from individuals with major psychosis and control individuals. Depending on the phenotypic category of the brain samples, 1.4%-7.5% of interrogated SNPs exhibited ASM effects. Next, we investigated ASM in the context of genetic studies of schizophrenia and detected that brain ASM SNPs were significantly overrepresented among sub-threshold SNPs from a schizophrenia genome-wide association study (GWAS). Brain ASM SNPs showed a much stronger enrichment in a schizophrenia GWAS than in 17 large GWASs of non-psychiatric diseases and traits, arguing that ASM effects are at least partially tissue specific. Studies of germline and control brain ASM SNPs supported a causal association between ASM and schizophrenia. Finally, significantly higher proportions of ASM SNPs than of non-ASM SNPs were detected at loci exhibiting epigenetic signatures of enhancers and promoters, and they were overrepresented within transcription factor binding regions and DNase I hypersensitive sites. All of these findings collectively indicate that ASM SNPs should be prioritized in follow-up GWASs.

Tumor heterogeneity is a limiting factor in cancer treatment and in the discovery of biomarkers to personalize it. We describe a computational purification tool, ISOpure, to directly address the effects of variable normal tissue contamination in clinical tumor specimens. ISOpure uses a set of tumor expression profiles and a panel of healthy tissue expression profiles to generate a purified cancer profile for each tumor sample and an estimate of the proportion of RNA originating from cancerous cells. Applying ISOpure before identifying gene signatures leads to significant improvements in the prediction of prognosis and other clinical variables in lung and prostate cancer.

The aneurysm clip impact-compression model of spinal cord injury (SCI) is a standard injury model in animals that closely mimics the primary mechanism of most human injuries: acute impact and persisting compression. Its histo-pathological and behavioural outcomes are extensively similar to human SCI. To understand the distinct molecular events underlying this injury model we analyzed global mRNA abundance changes during the acute, subacute and chronic stages of a moderate to severe injury to the rat spinal cord.

Most proteins destined for the extracellular space require disulfide bonds for folding and stability. Disulfide bonds are introduced co- and post-translationally in endoplasmic reticulum (ER) cargo in a redox relay that requires a terminal electron acceptor. Oxygen can serve as the electron acceptor in vitro, but its role in vivo remains unknown. Hypoxia causes ER stress, suggesting a role for oxygen in protein folding. Here we demonstrate the existence of two phases of disulfide bond formation in living mammalian cells, with differential requirements for oxygen. Disulfide bonds introduced rapidly during protein synthesis can occur without oxygen, whereas those introduced during post-translational folding or isomerization are oxygen dependent. Other protein maturation processes in the secretory pathway, including ER-localized N-linked glycosylation, glycan trimming, Golgi-localized complex glycosylation, and protein transport, occur independently of oxygen availability. These results suggest that an alternative electron acceptor is available transiently during an initial phase of disulfide bond formation and that post-translational oxygen-dependent disulfide bond formation causes hypoxia-induced ER stress.

2,3,7,8 Tetrachlorodibenzo-p-dioxin (TCDD) is an aromatic, long-lived environmental contaminant. While the pathogenesis of TCDD-induced toxicity is poorly understood, it has been shown that the aryl hydrocarbon receptor (AHR) is required. However, the specific transcriptomic changes that lead to toxic outcomes have not yet been identified. We previously identified a panel of 33 genes that respond to TCDD treatment in two TCDD-sensitive rodent species. To identify genes involved in the onset of hepatic toxicity, we explored 25 of these in-depth using liver from two rat strains: the TCDD-resistant Han/Wistar (H/W) and the TCDD-sensitive Long-Evans (L-E). Time course and dose-response analyses of mRNA abundance following TCDD insult indicate that eight genes are similarly regulated in livers of both strains of rat, suggesting that they are not central to the severe L-E-specific TCDD-induced toxicities. The remaining 17 genes exhibited various divergent mRNA abundances between L-E and H/W strains after TCDD treatment. Several genes displayed a biphasic response where the initial response to TCDD treatment was followed by a secondary response, usually of larger magnitude in L-E liver. This secondary response was most often an exaggeration of the original TCDD-induced response. Only cytochrome b5 type A (microsomal) (Cyb5a) had equivalent TCDD sensitivity to the prototypic AHR-responsive cytochrome P450, family 1, subfamily a, polypeptide 1 (Cyp1a1), while six genes were less sensitive. Four genes showed an early inter-strain difference that was sustained throughout most of the time course (atypical chemokine receptor 3 (Ackr3), collagen, type XVIII, alpha 1 (Col18a1), Cyb5a and glutamate dehydrogenase 1 (Glud1)), and of those genes examined in this study, are most likely to represent genes involved in the pathogenesis of TCDD-induced hepatotoxicity in L-E rats.

Recent evidence demonstrated CIN4 as a predictive marker of anthracycline benefit in early breast cancer. An analysis of the NCIC CTG MA.21 clinical trial was performed to test the role of existing CIN gene expression signatures as prognostic and predictive markers in the context of taxane based chemotherapy.RNA was extracted from patients in cyclophosphamide, epirubicin and flurouracil (CEF) and epirubicin, cyclophosphamide and paclitaxel (EC/T) arms of the NCIC CTG MA.21 trial and analysed using NanoString technology.After multivariate analysis both high CIN25 and CIN70 score was significantly associated with an increased in RFS (HR 1.76, 95%CI 1.07-2.86, p=0.0018 and HR 1.59, 95%CI 1.12-2.25, p=0.0096 respectively). Patients whose tumours had low CIN4 gene expression scores were associated with an increase in RFS (HR: 0.64, 95% CI 0.39-1.03, p=0.06) when treated with EC/T compared to patients treated with CEF.In conclusion we have demonstrated CIN25 and CIN70 as prognostic markers in breast cancer and that CIN4 is a potential predictive maker of benefit from taxane treatment.

The cancer transcriptome is remarkably complex, including low-abundance transcripts, many not polyadenylated. To fully characterize the transcriptome of localized prostate cancer, we performed ultra-deep total RNA-seq on 144 tumors with rich clinical annotation. This revealed a linear transcriptomic subtype associated with the aggressive intraductal carcinoma sub-histology and a fusion profile that differentiates localized from metastatic disease. Analysis of back-splicing events showed widespread RNA circularization, with the average tumor expressing 7,232 circular RNAs (circRNAs). The degree of circRNA production was correlated to disease progression in multiple patient cohorts. Loss-of-function screening identified 11.3% of highly abundant circRNAs as essential for cell proliferation; for ∼90% of these, their parental linear transcripts were not essential. Individual circRNAs can have distinct functions, with circCSNK1G3 promoting cell growth by interacting with miR-181. These data advocate for adoption of ultra-deep RNA-seq without poly-A selection to interrogate both linear and circular transcriptomes.

Recurrence of high-grade prostate cancer after radiotherapy is a significant clinical problem, resulting in increased morbidity and reduced patient survival. The molecular mechanisms of radiation resistance are being elucidated through the study of microRNA (miR) that negatively regulate gene expression. We performed bioinformatics analyses of The Cancer Genome Atlas (TCGA) dataset to evaluate the association between miR-106a and its putative target lipopolysaccharide-induced TNF-α factor (LITAF) in prostate cancer. We characterized the function of miR-106a through in vitro and in vivo experiments and employed transcriptomic analysis, western blotting, and 3'UTR luciferase assays to establish LITAF as a bona fide target of miR-106a. Using our well-characterized radiation-resistant cell lines, we identified that miR-106a was overexpressed in radiation-resistant cells compared to parental cells. In the TCGA, miR-106a was significantly elevated in high-grade human prostate tumors relative to intermediate- and low-grade specimens. An inverse correlation was seen with its target, LITAF. Furthermore, high miR-106a and low LITAF expression predict for biochemical recurrence at 5 years after radical prostatectomy. miR-106a overexpression conferred radioresistance by increasing proliferation and reducing senescence, and this was phenocopied by knockdown of LITAF. For the first time, we describe a role for miRNA in upregulating ATM expression. LITAF, not previously attributed to radiation response, mediates this interaction. This route of cancer radioresistance can be overcome using the specific ATM kinase inhibitor, KU-55933. Our research provides the first report of miR-106a and LITAF in prostate cancer radiation resistance and high-grade disease, and presents a viable therapeutic strategy that may ultimately improve patient outcomes.

Onco-proteogenomics aims to understand how changes in a cancer's genome influences its proteome. One challenge in integrating these molecular data is the identification of aberrant protein products from mass-spectrometry (MS) datasets, as traditional proteomic analyses only identify proteins from a reference sequence database.

The International Society of Nephrology/Renal Pathology Society (ISN/RPS) lupus nephritis (LN) classification is under reconsideration, given challenges with inter-rater reliability and resultant inconsistent relationship with treatment response. Integration of molecular classifiers into histologic evaluation can improve diagnostic precision and identify therapeutic targets. This study described the relationship between histological and molecular phenotypes and clinical responses in LN. Renal compartmental mRNA abundance was measured in 54 biopsy specimens from LN patients and correlated to ISN/RPS classification and individual histologic lesions. A subset of transcripts was also evaluated in sequential biopsies of a separate longitudinal cohort of 36 patients with paired samples obtained at the time of flare and at follow up. Unsupervised clustering based on mRNA abundance did not demonstrate a relationship with the (ISN/RPS) classification, nor did univariate statistical analysis. Exploratory analyses suggested a correlation with individual histologic lesions. Glomerular FN1 (fibronectin), SPP1 (secreted phosphoprotein 1), and LGALS3 (galectin 3) abundance correlated with disease activity and changed following treatment. Exploratory analyses suggested relationships between specific transcripts and individual histologic lesions, with the important representation of interferon-regulated genes. Our findings suggested that the current LN classification could be refined by the inclusion of molecular descriptors. Combining molecular and pathologic kidney biopsy phenotypes may hold promise to better classify disease and identify actionable treatment targets and merits further exploration in larger cohorts.

A comprehensive catalog of cancer driver mutations is essential for understanding tumorigenesis and developing therapies. Exome-sequencing studies have mapped many protein-coding drivers, yet few non-coding drivers are known because genome-wide discovery is challenging. We developed a driver discovery method, ActiveDriverWGS, and analyzed 120,788 cis-regulatory modules (CRMs) across 1,844 whole tumor genomes from the ICGC-TCGA PCAWG project. We found 30 CRMs with enriched SNVs and indels (FDR < 0.05). These frequently mutated regulatory elements (FMREs) were ubiquitously active in human tissues, showed long-range chromatin interactions and mRNA abundance associations with target genes, and were enriched in motif-rewiring mutations and structural variants. Genomic deletion of one FMRE in human cells caused proliferative deficiencies and transcriptional deregulation of cancer genes CCNB1IP1, CDH1, and CDKN2B, validating observations in FMRE-mutated tumors. Pathway analysis revealed further sub-significant FMREs at cancer genes and processes, indicating an unexplored landscape of infrequent driver mutations in the non-coding genome.

Many primary tumours have low levels of molecular oxygen (hypoxia), and hypoxic tumours respond poorly to therapy. Pan-cancer molecular hallmarks of tumour hypoxia remain poorly understood, with limited comprehension of its associations with specific mutational processes, non-coding driver genes and evolutionary features. Here, as part of the ICGC/TCGA Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium, which aggregated whole genome sequencing data from 2658 cancers across 38 tumour types, we quantify hypoxia in 1188 tumours spanning 27 cancer types. Elevated hypoxia associates with increased mutational load across cancer types, irrespective of underlying mutational class. The proportion of mutations attributed to several mutational signatures of unknown aetiology directly associates with the level of hypoxia, suggesting underlying mutational processes for these signatures. At the gene level, driver mutations in TP53, MYC and PTEN are enriched in hypoxic tumours, and mutations in PTEN interact with hypoxia to direct tumour evolutionary trajectories. Overall, hypoxia plays a critical role in shaping the genomic and evolutionary landscapes of cancer.

MYC is an oncogenic driver that regulates transcriptional activation and repression. Surprisingly, mechanisms by which MYC promotes malignant transformation remain unclear. We demonstrate that MYC interacts with the G9a H3K9-methyltransferase complex to control transcriptional repression. Inhibiting G9a hinders MYC chromatin binding at MYC-repressed genes and de-represses gene expression. By identifying the MYC box II region as essential for MYC-G9a interaction, a long-standing missing link between MYC transformation and gene repression is unveiled. Across breast cancer cell lines, the anti-proliferative response to G9a pharmacological inhibition correlates with MYC sensitivity and gene signatures. Consistently, genetically depleting G9a in vivo suppresses MYC-dependent tumor growth. These findings unveil G9a as an epigenetic regulator of MYC transcriptional repression and a therapeutic vulnerability in MYC-driven cancers.

IMA-08401 (C2) is a novel aryl hydrocarbon receptor (AHR) agonist and selective AHR modulator (SAHRM) that is structurally similar to laquinimod (LAQ). Both compounds are converted to the AHR-active metabolite DELAQ (IMA-06201) in vivo. SAHRMs have been proposed as therapeutic options for various autoimmune disorders. Clinical trials on LAQ have not reported any significant toxic outcomes and C2 has shown low toxicity in rats; however, their functional resemblance to the highly toxic AHR agonist 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) raises questions. Here, we characterize the hepatic transcriptomic changes induced by acute (single-dose) and subacute exposure (repeated dosing for 5 days followed by a 5-day recovery period) to C2 in Sprague-Dawley rats. Exposure to C2 leads to activation of the AHR, as shown by altered transcription of Cyp1a1. We identify a heightened response early after exposure that drops off by day 10. Acute exposure to C2 leads to changes to transcription of genes involved in antiviral and antibacterial responses, which highlights the immunomodulator effects of this AHR agonist. Subacute exposure causes an oxidative stress response in the liver, the consequences of which require further study on target tissues such as the CNS and immune system, both of which may be compromised in this patient population.

Driver gene mutations that are more prevalent in metastatic, castration-resistant prostate cancer (mCRPC) than localized disease represent candidate prognostic biomarkers. We analyze 1,844 localized (1,289) or mCRPC (555) tumors and quantify the prevalence of 113 somatic driver single nucleotide variants (SNVs), copy number aberrations (CNAs), and structural variants (SVs) in each state. One-third are significantly more prevalent in mCRPC than expected while a quarter are less prevalent. Mutations in AR and its enhancer are more prevalent in mCRPC, as are those in TP53, MYC, ZNRF3 and PRKDC. ZNRF3 loss is associated with decreased ZNRF3 mRNA abundance, WNT, cell cycle & PRC1/2 activity, and genomic instability. ZNRF3 loss, RNA downregulation and hypermethylation are prognostic of metastasis and overall survival, independent of clinical and pathologic indices. These data demonstrate a strategy for identifying biomarkers of localized cancer aggression, with ZNRF3 loss as a predictor of metastasis in prostate cancer.