Michael Wigler

Tumour evolution inferred by single-cell sequencing

Genomic analysis provides insights into the role of copy number variation in disease, but most methods are not designed to resolve mixed populations of cells. In tumours, where genetic heterogeneity is common, very important information may be lost that would be useful for reconstructing evolutionary history. Here we show that with flow-sorted nuclei, whole genome amplification and next generation sequencing we can accurately quantify genomic copy number within an individual nucleus. We apply single-nucleus sequencing to investigate tumour population structure and evolution in two human breast cancer cases. Analysis of 100 single cells from a polygenomic tumour revealed three distinct clonal subpopulations that probably represent sequential clonal expansions. Additional analysis of 100 single cells from a monogenomic primary tumour and its liver metastasis indicated that a single clonal expansion formed the primary tumour and seeded the metastasis. In both primary tumours, we also identified an unexpectedly abundant subpopulation of genetically diverse ‘pseudodiploid’ cells that do not travel to the metastatic site. In contrast to gradual models of tumour progression, our data indicate that tumours grow by punctuated clonal expansions with few persistent intermediates.

Transfer of purified herpes virus thymidine kinase gene to cultured mouse cells.

Treatment of Ltk−, mouse L cells deficient in thymidine kinase (tk), with Bam I restriction endonuclease cleaved DNA from herpes simplex virus-1 (HSV-1) produced tk+ clones with a frequency of 10−6/2 μg of HSV-1 DNA. Untreated cells or cells treated with Eco RI restriction endonuclease fragments produced no tk+ clones under the same conditions. The thymidine kinase activities of four independently derived clones were characterized by biochemical and serological techniques. By these criteria, the tk activities were found to be identical to HSV-1 tk and different from host wildtype tk. The tk+ phenotype was stable over several hundred cell generations, although the rate of reversion to the tk− phenotype, as judged by cloning efficiency in the presence of bromodeoxyuridine, was high (1–5 × 10−3). HSV-1 DNA Bam restriction fragments were separated by gel electrophoresis, and virtually all activity, as assayed by transfection, was found to reside in a 3.4 kb fragment. Transformation efficiency with the isolated fragment is 20 fold higher per gene equivalent than with the unfractionated total Bam digest. These results prove the usefulness of transfection assays as a means for the bioassay and isolation of restriction fragments carrying specific genetic information. Cells expressing HSV-1 tk may also provide a useful model system for the detailed analysis of eucaryotic and viral gene regulation.

Transformation of mammalian cells with genes from procaryotes and eucaryotes

Abstract We have stably transformed mammalian cells with precisely defined procaryotic and eucaryotic genes for which no selective criteria exist. The addition of a purified viral thymidine kinase (tk) gene to mouse cells lacking this enzyme results in the appearance of stable transformants which can be selected by their ability to grow in HAT. These biochemical transformants may represent a subpopulation of competent cells which are likely to integrate other unlinked genes at frequencies higher than the general population. Co-transformation experiments were therefore performed with the viral tk gene and bacteriophage ΦX174, plasmid pBR322 or the cloned chromosomal rabbit β-globin gene sequences. Tk + transformants were cloned and analyzed for co-transfer of additional DNA sequences by blot hybridization. In this manner, we have identified mouse cell lines which contain multiple copies of 4)X, pBR322 and the rabbit β-globin gene sequences. The ΦX co-transformants were studied in greatest detail. The frequency of co-transformation is high: 15 of 16 tk + transformants contain the ΦX sequences. Selective pressure was required to identify co-transformants. From one to more than fifty ΦX sequences are integrated into high molecular weight nuclear DNA isolated from independent clones. Analysis of subclones demonstrates that the ΦX genotype is stable through many generations in culture. This co-transformation system should allow the introduction and stable integration of virtually any defined gene into cultured cells. Ligation to either viral vectors or selectable biochemical markers is not required.

Biochemical transfer of single-copy eucaryotic genes using total cellular DNA as donor

Previous studies from our laboratories have demonstrated the feasibility of transferring the thymidine kinase (tk) gene from restriction endonuclease-generated fragments of herpes simplex virus (HSV) DNA to cultured mammalian cells. In this study, high molecular weight DNA from cells containing only one copy of the HSV gene coding for tk was successfully used to transform L+K-cells to the tk+ phenotype. The acquired phenotype was demonstrated to be donor-derived by analysis of the electrophoretic mobility of the tk activity, and the presence of HSV DNA sequences in the recipient cells was demonstrated. In companion experiments, we used high molecular weight DNA derived from tissues and cultured cells of a variety of species to transfer tk activity. The tk+ mouse cells transformed with human DNA were shown to express human type tk activity as determined by isoelectric focusing.

Identification and Validation of Oncogenes in Liver Cancer Using an Integrative Oncogenomic Approach

The heterogeneity and instability of human tumors hamper straightforward identification of cancer-causing mutations through genomic approaches alone. Herein we describe a mouse model of liver cancer initiated from progenitor cells harboring defined cancer-predisposing lesions. Genome-wide analyses of tumors in this mouse model and in human hepatocellular carcinomas revealed a recurrent amplification at mouse chromosome 9qA1, the syntenic region of human chromosome 11q22. Gene-expression analyses delineated cIAP1, a known inhibitor of apoptosis, and Yap, a transcription factor, as candidate oncogenes in the amplicon. In the genetic context of their amplification, both cIAP1 and Yap accelerated tumorigenesis and were required to sustain rapid growth of amplicon-containing tumors. Furthermore, cIAP1 and Yap cooperated to promote tumorigenesis. Our results establish a tractable model of liver cancer, identify two oncogenes that cooperate by virtue of their coamplification in the same genomic locus, and suggest an efficient strategy for the annotation of human cancer genes.

In yeast, RAS proteins are controlling elements of adenylate cyclase

S. cerevisiae strains containing RAS2val19, a RAS2 gene with a missense mutation analogous to one that activates the transforming potential of mammalian ras genes, have growth and biochemical properties strikingly similar to yeast strains carrying IAC or bcy1. Yeast strains carrying the IAC mutation have elevated levels of adenylate cyclase activity. bcy1 is a mutation that suppresses the lethality in adenylate cyclase deficient yeast. Yeast strains deficient in RAS function exhibit properties similar to adenylate cyclase deficient yeast. bcy1 suppresses lethality in ras1- ras2- yeast. Compared to wild-type yeast strains, intracellular cyclic AMP levels are significantly elevated in RAS2val19 strains, significantly depressed in ras2- strains, and virtually undetectable in ras1- ras2- bcy1 strains. Membranes from ras1- ras2- bcy1 yeast lack the GTP-stimulated adenylate cyclase activity present in membranes from wild-type cells, and membranes from RAS2val19 yeast strains have elevated levels of an apparently GTP-independent adenylate cyclase activity. Mixing membranes from ras1- ras2- yeast with membranes from adenylate cyclase deficient yeast reconstitutes a GTP-dependent adenylate cyclase.

Purification of a RAS-responsive adenylyl cyclase complex from Saccharomyces cerevisiae by use of an epitope addition method.

We developed a method for immunoaffinity purification of Saccharomyces cerevisiae adenylyl cyclase based on creating a fusion with a small peptide epitope. Using oligonucleotide technology to encode the peptide epitope we constructed a plasmid that expressed the fusion protein from the S. cerevisiae alcohol dehydrogenase promoter ADH1. A monoclonal antibody previously raised against the peptide was used to purify adenylyl cyclase by affinity chromatography. The purified enzyme appeared to be a multisubunit complex consisting of the 200-kilodalton adenylyl cyclase fusion protein and an unidentified 70-kilodalton protein. The purified protein could be activated by RAS proteins. Activation had an absolute requirement for a guanine nucleoside triphosphate.

The contribution of de novo coding mutations to autism spectrum disorder

Whole exome sequencing has proven to be a powerful tool for understanding the genetic architecture of human disease. Here we apply it to more than 2,500 simplex families, each having a child with an autistic spectrum disorder. By comparing affected to unaffected siblings, we show that 13% of de novo missense mutations and 43% of de novo likely gene-disrupting (LGD) mutations contribute to 12% and 9% of diagnoses, respectively. Including copy number variants, coding de novo mutations contribute to about 30% of all simplex and 45% of female diagnoses. Almost all LGD mutations occur opposite wild-type alleles. LGD targets in affected females significantly overlap the targets in males of lower intelligence quotient (IQ), but neither overlaps significantly with targets in males of higher IQ. We estimate that LGD mutation in about 400 genes can contribute to the joint class of affected females and males of lower IQ, with an overlapping and similar number of genes vulnerable to contributory missense mutation. LGD targets in the joint class overlap with published targets for intellectual disability and schizophrenia, and are enriched for chromatin modifiers, FMRP-associated genes and embryonically expressed genes. Most of the significance for the latter comes from affected females.

Multiple ras functions can contribute to mammalian cell transformation

We have developed a generalized approach, using two hybrid interactions, to isolate Ha-Ras effector loop mutations that separate the ability of Ha-Ras to interact with different downstream effectors. These mutations attenuate or eliminate Ha-ras(G12V) transformation of mammalian cells, but retain complementary activity, as demonstrated by synergistic induction of foci of growth-transformed cells, and by the ability to activate different downstream components. The transformation defect of Ha-ras(G12V, E37G) is rescued by a mutant, raf1, that restores interaction. These results indicate that multiple cellular components, including Raf1, are activated by Ha-Ras and contribute to Ha-Ras-induced mammalian cell transformation.

The NF1 locus encodes a protein functionally related to mammalian GAP and yeast IRA proteins

The von Recklinghausen neurofibromatosis locus, NF1, encodes a protein with homology restricted to the catalytic region of the RAS GTPase-activating protein, GAP, and with extensive homology to the IRA1 and IRA2 gene products of the yeast S. cerevisiae. A segment of the NF1 cDNA gene, expressed in yeast, can complement loss of IRA function and can inhibit both wild-type and mutant activated human H-ras genes that are coexpressed in yeast. Yeast expressing the NF1 segment have increased H-ras GTPase-stimulating activity. These studies indicate that the NF1 gene product can interact with RAS proteins and demonstrate structural and functional similarities and differences among the GAP, IRA1, IRA2, and NF1 proteins.