Clonal Decomposition and DNA Replication States Defined by Scaled Single-Cell Genome Sequencing.

Accurate measurement of clonal genotypes, mutational processes, and replication states from individual tumor-cell genomes will facilitate improved understanding of tumor evolution. We have developed DLP+, a scalable single-cell whole-genome sequencing platform implemented using commodity instruments, image-based object recognition, and open source computational methods. Using DLP+, we have generated a resource of 51,926 single-cell genomes and matched cell images from diverse cell types including cell lines, xenografts, and diagnostic samples with limited material. From this resource we have defined variation in mitotic mis-segregation rates across tissue types and genotypes. Analysis of matched genomic and image measurements revealed correlations between cellular morphology and genome ploidy states. Aggregation of cells sharing copy number profiles allowed for calculation of single-nucleotide resolution clonal genotypes and inference of clonal phylogenies and avoided the limitations of bulk deconvolution. Finally, joint analysis over the above features defined clone-specific chromosomal aneuploidy in polyclonal populations.

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.

Clonal fitness inferred from time-series modelling of single-cell cancer genomes.

Progress in defining genomic fitness landscapes in cancer, especially those defined by copy number alterations (CNAs), has been impeded by lack of time-series single-cell sampling of polyclonal populations and temporal statistical models1-7. Here we generated 42,000 genomes from multi-year time-series single-cell whole-genome sequencing of breast epithelium and primary triple-negative breast cancer (TNBC) patient-derived xenografts (PDXs), revealing the nature of CNA-defined clonal fitness dynamics induced by TP53 mutation and cisplatin chemotherapy. Using a new Wright-Fisher population genetics model8,9 to infer clonal fitness, we found that TP53 mutation alters the fitness landscape, reproducibly distributing fitness over a larger number of clones associated with distinct CNAs. Furthermore, in TNBC PDX models with mutated TP53, inferred fitness coefficients from CNA-based genotypes accurately forecast experimentally enforced clonal competition dynamics. Drug treatment in three long-term serially passaged TNBC PDXs resulted in cisplatin-resistant clones emerging from low-fitness phylogenetic lineages in the untreated setting. Conversely, high-fitness clones from treatment-naive controls were eradicated, signalling an inversion of the fitness landscape. Finally, upon release of drug, selection pressure dynamics were reversed, indicating a fitness cost of treatment resistance. Together, our findings define clonal fitness linked to both CNA and therapeutic resistance in polyclonal tumours.

Dynamics of genomic clones in breast cancer patient xenografts at single-cell resolution.

Human cancers, including breast cancers, comprise clones differing in mutation content. Clones evolve dynamically in space and time following principles of Darwinian evolution, underpinning important emergent features such as drug resistance and metastasis. Human breast cancer xenoengraftment is used as a means of capturing and studying tumour biology, and breast tumour xenografts are generally assumed to be reasonable models of the originating tumours. However, the consequences and reproducibility of engraftment and propagation on the genomic clonal architecture of tumours have not been systematically examined at single-cell resolution. Here we show, using deep-genome and single-cell sequencing methods, the clonal dynamics of initial engraftment and subsequent serial propagation of primary and metastatic human breast cancers in immunodeficient mice. In all 15 cases examined, clonal selection on engraftment was observed in both primary and metastatic breast tumours, varying in degree from extreme selective engraftment of minor (<5% of starting population) clones to moderate, polyclonal engraftment. Furthermore, ongoing clonal dynamics during serial passaging is a feature of tumours experiencing modest initial selection. Through single-cell sequencing, we show that major mutation clusters estimated from tumour population sequencing relate predictably to the most abundant clonal genotypes, even in clonally complex and rapidly evolving cases. Finally, we show that similar clonal expansion patterns can emerge in independent grafts of the same starting tumour population, indicating that genomic aberrations can be reproducible determinants of evolutionary trajectories. Our results show that measurement of genomically defined clonal population dynamics will be highly informative for functional studies using patient-derived breast cancer xenoengraftment.

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.

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.

Loss of MLL5 results in pleiotropic hematopoietic defects, reduced neutrophil immune function, and extreme sensitivity to DNA demethylation.

MLL5 is a divergent member of the Drosophila Trithorax-related (SET) domain and plant homeodomain (PHD) domain-containing chromatin regulators that are involved in the regulation of transcriptional "memory" during differentiation. Human MLL5 is located on chromosome 7q22, which frequently is deleted in myeloid leukemias, suggesting a possible role in hemopoiesis. To address this question, we generated a loss-of-function allele (Mll5(tm1Apa)) in the murine Mll5 locus. Unlike other Mll genes, Mll5(tm1Apa) homozygous mice are viable but display defects in immunity and hematopoiesis. First, Mll5(tm1Apa) homozygous mice show increased susceptibility to spontaneous eye infections, associated with a cell-autonomous impairment of neutrophil function. Second, Mll5(tm1Apa/tm1Apa) mice exhibit a mild impairment of erythropoiesis. Third, Mll5(tm1Apa/tm1Apa) hematopoietic stem cells (HSCs) have impaired competitive repopulating capacity both under normal conditions and when subjected to self-renewal stimulation by NUP98-HOXA10. Fourth, Mll5(tm1Apa) homozygous HSCs show a dramatic sensitivity to DNA demethylation-induced differentiation (5-azadeoxycytidine). Taken together, our data show that MLL5 is involved in terminal myeloid differentiation and the regulation of HSC self-renewal by a mechanism that involves DNA methylation. These data warrant investigation of MLL5 expression levels as a predictive marker of demethylating-agent response in patients with myelodysplastic syndromes and leukemias and identify MLL5 as a key regulator of normal hematopoiesis.

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

SAGE reveals expression of Wnt signalling pathway members during mouse prostate development.

To identify genes and pathways not previously implicated in the mesenchymal-epithelial (M/E) interactions that are critical for normal mouse prostate development, we constructed six serial analysis of gene expression (SAGE) libraries. Bioinformatic analyses revealed expression of various members of numerous signalling pathways and the differential expression of several members of the wingless-related MMTV integration site (Wnt) signalling pathway. This pathway has not been previously implicated in prostate development thus expression of selected Wnt pathway members in the developing prostate was confirmed by RT-qPCR. Of particular interest, an antagonist of the Wnt pathway, secreted frizzled related protein 2 (Sfrp2), was highly expressed in the early prostate libraries and down regulated at later developmental stages. The expression levels of four Wnt ligands reported to interact with Sfrp2 were, therefore, examined by RT-qPCR. We found that only Wnt4 transcripts were detectable in the developing prostate. Expression of Sfrp2 was validated using RT-qPCR and localization of Sfrp2 transcripts and protein was carried out using in situ hybridization and immunofluorescence, respectively. These studies provide the first evidence that Wnt pathway members are expressed in the developing prostate. Functional analyses are now required to establish the biological significance of this observation.

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.

Antimicrobial and antifungal activities of a novel cationic antimicrobial peptide, omiganan, in experimental skin colonisation models.

Omiganan pentahydrochloride is a novel, synthetic, cationic, antimicrobial peptide that is being developed for the prevention of catheter-related infections and the treatment of acne and rosacea. In this study, the efficacy of topical omiganan gel was evaluated in two skin colonisation models (ex vivo pig skin and in vivo guinea pig skin). When tested in the ex vivo pig skin colonisation model, omiganan 0.1-2% gels exhibited potent dose-dependent activity against gram-positive and gram-negative bacteria and yeasts; the maximum effect was observed at 1-2%. No significant difference was noted in activity toward meticillin-resistant and meticillin-sensitive Staphylococcus aureus, and drug activity was not affected by the inoculum size. The antimicrobial activity of omiganan 1% gel was rapid, with a 2.7 log(10)colony-forming unit (CFU)/site reduction in Staphylococcus epidermidis counts at 1 h post application and a 5.2 log(10)CFU/site reduction at 24 h. Additional studies in the guinea pig skin colonisation model confirmed the potent antimicrobial and antifungal activities of omiganan 1% gel. In conclusion, omiganan gels have been demonstrated to be rapidly bactericidal and fungicidal, with significant dose-dependent activity against a broad spectrum of infectious organisms. These results further confirm that the drug has the potential as a topical antimicrobial agent.

Identification of transcripts with enriched expression in the developing and adult pancreas.

BACKGROUND: Despite recent advances, the transcriptional hierarchy driving pancreas organogenesis remains largely unknown, in part due to the paucity of comprehensive analyses. To address this deficit we generated ten SAGE libraries from the developing murine pancreas spanning Theiler stages 17-26, making use of available Pdx1 enhanced green fluorescent protein (EGFP) and Neurog3 EGFP reporter strains, as well as tissue from adult islets and ducts. RESULTS: We used a specificity metric to identify 2,536 tags with pancreas-enriched expression compared to 195 other mouse SAGE libraries. We subsequently grouped co-expressed transcripts with differential expression during pancreas development using K-means clustering. We validated the clusters first using quantitative real time PCR and then by analyzing the Theiler stage 22 pancreas in situ hybridization staining patterns of over 600 of the identified genes using the GenePaint database. These were then categorized into one of the five expression domains within the developing pancreas. Based on these results we identified a cascade of transcriptional regulators expressed in the endocrine pancreas lineage and, from this, we developed a predictive regulatory network describing beta-cell development. CONCLUSION: Taken together, this work provides evidence that the SAGE libraries generated here are a valuable resource for continuing to elucidate the molecular mechanisms regulating pancreas development. Furthermore, our studies provide a comprehensive analysis of pancreas development, and insights into the regulatory networks driving this process are revealed.

Dynamics of expression of growth differentiation factor 15 in normal and PIN development in the mouse.

Growth differentiation factor (GDF15) is a distant member of the transforming growth factor-beta superfamily, a diverse group of structurally related proteins that exert multiple effects on cell fate such as on cell growth and differentiation but little is known about GDF15 in these processes. Previously we observed the mature GDF15 to be associated with human prostate carcinogenesis hence prompting us to study GDF15 further. Here we report gdf15 expression both at the RNA and protein levels, in normal prostatic tissues of wild type (wt) and prostatic intraepithelial neoplasia (PIN) of transgenic (Tg) 12T-7s model mice during embryonic, postnatal, and adult prostate formation up to 15 weeks after birth. Dynamic changes in expression, at both the mRNA and protein level, correlated with cell proliferation and differentiation during distinct phases of normal mouse prostate development and alterations in the dynamics of gdf15 expression correlated with the changes in development resulting in PIN formation. Most notably mature gdf15 protein was significantly elevated during hyperplasia and PIN development. Changes in the protein levels did not always correlate well with the mRNA levels. This was more prominent during PIN than during normal prostate development suggesting that this may also be an indicator of disturbed regulation of gdf15 in PIN. We propose that gdf15 is a growth factor with dual function either promoting proliferation or growth arrest and differentiation due most likely to differences in cellular differentiation. Because of the differentiation defect in PIN its epithelium no longer responds to gdf15 by cellular growth arrest as does the normal epithelium and gdf may even stimulate proliferation. The data supports our hypothesis that GDF15 plays a role in the early stages of human prostate cancer.

A mouse atlas of gene expression: large-scale digital gene-expression profiles from precisely defined developing C57BL/6J mouse tissues and cells.

We analyzed 8.55 million LongSAGE tags generated from 72 libraries. Each LongSAGE library was prepared from a different mouse tissue. Analysis of the data revealed extensive overlap with existing gene data sets and evidence for the existence of approximately 24,000 previously undescribed genomic loci. The visual cortex, pancreas, mammary gland, preimplantation embryo, and placenta contain the largest number of differentially expressed transcripts, 25% of which are previously undescribed loci.