Single-cell genomic variation induced by mutational processes in cancer

How cell-to-cell copy number alterations that underpin genomic instability¹ in human cancers drive genomic and phenotypic variation, and consequently the evolution of cancer², remains understudied. Here, by applying scaled single-cell whole-genome sequencing³ 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.

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.

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.

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.

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.

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.