Single-cell genomic variation induced by mutational processes in cancer.

How cell-to-cell copy number alterations that underpin genomic instability1 in human cancers drive genomic and phenotypic variation, and consequently the evolution of cancer2, remains understudied. Here, by applying scaled single-cell whole-genome sequencing3 to wild-type, TP53-deficient and TP53-deficient;BRCA1-deficient or TP53-deficient;BRCA2-deficient mammary epithelial cells (13,818 genomes), and to primary triple-negative breast cancer (TNBC) and high-grade serous ovarian cancer (HGSC) cells (22,057 genomes), we identify three distinct 'foreground' mutational patterns that are defined by cell-to-cell structural variation. Cell- and clone-specific high-level amplifications, parallel haplotype-specific copy number alterations and copy number segment length variation (serrate structural variations) had measurable phenotypic and evolutionary consequences. In TNBC and HGSC, clone-specific high-level amplifications in known oncogenes were highly prevalent in tumours bearing fold-back inversions, relative to tumours with homologous recombination deficiency, and were associated with increased clone-to-clone phenotypic variation. Parallel haplotype-specific alterations were also commonly observed, leading to phylogenetic evolutionary diversity and clone-specific mono-allelic expression. Serrate variants were increased in tumours with fold-back inversions and were highly correlated with increased genomic diversity of cellular populations. Together, our findings show that cell-to-cell structural variation contributes to the origins of phenotypic and evolutionary diversity in TNBC and HGSC, and provide insight into the genomic and mutational states of individual cancer cells.

Accurate determination of CRISPR-mediated gene fitness in transplantable tumours.

Assessing tumour gene fitness in physiologically-relevant model systems is challenging due to biological features of in vivo tumour regeneration, including extreme variations in single cell lineage progeny. Here we develop a reproducible, quantitative approach to pooled genetic perturbation in patient-derived xenografts (PDXs), by encoding single cell output from transplanted CRISPR-transduced cells in combination with a Bayesian hierarchical model. We apply this to 181 PDX transplants from 21 breast cancer patients. We show that uncertainty in fitness estimates depends critically on the number of transplant cell clones and the variability in clone sizes. We use a pathway-directed allelic series to characterize Notch signaling, and quantify TP53 / MDM2 drug-gene conditional fitness in outlier patients. We show that fitness outlier identification can be mirrored by pharmacological perturbation. Overall, we demonstrate that the gene fitness landscape in breast PDXs is dominated by inter-patient differences.

Ubiquitin-mediated DNA damage response is synthetic lethal with G-quadruplex stabilizer CX-5461.

CX-5461 is a G-quadruplex (G4) ligand currently in trials with initial indications of clinical activity in cancers with defects in homologous recombination repair. To identify more genetic defects that could sensitize tumors to CX-5461, we tested synthetic lethality for 480 DNA repair and genome maintenance genes to CX-5461, pyridostatin (PDS), a structurally unrelated G4-specific stabilizer, and BMH-21, which binds GC-rich DNA but not G4 structures. We identified multiple members of HRD, Fanconi Anemia pathways, and POLQ, a polymerase with a helicase domain important for G4 structure resolution. Significant synthetic lethality was observed with UBE2N and RNF168, key members of the DNA damage response associated ubiquitin signaling pathway. Loss-of-function of RNF168 and UBE2N resulted in significantly lower cell survival in the presence of CX-5461 and PDS but not BMH-21. RNF168 recruitment and histone ubiquitination increased with CX-5461 treatment, and nuclear ubiquitination response frequently co-localized with G4 structures. Pharmacological inhibition of UBE2N acted synergistically with CX-5461. In conclusion, we have uncovered novel genetic vulnerabilities to CX-5461 with potential significance for patient selection in future clinical trials.

Age-correlated protein and transcript expression in breast cancer and normal breast tissues is dominated by host endocrine effects.

The magnitude and scope of intrinsic age-correlated and host endocrine age-correlated gene expression in breast cancer is not well understood. From age-correlated gene expression in 3,071 breast cancer transcriptomes and epithelial protein expression of 42 markers in 5,001 breast cancers and 537 normal breast tissues, we identified a majority of age-correlated genes as putatively regulated by age-dependent estrogen signaling. Surprisingly, these included genes encoding the chromatin modifier EZH2 (which had a negative age correlation) and associated H3K27me3 (which had an inverse, positive age correlation). Among The Cancer Genome Atlas lung, thyroid, kidney and prostate transcriptomes, the largest overlap with breast cancer in age-correlated transcripts was lung cancer, for which about one-third of overlapping age-correlated transcripts appeared to be estrogen regulated. Age-quartile-stratified outcomes analysis of 3,500 breast cancers using EZH2, H3K27me3, FOXA1 and BCL2 proteins revealed distinct age-related prognostic significance. Age correlation in gene expression may thus be an important factor in ER, EZH2, H3K27me3 and other biomarker assessment and treatment strategies.

Dynamics of breast-cancer relapse reveal late-recurring ER-positive genomic subgroups.

The rates and routes of lethal systemic spread in breast cancer are poorly understood owing to a lack of molecularly characterized patient cohorts with long-term, detailed follow-up data. Long-term follow-up is especially important for those with oestrogen-receptor (ER)-positive breast cancers, which can recur up to two decades after initial diagnosis1-6. It is therefore essential to identify patients who have a high risk of late relapse7-9. Here we present a statistical framework that models distinct disease stages (locoregional recurrence, distant recurrence, breast-cancer-related death and death from other causes) and competing risks of mortality from breast cancer, while yielding individual risk-of-recurrence predictions. We apply this model to 3,240 patients with breast cancer, including 1,980 for whom molecular data are available, and delineate spatiotemporal patterns of relapse across different categories of molecular information (namely immunohistochemical subtypes; PAM50 subtypes, which are based on gene-expression patterns10,11; and integrative or IntClust subtypes, which are based on patterns of genomic copy-number alterations and gene expression12,13). We identify four late-recurring integrative subtypes, comprising about one quarter (26%) of tumours that are both positive for ER and negative for human epidermal growth factor receptor 2, each with characteristic tumour-driving alterations in genomic copy number and a high risk of recurrence (mean 47-62%) up to 20 years after diagnosis. We also define a subgroup of triple-negative breast cancers in which cancer rarely recurs after five years, and a separate subgroup in which patients remain at risk. Use of the integrative subtypes improves the prediction of late, distant relapse beyond what is possible with clinical covariates (nodal status, tumour size, tumour grade and immunohistochemical subtype). These findings highlight opportunities for improved patient stratification and biomarker-driven clinical trials.

Integrated structural variation and point mutation signatures in cancer genomes using correlated topic models.

Mutation signatures in cancer genomes reflect endogenous and exogenous mutational processes, offering insights into tumour etiology, features for prognostic and biologic stratification and vulnerabilities to be exploited therapeutically. We present a novel machine learning formalism for improved signature inference, based on multi-modal correlated topic models (MMCTM) which can at once infer signatures from both single nucleotide and structural variation counts derived from cancer genome sequencing data. We exemplify the utility of our approach on two hormone driven, DNA repair deficient cancers: breast and ovary (n = 755 samples total). We show how introducing correlated structure both within and between modes of mutation can increase accuracy of signature discovery, particularly in the context of sparse data. Our study emphasizes the importance of integrating multiple mutation modes for signature discovery and patient stratification, and provides a statistical modeling framework to incorporate additional features of interest for future studies.

CX-5461 is a DNA G-quadruplex stabilizer with selective lethality in BRCA1/2 deficient tumours.

G-quadruplex DNAs form four-stranded helical structures and are proposed to play key roles in different cellular processes. Targeting G-quadruplex DNAs for cancer treatment is a very promising prospect. Here, we show that CX-5461 is a G-quadruplex stabilizer, with specific toxicity against BRCA deficiencies in cancer cells and polyclonal patient-derived xenograft models, including tumours resistant to PARP inhibition. Exposure to CX-5461, and its related drug CX-3543, blocks replication forks and induces ssDNA gaps or breaks. The BRCA and NHEJ pathways are required for the repair of CX-5461 and CX-3543-induced DNA damage and failure to do so leads to lethality. These data strengthen the concept of G4 targeting as a therapeutic approach, specifically for targeting HR and NHEJ deficient cancers and other tumours deficient for DNA damage repair. CX-5461 is now in advanced phase I clinical trial for patients with BRCA1/2 deficient tumours (Canadian trial, NCT02719977, opened May 2016).

CLK-dependent exon recognition and conjoined gene formation revealed with a novel small molecule inhibitor.

CDC-like kinase phosphorylation of serine/arginine-rich proteins is central to RNA splicing reactions. Yet, the genomic network of CDC-like kinase-dependent RNA processing events remains poorly defined. Here, we explore the connectivity of genomic CDC-like kinase splicing functions by applying graduated, short-exposure, pharmacological CDC-like kinase inhibition using a novel small molecule (T3) with very high potency, selectivity, and cell-based stability. Using RNA-Seq, we define CDC-like kinase-responsive alternative splicing events, the large majority of which monotonically increase or decrease with increasing CDC-like kinase inhibition. We show that distinct RNA-binding motifs are associated with T3 response in skipped exons. Unexpectedly, we observe dose-dependent conjoined gene transcription, which is associated with motif enrichment in the last and second exons of upstream and downstream partners, respectively. siRNA knockdown of CLK2-associated genes significantly increases conjoined gene formation. Collectively, our results reveal an unexpected role for CDC-like kinase in conjoined gene formation, via regulation of 3'-end processing and associated splicing factors.The phosphorylation of serine/arginine-rich proteins by CDC-like kinase is a central regulatory mechanism for RNA splicing reactions. Here, the authors synthesize a novel small molecule CLK inhibitor and map CLK-responsive alternative splicing events and discover an effect on conjoined gene transcription.

Combined Use of Gene Expression Modeling and siRNA Screening Identifies Genes and Pathways Which Enhance the Activity of Cisplatin When Added at No Effect Levels to Non-Small Cell Lung Cancer Cells In Vitro.

Platinum-based combination chemotherapy is the standard treatment for advanced non-small cell lung cancer (NSCLC). While cisplatin is effective, its use is not curative and resistance often emerges. As a consequence of microenvironmental heterogeneity, many tumour cells are exposed to sub-lethal doses of cisplatin. Further, genomic heterogeneity and unique tumor cell sub-populations with reduced sensitivities to cisplatin play a role in its effectiveness within a site of tumor growth. Being exposed to sub-lethal doses will induce changes in gene expression that contribute to the tumour cell's ability to survive and eventually contribute to the selective pressures leading to cisplatin resistance. Such changes in gene expression, therefore, may contribute to cytoprotective mechanisms. Here, we report on studies designed to uncover how tumour cells respond to sub-lethal doses of cisplatin. A microarray study revealed changes in gene expressions that occurred when A549 cells were exposed to a no-observed-effect level (NOEL) of cisplatin (e.g. the IC10). These data were integrated with results from a genome-wide siRNA screen looking for novel therapeutic targets that when inhibited transformed a NOEL of cisplatin into one that induced significant increases in lethality. Pathway analyses were performed to identify pathways that could be targeted to enhance cisplatin activity. We found that over 100 genes were differentially expressed when A549 cells were exposed to a NOEL of cisplatin. Pathways associated with apoptosis and DNA repair were activated. The siRNA screen revealed the importance of the hedgehog, cell cycle regulation, and insulin action pathways in A549 cell survival and response to cisplatin treatment. Results from both datasets suggest that RRM2B, CABYR, ALDH3A1, and FHL2 could be further explored as cisplatin-enhancing gene targets. Finally, pathways involved in repairing double-strand DNA breaks and INO80 chromatin remodeling were enriched in both datasets, warranting further research into combinations of cisplatin and therapeutics targeting these pathways.

Translational Activation of HIF1α by YB-1 Promotes Sarcoma Metastasis.

Metastatic dissemination is the leading cause of death in cancer patients, which is particularly evident for high-risk sarcomas such as Ewing sarcoma, osteosarcoma, and rhabdomyosarcoma. Previous research identified a crucial role for YB-1 in the epithelial-to-mesenchymal transition (EMT) and metastasis of epithelial malignancies. Based on clinical data and two distinct animal models, we now report that YB-1 is also a major metastatic driver in high-risk sarcomas. Our data establish YB-1 as a critical regulator of hypoxia-inducible factor 1α (HIF1α) expression in sarcoma cells. YB-1 enhances HIF1α protein expression by directly binding to and activating translation of HIF1A messages. This leads to HIF1α-mediated sarcoma cell invasion and enhanced metastatic capacity in vivo, highlighting a translationally regulated YB-1-HIF1α axis in sarcoma metastasis.

A co-culture genome-wide RNAi screen with mammary epithelial cells reveals transmembrane signals required for growth and differentiation.

INTRODUCTION: The extracellular signals regulating mammary epithelial cell growth are of relevance to understanding the pathophysiology of mammary epithelia, yet they remain poorly characterized. In this study, we applied an unbiased approach to understanding the functional role of signalling molecules in several models of normal physiological growth and translated these results to the biological understanding of breast cancer subtypes. METHODS: We developed and utilized a cytogenetically normal clonal line of hTERT immortalized human mammary epithelial cells in a fibroblast-enhanced co-culture assay to conduct a genome-wide small interfering RNA (siRNA) screen for evaluation of the functional effect of silencing each gene. Our selected endpoint was inhibition of growth. In rigorous postscreen validation processes, including quantitative RT-PCR, to ensure on-target silencing, deconvolution of pooled siRNAs and independent confirmation of effects with lentiviral short-hairpin RNA constructs, we identified a subset of genes required for mammary epithelial cell growth. Using three-dimensional Matrigel growth and differentiation assays and primary human mammary epithelial cell colony assays, we confirmed that these growth effects were not limited to the 184-hTERT cell line. We utilized the METABRIC dataset of 1,998 breast cancer patients to evaluate both the differential expression of these genes across breast cancer subtypes and their prognostic significance. RESULTS: We identified 47 genes that are critically important for fibroblast-enhanced mammary epithelial cell growth. This group was enriched for several axonal guidance molecules and G protein-coupled receptors, as well as for the endothelin receptor PROCR. The majority of genes (43 of 47) identified in two dimensions were also required for three-dimensional growth, with HSD17B2, SNN and PROCR showing greater than tenfold reductions in acinar formation. Several genes, including PROCR and the neuronal pathfinding molecules EFNA4 and NTN1, were also required for proper differentiation and polarization in three-dimensional cultures. The 47 genes identified showed a significant nonrandom enrichment for differential expression among 10 molecular subtypes of breast cancer sampled from 1,998 patients. CD79A, SERPINH1, KCNJ5 and TMEM14C exhibited breast cancer subtype-independent overall survival differences. CONCLUSION: Diverse transmembrane signals are required for mammary epithelial cell growth in two-dimensional and three-dimensional conditions. Strikingly, we define novel roles for axonal pathfinding receptors and ligands and the endothelin receptor in both growth and differentiation.

A tumor DNA complex aberration index is an independent predictor of survival in breast and ovarian cancer.

Complex focal chromosomal rearrangements in cancer genomes, also called "firestorms", can be scored from DNA copy number data. The complex arm-wise aberration index (CAAI) is a score that captures DNA copy number alterations that appear as focal complex events in tumors, and has potential prognostic value in breast cancer. This study aimed to validate this DNA-based prognostic index in breast cancer and test for the first time its potential prognostic value in ovarian cancer. Copy number alteration (CNA) data from 1950 breast carcinomas (METABRIC cohort) and 508 high-grade serous ovarian carcinomas (TCGA dataset) were analyzed. Cases were classified as CAAI positive if at least one complex focal event was scored. Complex alterations were frequently localized on chromosome 8p (n = 159), 17q (n = 176) and 11q (n = 251). CAAI events on 11q were most frequent in estrogen receptor positive (ER+) cases and on 17q in estrogen receptor negative (ER-) cases. We found only a modest correlation between CAAI and the overall rate of genomic instability (GII) and number of breakpoints (r = 0.27 and r = 0.42, p < 0.001). Breast cancer specific survival (BCSS), overall survival (OS) and ovarian cancer progression free survival (PFS) were used as clinical end points in Cox proportional hazard model survival analyses. CAAI positive breast cancers (43%) had higher mortality: hazard ratio (HR) of 1.94 (95%CI, 1.62-2.32) for BCSS, and of 1.49 (95%CI, 1.30-1.71) for OS. Representations of the 70-gene and the 21-gene predictors were compared with CAAI in multivariable models and CAAI was independently significant with a Cox adjusted HR of 1.56 (95%CI, 1.23-1.99) for ER+ and 1.55 (95%CI, 1.11-2.18) for ER- disease. None of the expression-based predictors were prognostic in the ER- subset. We found that a model including CAAI and the two expression-based prognostic signatures outperformed a model including the 21-gene and 70-gene signatures but excluding CAAI. Inclusion of CAAI in the clinical prognostication tool PREDICT significantly improved its performance. CAAI positive ovarian cancers (52%) also had worse prognosis: HRs of 1.3 (95%CI, 1.1-1.7) for PFS and 1.3 (95%CI, 1.1-1.6) for OS. This study validates CAAI as an independent predictor of survival in both ER+ and ER- breast cancer and reveals a significant prognostic value for CAAI in high-grade serous ovarian cancer.

Up-regulation of the interferon-related genes in BRCA2 knockout epithelial cells.

BRCA2 mutations are significantly associated with early-onset breast cancer, and the tumour-suppressing function of BRCA2 has been attributed to its involvement in homologous recombination (HR)-mediated DNA repair. In order to identify additional functions of BRCA2, we generated BRCA2-knockout HCT116 human colorectal carcinoma cells. Using genome-wide microarray analyses, we have discovered a link between the loss of BRCA2 and the up-regulation of a subset of interferon (IFN)-related genes, including APOBEC3F and APOBEC3G. The over-expression of IFN-related genes was confirmed in different human BRCA2(-/-) and mouse Brca2(-/-) tumour cell lines, and was independent of senescence and apoptosis. In isogenic wild-type BRCA2 cells, we observed over-expression of IFN-related genes after treatment with DNA-damaging agents, and following ionizing radiation. Cells with endogenous DNA damage because of defective BRCA1 or RAD51 also exhibited over-expression of IFN-related genes. Transcriptional activity of the IFN-stimulated response element (ISRE) was increased in BRCA2 knockout cells, and the expression of BRCA2 greatly decreased IFNα-stimulated ISRE reporter activity, suggesting that BRCA2 directly represses the expression of IFN-related genes through the ISRE. Finally, the colony-forming capacity of BRCA2 knockout cells was significantly reduced in the presence of either IFNß or IFNγ, suggesting that IFNs may have potential as therapeutic agents in cancer cells with BRCA2 mutations. The GEO Accession No. for microarray analysis is GSE54830.

TP53 mutation spectrum in breast cancer is subtype specific and has distinct prognostic relevance.

PURPOSE: In breast cancer, the TP53 gene is frequently mutated and the mutations have been associated with poor prognosis. The prognostic impact of the different types of TP53 mutations across the different molecular subtypes is still poorly understood. Here, we characterize the spectrum and prognostic significance of TP53 mutations with respect to the PAM50 subtypes and integrative clusters (IC). EXPERIMENTAL DESIGN: TP53 mutation status was obtained for 1,420 tumor samples from the METABRIC cohort by sequencing all coding exons using the Sanger method. RESULTS: TP53 mutations were found in 28.3% of the tumors, conferring a worse overall and breast cancer-specific survival [HR = 2.03; 95% confidence interval (CI), 1.65-2.48, P < 0.001], and were also found to be an independent marker of poor prognosis in estrogen receptor-positive cases (HR = 1.86; 95% CI, 1.39-2.49, P < 0.001). The mutation spectrum of TP53 varied between the breast cancer subtypes, and individual alterations showed subtype-specific association. TP53 mutations were associated with increased mortality in patients with luminal B, HER2-enriched, and normal-like tumors, but not in patients with luminal A and basal-like tumors. Similar observations were made in ICs, where mutation associated with poorer outcome in IC1, IC4, and IC5. The combined effect of TP53 mutation, TP53 LOH, and MDM2 amplification on mortality was additive. CONCLUSION: This study reveals that TP53 mutations have different clinical relevance in molecular subtypes of breast cancer, and suggests diverse roles for TP53 in the biology underlying breast cancer development.

Mutation of the salt bridge-forming residues in the ETV6-SAM domain interface blocks ETV6-NTRK3-induced cellular transformation.

The ETV6-NTRK3 (EN) chimeric oncogene is expressed in diverse tumor types. EN is generated by a t(12;15) translocation, which fuses the N-terminal SAM (sterile α-motif) domain of the ETV6 (or TEL) transcription factor to the C-terminal PTK (protein-tyrosine kinase) domain of the neurotrophin-3 receptor NTRK3. SAM domain-mediated polymerization of EN leads to constitutive activation of the PTK domain and constitutive signaling of the Ras-MAPK and PI3K-Akt pathways, which are essential for EN oncogenesis. Here we show through complementary biophysical and cellular biological techniques that mutation of Lys-99, which participates in a salt bridge at the SAM polymer interface, reduces self-association of the isolated SAM domain as well as high molecular mass complex formation of EN and abrogates the transformation activity of EN. We also show that mutation of Asp-101, the intermolecular salt bridge partner of Lys-99, similarly blocks transformation of NIH3T3 cells by EN, reduces EN tyrosine phosphorylation, inhibits Akt and Mek1/2 signaling downstream of EN, and abolishes tumor formation in nude mice. In contrast, mutations of Glu-100 and Arg-103, residues in the vicinity of the interdomain Lys-99-Asp-101 salt bridge, have little or no effect on these oncogenic characteristics of EN. Our results underscore the importance of specific electrostatic interactions for SAM polymerization and EN transformation.

The genomic and transcriptomic architecture of 2,000 breast tumours reveals novel subgroups.

The elucidation of breast cancer subgroups and their molecular drivers requires integrated views of the genome and transcriptome from representative numbers of patients. We present an integrated analysis of copy number and gene expression in a discovery and validation set of 997 and 995 primary breast tumours, respectively, with long-term clinical follow-up. Inherited variants (copy number variants and single nucleotide polymorphisms) and acquired somatic copy number aberrations (CNAs) were associated with expression in ~40% of genes, with the landscape dominated by cis- and trans-acting CNAs. By delineating expression outlier genes driven in cis by CNAs, we identified putative cancer genes, including deletions in PPP2R2A, MTAP and MAP2K4. Unsupervised analysis of paired DNA­RNA profiles revealed novel subgroups with distinct clinical outcomes, which reproduced in the validation cohort. These include a high-risk, oestrogen-receptor-positive 11q13/14 cis-acting subgroup and a favourable prognosis subgroup devoid of CNAs. Trans-acting aberration hotspots were found to modulate subgroup-specific gene networks, including a TCR deletion-mediated adaptive immune response in the 'CNA-devoid' subgroup and a basal-specific chromosome 5 deletion-associated mitotic network. Our results provide a novel molecular stratification of the breast cancer population, derived from the impact of somatic CNAs on the transcriptome.

Nucleic acid quantity and quality from paraffin blocks: defining optimal fixation, processing and DNA/RNA extraction techniques.

Although the extraction and analysis of nucleic acids from formalin-fixed paraffin-embedded tissues is a routine and growing part of pathology practice, no generally accepted recommendations exist to guide laboratories in their selection of tissue fixation, processing and DNA/RNA extraction techniques. The aim of this study was to determine how fixation method and length, paraffin embedding, processing conditions and nucleic acid extraction methods affect quality and quantity of DNA and RNA, and their performance in downstream applications. Nine tissue samples were subjected to freezing, fixation in formalin for <24 h and 7 days followed by conventional processing, and fixation in molecular fixative for <24 h and 7 days followed by rapid processing. DNA and RNA were isolated using in-house extraction and commercial kits, and assessed by PCR reactions for amplicons with varying sizes ranging from 268 to 1327 bp and one-step RT-PCR for 621 bp and 816 bp amplicons of housekeeping genes. Molecular fixative (MF) appeared to perform well under nearly all circumstances (extraction methods, fixation lengths and longer amplicons), often performing as well as frozen samples. Formalin fixation generally performed well only for shorter length amplicons and short fixation (<24 h). WaxFree kit showed consistently higher success rates for DNA and poorer rates for RNA. RecoverAll kit generally performed suboptimally in combination with prolonged formalin fixation. In conclusion, the Molecular Fixative regardless of fixation length, and the rapid tissue processing system were able to preserve large DNA and RNA fragments in paraffin blocks, making these techniques preferable for use in downstream molecular diagnostic assays.

Mll5 is required for normal spermatogenesis.

BACKGROUND: Mll5 is currently a member of the Mll family of SET domain histone methyltransferase proteins but studies have also showed that it could be part of the SET3 branch of proteins. Recently, constitutive knock out animal studies have shown that Mll5 is required for proper haematopoietic stem cell differentiation, and loss of Mll5 results in synthetic lethality for genome de-methylation. Mll5 deficient male mice are infertile and here we analyse the consequences of Mll5 deficiency for spermatogenesis. METHODOLOGY/PRINCIPAL FINDINGS: Mll5 deficient male mice, but not female mice, are infertile. Here we show using RNA in-situ hybridization that Mll5 is expressed in the germ cells of the testes of wild type mice. Consistent with the expression of Mll5, we demonstrate by electron microscopy, video microscopy and in vitro fertilisation techniques that Mll5 deficient mice have defects in terminal maturation and packaging of sperm. The defects seen include detachment of the acrosomal cap and impaired excess cytoplasm removal. Functional tests of sperm motility show a lack of progressive motility of spermatozoa from Mll5 deficient animals. None of these defects could be rescued by in vitro fertilization. Using microarray analysis we show that transcripts implicated in spermatogenesis are dysregulated. CONCLUSIONS/SIGNIFICANCE: Our data demonstrate a clear role of Mll5 in mammalian spermatogenesis at the level of terminal differentiation providing further support for its classification in the SET3 branch of proteins. Moreover, this study identifies Tlk2, Utx, Gpr64, Sult4a1, Rap2ip, Vstm2 and HoxA10 as possible Mll5 targets that together may account for the observed spermatozoa maturation defects.

The testosterone-dependent and independent transcriptional networks in the hypothalamus of Gpr54 and Kiss1 knockout male mice are not fully equivalent.

BACKGROUND: Humans and mice with loss of function mutations in GPR54 (KISS1R) or kisspeptin do not progress through puberty, caused by a failure to release GnRH. The transcriptional networks regulated by these proteins in the hypothalamus have yet to be explored by genome-wide methods. RESULTS: We show here, using 1 million exon mouse arrays (Exon 1.0 Affymetrix) and quantitative polymerase chain reaction (QPCR) validation to analyse microdissected hypothalamic tissue from Gpr54 and Kiss1 knockout mice, the extent of transcriptional regulation in the hypothalamus. The sensitivity to detect important transcript differences in microdissected RNA was confirmed by the observation of counter-regulation of Kiss1 expression in Gpr54 knockouts and confirmed by immunohistochemistry (IHC). Since Gpr54 and Kiss1 knockout animals are effectively pre-pubertal with low testosterone (T) levels, we also determined which of the validated transcripts were T-responsive and which varied according to genotype alone. We observed four types of transcriptional regulation (i) genotype only dependent regulation, (ii) T only dependent regulation, (iii) genotype and T-dependent regulation with interaction between these variables, (iv) genotype and T-dependent regulation with no interaction between these variables. The results implicate for the first time several transcription factors (e.g. Npas4, Esr2), proteases (Klk1b22), and the orphan 10-transmembrane transporter TMEM144 in the biology of GPR54/kisspeptin function in the hypothalamus. We show for the neuronal activity regulated transcription factor NPAS4, that distinct protein over-expression is seen in the hypothalamus and hippocampus in Gpr54 knockout mice. This links for the first time the hypothalamic-gonadal axis with this important regulator of inhibitory synapse formation. Similarly we confirm TMEM144 up-regulation in the hypothalamus by RNA in situ hybridization and western blot. CONCLUSIONS: Taken together, global transcriptional profiling shows that loss of GPR54 and kisspeptin are not fully equivalent in the mouse hypothalamus.

ZNF703 is a common Luminal B breast cancer oncogene that differentially regulates luminal and basal progenitors in human mammary epithelium.

The telomeric amplicon at 8p12 is common in oestrogen receptor-positive (ER+) breast cancers. Array-CGH and expression analyses of 1172 primary breast tumours revealed that ZNF703 was the single gene within the minimal amplicon and was amplified predominantly in the Luminal B subtype. Amplification was shown to correlate with increased gene and protein expression and was associated with a distinct expression signature and poor clinical outcome. ZNF703 transformed NIH 3T3 fibroblasts, behaving as a classical oncogene, and regulated proliferation in human luminal breast cancer cell lines and immortalized human mammary epithelial cells. Manipulation of ZNF703 expression in the luminal MCF7 cell line modified the effects of TGFß on proliferation. Overexpression of ZNF703 in normal human breast epithelial cells enhanced the frequency of in vitro colony-forming cells from luminal progenitors. Taken together, these data strongly point to ZNF703 as a novel oncogene in Luminal B breast cancer.