IMPROVE-DD: Integrating multiple phenotype resources optimizes variant evaluation in genetically determined developmental disorders.

Diagnosing rare developmental disorders using genome-wide sequencing data commonly necessitates review of multiple plausible candidate variants, often using ontologies of categorical clinical terms. We show that Integrating Multiple Phenotype Resources Optimizes Variant Evaluation in Developmental Disorders (IMPROVE-DD) by incorporating additional classes of data commonly available to clinicians and recorded in health records. In doing so, we quantify the distinct contributions of sex, growth, and development in addition to Human Phenotype Ontology (HPO) terms and demonstrate added value from these readily available information sources. We use likelihood ratios for nominal and quantitative data and propose a classifier for HPO terms in this framework. This Bayesian framework results in more robust diagnoses. Using data systematically collected in the Deciphering Developmental Disorders study, we considered 77 genes with pathogenic/likely pathogenic variants in ≥10 individuals. All genes showed at least a satisfactory prediction by receiver operating characteristic when testing on training data (AUC ≥ 0.6), and HPO terms were the best predictor for the majority of genes, though a minority (13/77) of genes were better predicted by other phenotypic data types. Overall, classifiers based upon multiple integrated phenotypic data sources performed better than those based upon any individual source, and importantly, integrated models produced notably fewer false positives. Finally, we show that IMPROVE-DD models with good predictive performance on cross-validation can be constructed from relatively few individuals. This suggests new strategies for candidate gene prioritization and highlights the value of systematic clinical data collection to support diagnostic programs.

Multiomic Characterization of High-Grade Serous Ovarian Carcinoma Enables High-Resolution Patient Stratification.

PURPOSE: High-grade serous ovarian carcinoma (HGSOC) is the most common ovarian cancer type; most patients experience disease recurrence that accumulates chemoresistance, leading to treatment failure. Genomic and transcriptomic features have been associated with differential outcome and treatment response. However, the relationship between events at the gene sequence, copy number, and gene-expression levels remains poorly defined. EXPERIMENTAL DESIGN: We perform multiomic characterization of a large HGSOC cohort (n = 362) with detailed clinical annotation to interrogate the relationship between patient subgroups defined by specific molecular events. RESULTS: BRCA2-mutant (BRCA2m) and EMSY-overexpressing cases demonstrated prolonged survival [multivariable hazard ratios (HR) 0.40 and 0.51] and significantly higher first- and second-line chemotherapy response rate. CCNE1-gained (CCNE1g) cases demonstrated underrepresentation of FIGO stage IV cases, with shorter survival but no significant difference in treatment response. We demonstrate marked overlap between the TCGA- and Tothill-derived subtypes. IMR/C2 cases displayed higher BRCA1/2m frequency (25.5%, 32.5%) and significantly greater immune cell infiltration, whereas PRO/C5 cases had the highest CCNE1g rate (23.9%, 22.2%) and were uniformly low in immune cell infiltration. The survival benefit for cases with aberrations in homologous recombination repair (HRR) genes was apparent across all transcriptomic subtypes (HR range, 0.48-0.68). There was significant co-occurrence of RB loss and HRR gene aberrations; RB loss was further associated with favorable survival within HRR-aberrant cases (multivariable HR, 0.50). CONCLUSIONS: These data paint a high-resolution picture of the molecular landscape in HGSOC, better defining patients who may benefit most from specific molecular therapeutics and highlighting those for whom novel treatment strategies are needed to improve outcomes.

ISGylation drives basal breast tumour progression by promoting EGFR recycling and Akt signalling.

ISG15 is an ubiquitin-like modifier that is associated with reduced survival rates in breast cancer patients. The mechanism by which ISG15 achieves this however remains elusive. We demonstrate that modification of Rab GDP-Dissociation Inhibitor Beta (GDI2) by ISG15 (ISGylation) alters endocytic recycling of the EGF receptor (EGFR) in non-interferon stimulated cells using CRISPR-knock out models for ISGylation. By regulating EGFR trafficking, ISGylation enhances EGFR recycling and sustains Akt-signalling. We further show that Akt signalling positively correlates with levels of ISG15 and its E2-ligase in basal breast cancer cohorts, confirming the link between ISGylation and Akt signalling in human tumours. Persistent and enhanced Akt activation explains the more aggressive tumour behaviour observed in human breast cancers. We show that ISGylation can act as a driver of tumour progression rather than merely being a bystander.

Mutational bias in spermatogonia impacts the anatomy of regulatory sites in the human genome.

Mutation in the germline is the ultimate source of genetic variation, but little is known about the influence of germline chromatin structure on mutational processes. Using ATAC-seq, we profile the open chromatin landscape of human spermatogonia, the most proliferative cell type of the germline, identifying transcription factor binding sites (TFBSs) and PRDM9 binding sites, a subset of which will initiate meiotic recombination. We observe an increase in rare structural variant (SV) breakpoints at PRDM9-bound sites, implicating meiotic recombination in the generation of structural variation. Many germline TFBSs, such as NRF1, are also associated with increased rates of SV breakpoints, apparently independent of recombination. Singleton short insertions (≥5 bp) are highly enriched at TFBSs, particularly at sites bound by testis active TFs, and their rates correlate with those of structural variant breakpoints. Short insertions often duplicate the TFBS motif, leading to clustering of motif sites near regulatory regions in this male-driven evolutionary process. Increased mutation loads at germline TFBSs disproportionately affect neural enhancers with activity in spermatogonia, potentially altering neurodevelopmental regulatory architecture. Local chromatin structure in spermatogonia is thus pervasive in shaping both evolution and disease.

Integrated molecular characterisation of endometrioid ovarian carcinoma identifies opportunities for stratification.

Endometrioid ovarian carcinoma (EnOC) is an under-investigated ovarian cancer type. Recent studies have described disease subtypes defined by genomics and hormone receptor expression patterns; here, we determine the relationship between these subtyping layers to define the molecular landscape of EnOC with high granularity and identify therapeutic vulnerabilities in high-risk cases. Whole exome sequencing data were integrated with progesterone and oestrogen receptor (PR and ER) expression-defined subtypes in 90 EnOC cases following robust pathological assessment, revealing dominant clinical and molecular features in the resulting integrated subtypes. We demonstrate significant correlation between subtyping approaches: PR-high (PR + /ER + , PR + /ER-) cases were predominantly CTNNB1-mutant (73.2% vs 18.4%, P < 0.001), while PR-low (PR-/ER + , PR-/ER-) cases displayed higher TP53 mutation frequency (38.8% vs 7.3%, P = 0.001), greater genomic complexity (P = 0.007) and more frequent copy number alterations (P = 0.001). PR-high EnOC patients experience favourable disease-specific survival independent of clinicopathological and genomic features (HR = 0.16, 95% CI 0.04-0.71). TP53 mutation further delineates the outcome of patients with PR-low tumours (HR = 2.56, 95% CI 1.14-5.75). A simple, routinely applicable, classification algorithm utilising immunohistochemistry for PR and p53 recapitulated these subtypes and their survival profiles. The genomic profile of high-risk EnOC subtypes suggests that inhibitors of the MAPK and PI3K-AKT pathways, alongside PARP inhibitors, represent promising candidate agents for improving patient survival. Patients with PR-low TP53-mutant EnOC have the greatest unmet clinical need, while PR-high tumours-which are typically CTNNB1-mutant and TP53 wild-type-experience excellent survival and may represent candidates for trials investigating de-escalation of adjuvant chemotherapy to agents such as endocrine therapy.

Structural Variants at the BRCA1/2 Loci are a Common Source of Homologous Repair Deficiency in High-grade Serous Ovarian Carcinoma.

PURPOSE: The abundance and effects of structural variation at BRCA1/2 in tumors are not well understood. In particular, the impact of these events on homologous recombination repair deficiency (HRD) has yet to be demonstrated. EXPERIMENTAL DESIGN: Exploiting a large collection of whole-genome sequencing data from high-grade serous ovarian carcinoma (N = 205) together with matched RNA sequencing for the majority of tumors (N = 150), we have comprehensively characterized mutation and expression at BRCA1/2. RESULTS: In addition to the known spectrum of short somatic mutations (SSM), we discovered that multi-megabase structural variants (SV) were a frequent, unappreciated source of BRCA1/2 disruption in these tumors, and we found a genome-wide enrichment for large deletions at the BRCA1/2 loci across the cohort. These SVs independently affected a substantial proportion of patients (16%) in addition to those affected by SSMs (24%), conferring HRD and impacting patient survival. We also detail compound deficiencies involving SSMs and SVs at both loci, demonstrating that the strongest risk of HRD emerges from combined SVs at both BRCA1 and BRCA2 in the absence of SSMs. Furthermore, these SVs are abundant and disruptive in other cancer types. CONCLUSIONS: These results extend our understanding of the mutational landscape underlying HRD, increase the number of patients predicted to benefit from therapies exploiting HRD, and suggest there is currently untapped potential in SV detection for patient stratification.

Functional brain defects in a mouse model of a chromosomal t(1;11) translocation that disrupts DISC1 and confers increased risk of psychiatric illness.

A balanced t(1;11) translocation that directly disrupts DISC1 is linked to schizophrenia and affective disorders. We previously showed that a mutant mouse, named Der1, recapitulates the effect of the translocation upon DISC1 expression. Here, RNAseq analysis of Der1 mouse brain tissue found enrichment for dysregulation of the same genes and molecular pathways as in neuron cultures generated previously from human t(1;11) translocation carriers via the induced pluripotent stem cell route. DISC1 disruption therefore apparently accounts for a substantial proportion of the effects of the t(1;11) translocation. RNAseq and pathway analysis of the mutant mouse predicts multiple Der1-induced alterations converging upon synapse function and plasticity. Synaptosome proteomics confirmed that the Der1 mutation impacts synapse composition, and electrophysiology found reduced AMPA:NMDA ratio in hippocampal neurons, indicating changed excitatory signalling. Moreover, hippocampal parvalbumin-positive interneuron density is increased, suggesting that the Der1 mutation affects inhibitory control of neuronal circuits. These phenotypes predict that neurotransmission is impacted at many levels by DISC1 disruption in human t(1;11) translocation carriers. Notably, genes implicated in schizophrenia, depression and bipolar disorder by large-scale genetic studies are enriched among the Der1-dysregulated genes, just as we previously observed for the t(1;11) translocation carrier-derived neurons. Furthermore, RNAseq analysis predicts that the Der1 mutation primarily targets a subset of cell types, pyramidal neurons and interneurons, previously shown to be vulnerable to the effects of common schizophrenia-associated genetic variants. In conclusion, DISC1 disruption by the t(1;11) translocation may contribute to the psychiatric disorders of translocation carriers through commonly affected pathways and processes in neurotransmission.

Gene Co-Expression Network Analysis Identifies Vitamin D-Associated Gene Modules in Adult Normal Rectal Epithelium Following Supplementation.

Colorectal cancer (CRC) is a common, multifactorial disease. While observational studies have identified an association between lower vitamin D and higher CRC risk, supplementation trials have been inconclusive and the mechanisms by which vitamin D may modulate CRC risk are not well understood. We sought to perform a weighted gene co-expression network analysis (WGCNA) to identify modules present after vitamin D supplementation (when plasma vitamin D level was sufficient) which were absent before supplementation, and then to identify influential genes in those modules. The transcriptome from normal rectal mucosa biopsies of 49 individuals free from CRC were assessed before and after 12 weeks of 3200IU/day vitamin D (Fultium-D3) supplementation using paired-end total RNAseq. While the effects on expression patterns following vitamin D supplementation were subtle, WGCNA identified highly correlated genes forming gene modules. Four of the 17 modules identified in the post-vitamin D network were not preserved in the pre-vitamin D network, shedding new light on the biochemical impact of supplementation. These modules were enriched for GO terms related to the immune system, hormone metabolism, cell growth and RNA metabolism. Across the four treatment-associated modules, 51 hub genes were identified, with enrichment of 40 different transcription factor motifs in promoter regions of those genes, including VDR:RXR. Six of the hub genes were nominally differentially expressed in studies of vitamin D effects on adult normal mucosa organoids: LCN2, HLA-C, AIF1L, PTPRU, PDE4B and IFI6. By taking a gene-correlation network approach, we have described vitamin D induced changes to gene modules in normal human rectal epithelium in vivo, the target tissue from which CRC develops.

Molecular stratification of endometrioid ovarian carcinoma predicts clinical outcome.

Endometrioid ovarian carcinoma (EnOC) demonstrates substantial clinical and molecular heterogeneity. Here, we report whole exome sequencing of 112 EnOC cases following rigorous pathological assessment. We detect a high frequency of mutation in CTNNB1 (43%), PIK3CA (43%), ARID1A (36%), PTEN (29%), KRAS (26%), TP53 (26%) and SOX8 (19%), a recurrently-mutated gene previously unreported in EnOC. POLE and mismatch repair protein-encoding genes were mutated at lower frequency (6%, 18%) with significant co-occurrence. A molecular taxonomy is constructed, identifying clinically distinct EnOC subtypes: cases with TP53 mutation demonstrate greater genomic complexity, are commonly FIGO stage III/IV at diagnosis (48%), are frequently incompletely debulked (44%) and demonstrate inferior survival; conversely, cases with CTNNB1 mutation, which is mutually exclusive with TP53 mutation, demonstrate low genomic complexity and excellent clinical outcome, and are predominantly stage I/II at diagnosis (89%) and completely resected (87%). Moreover, we identify the WNT, MAPK/RAS and PI3K pathways as good candidate targets for molecular therapeutics in EnOC.

Pervasive lesion segregation shapes cancer genome evolution.

Cancers arise through the acquisition of oncogenic mutations and grow by clonal expansion1,2. Here we reveal that most mutagenic DNA lesions are not resolved into a mutated DNA base pair within a single cell cycle. Instead, DNA lesions segregate, unrepaired, into daughter cells for multiple cell generations, resulting in the chromosome-scale phasing of subsequent mutations. We characterize this process in mutagen-induced mouse liver tumours and show that DNA replication across persisting lesions can produce multiple alternative alleles in successive cell divisions, thereby generating both multiallelic and combinatorial genetic diversity. The phasing of lesions enables accurate measurement of strand-biased repair processes, quantification of oncogenic selection and fine mapping of sister-chromatid-exchange events. Finally, we demonstrate that lesion segregation is a unifying property of exogenous mutagens, including UV light and chemotherapy agents in human cells and tumours, which has profound implications for the evolution and adaptation of cancer genomes.

Increased ultra-rare variant load in an isolated Scottish population impacts exonic and regulatory regions.

Human population isolates provide a snapshot of the impact of historical demographic processes on population genetics. Such data facilitate studies of the functional impact of rare sequence variants on biomedical phenotypes, as strong genetic drift can result in higher frequencies of variants that are otherwise rare. We present the first whole genome sequencing (WGS) study of the VIKING cohort, a representative collection of samples from the isolated Shetland population in northern Scotland, and explore how its genetic characteristics compare to a mainland Scottish population. Our analyses reveal the strong contributions played by the founder effect and genetic drift in shaping genomic variation in the VIKING cohort. About one tenth of all high-quality variants discovered are unique to the VIKING cohort or are seen at frequencies at least ten fold higher than in more cosmopolitan control populations. Multiple lines of evidence also suggest relaxation of purifying selection during the evolutionary history of the Shetland isolate. We demonstrate enrichment of ultra-rare VIKING variants in exonic regions and for the first time we also show that ultra-rare variants are enriched within regulatory regions, particularly promoters, suggesting that gene expression patterns may diverge relatively rapidly in human isolates.

Zebrafish MITF-Low Melanoma Subtype Models Reveal Transcriptional Subclusters and MITF-Independent Residual Disease.

The melanocyte-inducing transcription factor (MITF)-low melanoma transcriptional signature is predictive of poor outcomes for patients, but little is known about its biological significance, and animal models are lacking. Here, we used zebrafish genetic models with low activity of Mitfa (MITF-low) and established that the MITF-low state is causal of melanoma progression and a predictor of melanoma biological subtype. MITF-low zebrafish melanomas resembled human MITF-low melanomas and were enriched for stem and invasive (mesenchymal) gene signatures. MITF-low activity coupled with a p53 mutation was sufficient to promote superficial growth melanomas, whereas BRAFV600E accelerated MITF-low melanoma onset and further promoted the development of MITF-high nodular growth melanomas. Genetic inhibition of MITF activity led to rapid regression; recurrence occurred following reactivation of MITF. At the regression site, there was minimal residual disease that was resistant to loss of MITF activity (termed MITF-independent cells) with very low-to-no MITF activity or protein. Transcriptomic analysis of MITF-independent residual disease showed enrichment of mesenchymal and neural crest stem cell signatures similar to human therapy-resistant melanomas. Single-cell RNA sequencing revealed MITF-independent residual disease was heterogeneous depending on melanoma subtype. Further, there was a shared subpopulation of residual disease cells that was enriched for a neural crest G0-like state that preexisted in the primary tumor and remained present in recurring melanomas. These findings suggest that invasive and stem-like programs coupled with cellular heterogeneity contribute to poor outcomes for MITF-low melanoma patients and that MITF-independent subpopulations are an important therapeutic target to achieve long-term survival outcomes. SIGNIFICANCE: This study provides a useful model for MITF-low melanomas and MITF-independent cell populations that can be used to study the mechanisms that drive these tumors as well as identify potential therapeutic options.Graphical Abstract: http://cancerres.aacrjournals.org/content/canres/79/22/5769/F1.large.jpg.

An actionable KCNH2 Long QT Syndrome variant detected by sequence and haplotype analysis in a population research cohort.

The Viking Health Study Shetland is a population-based research cohort of 2,122 volunteer participants with ancestry from the Shetland Isles in northern Scotland. The high kinship and detailed phenotype data support a range of approaches for associating rare genetic variants, enriched in this isolate population, with quantitative traits and diseases. As an exemplar, the c.1750G > A; p.Gly584Ser variant within the coding sequence of the KCNH2 gene implicated in Long QT Syndrome (LQTS), which occurred once in 500 whole genome sequences from this population, was investigated. Targeted sequencing of the KCNH2 gene in family members of the initial participant confirmed the presence of the sequence variant and identified two further members of the same family pedigree who shared the variant. Investigation of these three related participants for whom single nucleotide polymorphism (SNP) array genotypes were available allowed a unique shared haplotype of 1.22 Mb to be defined around this locus. Searching across the full cohort for this haplotype uncovered two additional apparently unrelated individuals with no known genealogical connection to the original kindred. All five participants with the defined haplotype were shown to share the rare variant by targeted Sanger sequencing. If this result were verified in a healthcare setting, it would be considered clinically actionable, and has been actioned in relatives ascertained independently through clinical presentation. The General Practitioners of four study participants with the rare variant were alerted to the research findings by letters outlining the phenotype (prolonged electrocardiographic QTc interval). A lack of detectable haplotype sharing between c.1750G > A; p.Gly584Ser chromosomes from previously reported individuals from Finland and those in this study from Shetland suggests that this mutation has arisen more than once in human history. This study showcases the potential value of isolate population-based research resources for genomic medicine. It also illustrates some challenges around communication of actionable findings in research participants in this context.

Clinical and molecular characterization of ovarian carcinoma displaying isolated lymph node relapse.

BACKGROUND: Disease relapse is the primary cause of death from ovarian carcinoma. Isolated lymph node relapse is a rare pattern of ovarian carcinoma recurrence, with a reported median postrelapse survival of 2.5 to 4 years. To date, investigations have not compared isolated lymph node relapse ovarian carcinoma directly to a matched extranodal relapse cohort or performed molecular characterization of cases that subsequently experience isolated lymph node relapse. OBJECTIVE: Here we seek to compare the clinical outcome, tumor-infiltrating lymphocyte burden, and frequency of known prognostic genomic events in isolated lymph node relapse ovarian carcinoma vs extranodal relapse ovarian carcinoma. STUDY DESIGN: Forty-nine isolated lymph node relapse ovarian carcinoma patients were identified and matched to 49 extranodal relapse cases using the Edinburgh Ovarian Cancer Database, from which the clinical data for identified patients were retrieved. Matching criteria were disease stage, histologic subtype and grade, extent of residual disease following surgical debulking, and age at diagnosis. Clinicopathologic factors and survival data were compared between the isolated lymph node relapse and extranodal relapse cohorts. Genomic characterization of tumor material from diagnosis was performed using panel-based high-throughput sequencing and tumor-infiltrating T cell burden was assessed using immunohistochemistry for CD3+ and CD8+ cells. RESULTS: Isolated lymph node relapse cases demonstrated significantly prolonged postrelapse survival and overall survival vs extranodal relapse upon multivariable analysis (HRmulti = 0.52 [0.33-0.84] and 0.51 [0.31-0.84]). Diagnostic specimens from high-grade serous ovarian carcinomas that subsequently displayed isolated lymph node relapse harbored significantly greater CD3+ and CD8+ cell infiltration compared to extranodal relapse cases (P = .001 and P = .009, Bonferroni-adjusted P = .003 and P = .019). Isolated lymph node relapse high-grade serous ovarian carcinoma cases did not show marked enrichment or depletion of cases with BRCA1/2 mutation or CCNE1 copy number gain when compared to their extranodal relapse counterparts (24.4% vs 19.4% and 18.2% vs 22.6%, P = .865 and P = .900). CONCLUSION: Isolated lymph node relapse ovarian carcinoma represents a distinct clinical entity with favorable outcome compared to extranodal relapse. There was no clear enrichment or depletion of BRCA1/2 mutation or CCNE1 gain in the isolated lymph node relapse ovarian carcinoma cohort compared with extranodal relapse cases, suggesting that these known prognostic genomically defined subtypes of disease do not display markedly altered propensity for isolated lymph node relapse. Diagnostic tumor material from isolated lymph node relapse patients demonstrated greater CD3+ and CD8+ cell infiltration, indicating stronger tumor engagement by T cell populations, which may contribute to the more indolent disease course of isolated lymph node relapse.

Modeling double strand break susceptibility to interrogate structural variation in cancer.

BACKGROUND: Structural variants (SVs) are known to play important roles in a variety of cancers, but their origins and functional consequences are still poorly understood. Many SVs are thought to emerge from errors in the repair processes following DNA double strand breaks (DSBs). RESULTS: We used experimentally quantified DSB frequencies in cell lines with matched chromatin and sequence features to derive the first quantitative genome-wide models of DSB susceptibility. These models are accurate and provide novel insights into the mutational mechanisms generating DSBs. Models trained in one cell type can be successfully applied to others, but a substantial proportion of DSBs appear to reflect cell type-specific processes. Using model predictions as a proxy for susceptibility to DSBs in tumors, many SV-enriched regions appear to be poorly explained by selectively neutral mutational bias alone. A substantial number of these regions show unexpectedly high SV breakpoint frequencies given their predicted susceptibility to mutation and are therefore credible targets of positive selection in tumors. These putatively positively selected SV hotspots are enriched for genes previously shown to be oncogenic. In contrast, several hundred regions across the genome show unexpectedly low levels of SVs, given their relatively high susceptibility to mutation. These novel coldspot regions appear to be subject to purifying selection in tumors and are enriched for active promoters and enhancers. CONCLUSIONS: We conclude that models of DSB susceptibility offer a rigorous approach to the inference of SVs putatively subject to selection in tumors.

Breaking point: the genesis and impact of structural variation in tumours.

Somatic structural variants undoubtedly play important roles in driving tumourigenesis. This is evident despite the substantial technical challenges that remain in accurately detecting structural variants and their breakpoints in tumours and in spite of our incomplete understanding of the impact of structural variants on cellular function. Developments in these areas of research contribute to the ongoing discovery of structural variation with a clear impact on the evolution of the tumour and on the clinical importance to the patient. Recent large whole genome sequencing studies have reinforced our impression of each tumour as a unique combination of mutations but paradoxically have also discovered similar genome-wide patterns of single-nucleotide and structural variation between tumours. Statistical methods have been developed to deconvolute mutation patterns, or signatures, that recur across samples, providing information about the mutagens and repair processes that may be active in a given tumour. These signatures can guide treatment by, for example, highlighting vulnerabilities in a particular tumour to a particular chemotherapy. Thus, although the complete reconstruction of the full evolutionary trajectory of a tumour genome remains currently out of reach, valuable data are already emerging to improve the treatment of cancer.

DISC1 regulates N-methyl-D-aspartate receptor dynamics: abnormalities induced by a Disc1 mutation modelling a translocation linked to major mental illness.

The neuromodulatory gene DISC1 is disrupted by a t(1;11) translocation that is highly penetrant for schizophrenia and affective disorders, but how this translocation affects DISC1 function is incompletely understood. N-methyl-D-aspartate receptors (NMDAR) play a central role in synaptic plasticity and cognition, and are implicated in the pathophysiology of schizophrenia through genetic and functional studies. We show that the NMDAR subunit GluN2B complexes with DISC1-associated trafficking factor TRAK1, while DISC1 interacts with the GluN1 subunit and regulates dendritic NMDAR motility in cultured mouse neurons. Moreover, in the first mutant mouse that models DISC1 disruption by the translocation, the pool of NMDAR transport vesicles and surface/synaptic NMDAR expression are increased. Since NMDAR cell surface/synaptic expression is tightly regulated to ensure correct function, these changes in the mutant mouse are likely to affect NMDAR signalling and synaptic plasticity. Consistent with these observations, RNASeq analysis of the translocation carrier-derived human neurons indicates abnormalities of excitatory synapses and vesicle dynamics. RNASeq analysis of the human neurons also identifies many differentially expressed genes previously highlighted as putative schizophrenia and/or depression risk factors through large-scale genome-wide association and copy number variant studies, indicating that the translocation triggers common disease pathways that are shared with unrelated psychiatric patients. Altogether, our findings suggest that translocation-induced disease mechanisms are likely to be relevant to mental illness in general, and that such disease mechanisms include altered NMDAR dynamics and excitatory synapse function. This could contribute to the cognitive disorders displayed by translocation carriers.

Chromatin loop anchors are associated with genome instability in cancer and recombination hotspots in the germline.

BACKGROUND: Chromatin loops form a basic unit of interphase nuclear organization, with chromatin loop anchor points providing contacts between regulatory regions and promoters. However, the mutational landscape at these anchor points remains under-studied. Here, we describe the unusual patterns of somatic mutations and germline variation associated with loop anchor points and explore the underlying features influencing these patterns. RESULTS: Analyses of whole genome sequencing datasets reveal that anchor points are strongly depleted for single nucleotide variants (SNVs) in tumours. Despite low SNV rates in their genomic neighbourhood, anchor points emerge as sites of evolutionary innovation, showing enrichment for structural variant (SV) breakpoints and a peak of SNVs at focal CTCF sites within the anchor points. Both CTCF-bound and non-CTCF anchor points harbour an excess of SV breakpoints in multiple tumour types and are prone to double-strand breaks in cell lines. Common fragile sites, which are hotspots for genome instability, also show elevated numbers of intersecting loop anchor points. Recurrently disrupted anchor points are enriched for genes with functions in cell cycle transitions and regions associated with predisposition to cancer. We also discover a novel class of CTCF-bound anchor points which overlap meiotic recombination hotspots and are enriched for the core PRDM9 binding motif, suggesting that the anchor points have been foci for diversity generated during recent human evolution. CONCLUSIONS: We suggest that the unusual chromatin environment at loop anchor points underlies the elevated rates of variation observed, marking them as sites of regulatory importance but also genomic fragility.

Conserved temporal ordering of promoter activation implicates common mechanisms governing the immediate early response across cell types and stimuli.

The promoters of immediate early genes (IEGs) are rapidly activated in response to an external stimulus. These genes, also known as primary response genes, have been identified in a range of cell types, under diverse extracellular signals and using varying experimental protocols. Whereas genomic dissection on a case-by-case basis has not resulted in a comprehensive catalogue of IEGs, a rigorous meta-analysis of eight genome-wide FANTOM5 CAGE (cap analysis of gene expression) time course datasets reveals successive waves of promoter activation in IEGs, recapitulating known relationships between cell types and stimuli: we obtain a set of 57 (42 protein-coding) candidate IEGs possessing promoters that consistently drive a rapid but transient increase in expression over time. These genes show significant enrichment for known IEGs reported previously, pathways associated with the immediate early response, and include a number of non-coding RNAs with roles in proliferation and differentiation. Surprisingly, we also find strong conservation of the ordering of activation for these genes, such that 77 pairwise promoter activation orderings are conserved. Using the leverage of comprehensive CAGE time series data across cell types, we also document the extensive alternative promoter usage by such genes, which is likely to have been a barrier to their discovery until now. The common activation ordering of the core set of early-responding genes we identify may indicate conserved underlying regulatory mechanisms. By contrast, the considerably larger number of transiently activated genes that are specific to each cell type and stimulus illustrates the breadth of the primary response.

Enhanced response rate to pegylated liposomal doxorubicin in high grade serous ovarian carcinomas harbouring BRCA1 and BRCA2 aberrations.

BACKGROUND: Approximately 10-15% of ovarian carcinomas (OC) are attributed to inherited susceptibility, the majority of which are due to mutations in BRCA1 or BRCA2 (BRCA1/2). These patients display superior clinical outcome, including enhanced sensitivity to platinum-based chemotherapy. Here, we seek to investigate whether BRCA1/2 status influences the response rate to single-agent pegylated liposomal doxorubicin (PLD) in high grade serous (HGS) OC. METHODS: One hundred and forty-eight patients treated with single-agent PLD were identified retrospectively from the Edinburgh Ovarian Cancer Database. DNA was extracted from formalin-fixed paraffin-embedded (FFPE) archival tumour material and sequenced using the Ion Ampliseq BRCA1 and BRCA2 panel. A minimum variant allele frequency threshold was applied to correct for sequencing artefacts associated with formalin fixation. RESULTS: A superior response rate to PLD was observed in patients with HGS OC who harboured variants likely to affect BRCA1 or BRCA2 function compared to the BRCA1/2 wild-type population (36%, 9 of 25 patients versus 12.1%, 7 of 58 patients; p = 0.016). An enhanced response rate was also seen in patients harbouring only the BRCA1 SNP rs1799950, predicted to be detrimental to BRCA1 function (50%, 3 of 6 patients versus 12.1%, 7 of 58 patients; p = 0.044). CONCLUSIONS: These data demonstrate that HGS OC patients with BRCA1/2 variants predicted damaging to protein function experience superior sensitivity to PLD, consistent with impaired DNA repair. Further characterisation of rs1799950 is now warranted in relation to chemosensitivity and susceptibility to developing ovarian carcinoma.