Population Of DNA Molecules Research Articles

MS1-3: Tumor Heterogeneity and Metastasis.

Abstract Massively parallel sequencing using Next Gen sequencing approaches have opened a new window into our understanding of the molecular architecture of breast cancer. A significant advantage of these technologies are that they are digital in nature because the measurement technology begins with single DNA molecules attached to a substrate. As a result, not only do these technologies provide an efficient means to identify somatic mutations, information is also provided on the frequency with which any given mutation is observed within the population of DNA molecules in the sample under analysis. The first and most striking feature of Next Gen analysis of breast cancer is the large numbers of somatic mutations that have been identified. Genome wide, thousands of mutations can be present in a single tumor. While this finding generates immense complexity from the perspective of defining biological drivers versus passengers, it does open up a new opportunity to study variations in mutation frequency. In a perfectly monoclonal tumor with a diploid genome and no contaminating normal cell DNA the mutation frequencies generated by digital sequencing will be 50% for a heterozygous mutation or 100% for homozygous mutations. Very few breast cancers subjected to whole genome sequencing exhibit this pattern, rather mutations occur with a wide variation of frequencies from the detection limit, around 5%, through to 100%. Furthermore mutation frequency clustering often occurs, a phenomenon best explained by the presence of a repertoire of mutations present in a founder clone that are present in all tumors, combined with additional less common mutations that represent a subdominant or minority population of cells that have arisen through additional clonal outgrowth. In a multiclonal model one of the most critical questions is “which clone determines the prognosis?”. One way to begin to answer this question is to sequence tumors before and after therapy and to compare primary tumors versus metastases. In our initial experience of a patient with a basal-like breast cancer in which a whole cancer genome was generated from her primary, her subsequent brain metastasis and a mouse xenograft generated from the breast primary, there was clear evidence of a shift in mutation frequency primary to metastasis, suggesting the brain metastasis arose from a minor sub-clone. The mutational profile of the xenograft more closely resembled the metastasis rather than the primary, suggesting the grafting process captured the metastatic clone (1). In this invited lecture we will present an extension of our sequencing experiences with further human progenitor-mouse xenograft comparisons, comparisons between human primaries and paired metastases and also examples of ER+ breast cancer genomes obtained before and after neoadjuvant endocrine therapy. In all these experiments multiclonality appears to be the rule rather than the exception. The clinical and biological implications of these findings will be discussed. 1. Ding L, Ellis MJ, Li S, Larson DE, Chen K, Wallis JW, et al. Genome remodelling in a basal-like breast cancer metastasis and xenograft. Nature. 2010;464(7291):999–1005. PMCID: 2872544. Citation Information: Cancer Res 2011;71(24 Suppl):Abstract nr MS1-3.

Novel subfamily of mitochondrial HMG box-containing proteins: functional analysis of Gcf1p from Candida albicans

Mitochondria of eukaryotic organisms contain populations of DNA molecules that are packed into higher-order structures called mitochondrial nucleoids (mt-nucleoids). In Saccharomyces cerevisiae, the compaction of mitochondrial DNA (mtDNA) into mt-nucleoids is mediated primarily by the high-mobility group (HMG) box-containing protein Abf2, which is an important player in stabilization and metabolism of mtDNA. Although it is evident that analogous proteins must exist in other yeast species, an apparently fast divergence rate has precluded their identification, characterization and comparative analysis. Using in silico analysis of the complete genome sequence of the pathogenic yeast Candida albicans we predicted that the ORF 19.400/19.8030 assigned as GCF1 encodes a putative mitochondrial HMG box-containing protein. In contrast to Abf2p, which contains two HMG boxes, Gcf1p contains only one C-terminal HMG box. In addition, it contains one putative coiled-coil domain with a potential role in protein dimerization. Fluorescence microscopy analysis of a C-terminally tagged Gcf1p with green fluorescent protein (GFP) revealed its mitochondrial localization in both heterologous (S. cerevisiae) and native (C. albicans) hosts. Biochemical analyses of DNA-binding properties indicate that Gcf1p is, similarly to Abf2p, a non-specific DNA-binding protein. To analyse the role of Gcf1p in mtDNA metabolism, we constructed strains lacking one functional allele of the GCF1 gene and carrying one GCF1 allele under the control of the MET3 promoter. Under repressible conditions this strain exhibited a more than 3000-fold decrease in levels of GCF1 mRNA, which was correlated with a substantial decrease in the number of mtDNA copies as well as recombination intermediates. The dramatic effect of reduced levels of Gcf1p on mtDNA metabolism indicates that the protein is involved in essential molecular transactions that relate to the mitochondrial genome.

Prospects for applying genotypic selection of somatic oncomutation to chemical risk assessment.

Genotypic selection methods detect rare sequence changes in populations of DNA molecules. These methods have been used to investigate the chemical induction of mutation and for the detection and diagnosis of cancer. The possible use of genotypic selection for improving current risk assessment practices is based on the premise that the frequency of somatic mutation is of critical importance in understanding and modeling carcinogenesis. If genotypic selection can measure the induction of specific mutations that disrupt normal cell/tissue homeostasis, then it could provide key mechanistic information for cancer risk assessment. For example, genotypic selection data might support a particular low-dose extrapolation method or characterize the relationship between rodent and human cancer risk. Strategies for evaluating the use of genotypic selection in cancer risk assessment include the concept of developing a battery of targets that detect a range of agent-specific effects. Ideal targets to examine by genotypic selection are the oncogene and tumor suppressor gene mutations frequently detected in human tumors because these are thought to represent tumor-initiating events. The most commonly occurring basepair (bp) substitutions within the ras and p53 genes are identified. Also, the battery of genotypic selection methods is defined in terms of the most important mutational specificities to include. In theory, the major basepair substitution mutations induced by 29 of 31 chemical carcinogens could be detected by analyzing three different mutations: G:C-->T:A, G:C-->A:T, and A:T-->T:A. Genotypic selection will have the greatest impact on risk assessment if measurement of spontaneous mutation is possible. Data from phenotypic selection assays suggest this corresponds to detection of mutant fractions of approximately 10(-7), and this would necessitate examining DNA samples containing >10(7) target molecules. Despite its apparent potential, considerable development and validation is needed before genotypic selection data can be applied to cancer risk assessment.

Biological properties of poly-L-lysine-DNA complexes generated by cooperative binding of the polycation.

We have evaluated the effect of NaCl concentration on the mode of binding of poly-L-lysine to DNA and the resulting structural and functional features of the condensed DNA particles using DNA precipitation, DNase I resistance, electron microscopy, and receptor-mediated gene transfer assays. At a high concentration of NaCl and in the presence of excess DNA, poly-L-lysine interacted with DNA cooperatively, fully condensing some of the DNA and leaving the rest of the DNA unbound. At low NaCl concentrations, poly-L-lysine molecules interacted with DNA in a noncooperative fashion, i.e. they bind randomly to the whole population of DNA molecules. Cooperative binding of poly-L-lysine to DNA occurred over a narrow range of NaCl concentrations, and the specific salt concentration depended on the length of the poly-L-lysine. The ability of condensed DNA to withstand digestion by DNase I was correlated with the structural features of the condensed DNA as determined by electron microscopy. Using our condensation procedure, cooperative binding of poly-L-lysine to DNA is a necessary prerequisite for the preparation of condensed DNA having a spherical shape and a diameter of 15-30 nm. Condensed DNA, containing galactosylated poly-L-lysine, was evaluated further for the extent and specificity of receptor-mediated gene transfer into HuH-7 human hepatoma cells via the asialoglycoprotein receptor. Efficient receptor-mediated transfection occurred only when condensed DNA complexes had a spherical shape with a diameter of 15-30 nm; asialofetuin, a natural ligand for the asialoglycoprotein receptor, inhibited this process by up to 90%. Our results support the importance of appropriate DNA condensation for the uptake and ultimate expression of DNA in hepatic cells.

Use of tagged random hexamer amplification (TRHA) to clone and sequence minute quantities of DNA--application to a 180 kb plasmid isolated from Sphingomonas F199.

We have developed a novel method to clone and sequence minute quantities of DNA. The method was applied to sequence a 180 kb plasmid pNL1. The first step was the production of a size distributed population of DNA molecules that were derived from the 180 kb plasmid pNL1. The first step was accomplished by a random synthesis reaction using Klenow fragment and random hexamers tagged with a T7 primer at the primer 5'-end (T7-dN6, 5'-GTAATACGACTCACTATAGGGCNNNNNN-3'. In the second step, Klenow-synthesized molecules were amplified by PCR using T7 primer (5'-GTAATACGACTCACTATAGGGC-3'). With a hundred nanograms starting plasmid DNA from pNL1, we were able to generate Klenow-synthesized molecules with sizes ranging from 28 bp to >23 kb which were detectable on an agarose gel. The Klenow-synthesized molecules were then used as templates for standard PCR with T7 primer. PCR products of sizes ranging from 0.3 to 1.3 kb were obtained for cloning and sequencing. From the same Klenow-synthesized molecules, we were also able to generate PCR products with sizes up to 23 kb by long range PCR. A total 232.5 kb sequences were obtained from 593 plasmid clones and over twenty putative genes were identified. Sequences from these 593 clones were assembled into 62 contigs and 99 individual sequence fragments with a total unique sequence of 86.3 kb.

Quantitative multiplex sequence analysis of mutational hot spots. Frequency and specificity of mutations induced by a site-specific ethenocytosine in M13 viral DNA.

We describe an assay for determining the frequency and specificity of mutations occurring at hot spots within a population of DNA molecules. The procedure consists of (a) annealing the DNA population with a labeled oligonucleotide designed to prime DNA synthesis at the mutational hot spot; (b) DNA elongation in the presence of a single dideoxynucleoside triphosphate together with 1-3 deoxynucleoside triphosphates, and (c) quantitation of all limit elongation products by high-resolution gel electrophoresis followed by autoradiography and computing densitometry. Derivation of mutational frequency and specificity over a wide range of values is demonstrated for M13 viral DNA mixtures containing defined proportions of wild-type and mutant DNAs, as well as for M13 viral DNA populations obtained by transfection of DNA bearing a defined site-specific ethenocytosine lesion. The assay is shown to yield results similar to those obtained by laborious clone-by-clone sequencing of viral progeny. The method is not affected significantly by several tested variables and appears to be suitable for use as a quantitative assay for sequence microheterogeneity at defined positions within DNA populations. Application of the methodology demonstrates that ethenocytosine, an exocyclic DNA lesion induced by carcinogens such as vinyl chloride and urethane, is a highly efficient mutagenic lesion with a mutational specificity expected for noninstructive lesions.

Replication of an adenovirus type 34 mutant dna containing tandem reiterations of the inverted terminal repeat

A mutant of human adenovirus type 34 (Ad34) has been isolated which contains DNA molecules with tandem reiterations of from two to eight copies of a 131-bp sequence within the right-sided inverted terminal repetition. Terminal heterogeneity was not eliminated by repeated plaque purifications indicating that the population of DNA molecules with various numbers of reiterations could rapidly evolve from the DNA of a single virus particle. These enlarged DNA molecules were capable of replication both in vivo and in vitro. The nucleotide sequence of the mutant Ad34 inverted terminal repetitions contained most of the essential features of the Ad origin of DNA replication. These features include the ATAATATACC sequence which is present between the highly conserved bases 9–18 in all human adenoviruses, as well as the consensus sequences for the binding of nuclear factor I and nuclear factor III. However, the reiterated sequences lacked a dG appropriately placed on the template strand to serve as a potential site for internal initiation. It appears that the rapid amplification of two to eight copies of the reiterated terminal sequences does not arise from internal initiation during replication but probably from homologous recombination.

A new method for random mutagenesis of complete genes: enzymatic generation of mutant libraries in vitro.

A new efficient in vitro mutagenesis method for the generation of complete random mutant libraries, containing all possible single base substitution mutations in a cloned gene is described. The method is based on controlled use of polymerases. Four populations of DNA molecules are first generated by primer elongation so that they terminate randomly, but always just before a known type of base (before A, C, G or T respectively). Each of the four populations is then mutagenized in a separate misincorporation reaction, where the correct base can now be omitted. The regeneration of wild-type sequences can thus be efficiently avoided. Also, the misincorporating nucleotide concentrations can be optimized to give the three possible single mutations in close to equal ratio. The mutagenesis can be precisely localized within a predetermined target region of any size, and vector sequences remain intact. We have mutagenized the DNA coding for the alpha-fragment of Escherichia coli beta-galactosidase, and identified 176 different base substitution mutations by sequencing. The present method gives mutant yields of 40-60%, when the mutants contain about one amino acid change per protein molecule. All types of base substitution mutations can be generated and deletions are rare. The efficiency of this method permits the use of relatively elaborate screening systems to isolate mutants of either structural genes or regulatory regions.

DNA ligase from Drosophila melanogaster embryos. Substrate specificity and mechanism of action.

DNA ligase has been purified to homogeneity from 6-12 h Drosophila melanogaster embryos (Rabin, B. A., Hawley, R. S., and Chase, J. W. (1986) J. Biol. Chem. 261, 10637-10645). This enzyme had an apparent Km for ATP of 1.6 microM. Of a variety of nucleotides tested, only adenosine 5'-O-(3-thio)triphosphate could substitute for ATP in the joining reaction. The enzyme was competitively inhibited by dATP, with an apparent Ki of 2.3 microM. The apparent Km for DNA using p(dT)20 annealed with poly(dA) as substrate was 1.0 microM. Studies utilizing synthetic homopolymers showed that in addition to joining DNA to DNA, this enzyme could join the 5'-phosphoryl termini of RNA to the 3'-hydroxyl termini of DNA or RNA, when they were annealed with DNA. In addition, p(dT)7U could be joined when annealed with poly(dA). No joining was detected when RNA served as the template. Drosophila DNA ligase also catalyzed the joining of oligonucleotides containing a single mismatched nucleotide at their 3'-hydroxyl termini, as well as DNA containing short, complementary 5'-protruding ends, and in the presence of polyethylene glycol 6000, blunt-ended duplex DNA. The overall reaction mechanism was shown to be identical to that of the homologous prokaryotic DNA ligases. The joining reactions catalyzed by the Drosophila and T4 DNA ligases were shown to be reversible. Incubation of superhelical closed circular DNA molecules with the purified enzymes and AMP resulted in the production of a population of DNA molecules which had lost most, if not all, of their superhelical density.

Construction of the physical map of the chloroplast DNA of Phaseolus vulgaris and localization of ribosomal and transfer RNA genes

Construction of a physical map of the chloroplast DNA from Phaseolus vulgaris showed that this circular molecule is segmentally organized into four regions. Unlike other chloroplast DNAs which have analogous organization, two single-copy regions that separate two inverted repeats have been demonstrated to exist in both relative orientations, giving rise to two populations of DNA molecules. Hybridization studies using individual rRNA and tRNA species revealed the location of a set of rRNA genes and at least seven tRNA genes in each inverted repeat region, a minimum of 17 tRNA genes in the large single-copy region and one tRNA gene in the small single-copy region. The tRNA genes code for 24 tRNA species corresponding to 16 amino acids. Comparison of this gene map with those of other chloroplast DNAs suggests that DNA sequence rearrangements, involving some tRNA genes, have occurred.

Instability and reiteration of DNA sequences within the vaccinia virus genome.

The sequence arrangement within the nontranscribed portion of the inverted terminal repetition of the vaccinia virus genome exists in quasi-stable and unstable forms that are not distinguishable on the basis of viral infectivity. The unstable forms, which composed about 20% of a serially passaged stock of virus, were recognized by terminal heterogeneity on restriction endonuclease analysis. Instead of a single terminal fragment from each end of the genome, an array of eight or more fragments differing in size by 1650-base-pair increments was detected. This feature was not eliminated by repeated plaque purification, indicating that the population of DNA molecules with various numbers of reiterations can rapidly evolve from the DNA of a single virus particle. However, at each successive round of plaque purification, about 20% of the unstable isolates revert back to the more stable form. Stable forms are characterized by the presence of a set of 13-17 tandem 70-base-pair repeats on each side of a 435-base-pair intervening sequence near both ends of the genome. In contrast, the unstable forms possess sets of tandem repeats and intervening sequences that alternate many times in series. The transition between the two genomic forms and the evolution of the unstable form appear to be mediated by recombinational events.

The structure of the termini of the DNA of Epstein-Barr virus

We have studied the DNA of Epstein-Barr virus (EBV) isolated from the B95-8 strain of that virus ( Miller and Lipman, 1973). When EBV DNA is partially digested with λ-exonuclease and allowed to reanneal, up to 50% of the full-length molecules circularize. The arrangements of nucleotide sequences containing the terminal repeats identified in this circularization experiment have been determined. Those fragments of viral DNA generated by digestion with restriction endonucleases which are terminal and contain the terminal repeats have been identified by their sensitivity to digestion of full-length DNA by λ-exonuclease and by virtue of their being partially homologous to one another. The population of DNA molecules in the B95-8 strain of EBV was found to be nonuniform. The nonuniformity results from different molecules having different numbers of a 0.37 megadalton terminal repeat at each end. About 70% of molecules have four terminal repeats at one end, while four equal classes, each comprising approximately 25% of the population, have one, two, three or four repeats at the other end. The arrangements of nucleotide sequences identified as being terminal in virion DNA were studied in the intracellular circular viral DNA of cells transformed by a single particle on EBV. All fragments produced by digestion with endonucleases and scored as being terminal in virion DNA were absent from intracellular circular DNA. An additional fragment was identified in the digests of intracellular DNA of each transformed clone. The molecular weights of the new fragments equal the sum of the molecular weights of two terminal fragments which are joined upon intracellular circularization of viral DNA.

DNA of Epstein-Barr Virus. II. Comparison of the Molecular Weights of Restriction Endonuclease Fragments of the DNA of Epstein-Barr Virus Strains and Identification of End Fragments of the B95-8 Strain

Incubation of the DNA of the B95-8 strain of Epstein-Barr virus [EBV (B95-8) DNA] with EcoRI, Hsu I, Sal I, or Kpn I restriction endonuclease yielded 8 to 15 fragments separable on 0.4% agarose gels and ranging in molecular weight from less than 1 to more than 30 x 10(6). Bam I and Bgl II yielded fragments smaller than 11 x 10(6). Preincubation of EBV (B95-8) DNA with lambda exonuclease resulted in a decrease in the Hsu I A and Sal I A and D fragments, indicating that these fragments are positioned near termini. The electrophoretic profiles of the fragments produced by cleavage of the DNA of the B95-8, HR-1, and Jijoye strains of EBV were each distinctive. The molecular weights of some EcoRI, Hsu I, and Sal I fragments from the DNA of the HR-1 strain of EBV [EBV (HR-1) DNA] and of EcoRI fragments of the DNA of the Jijoye strain of EBV were identical to that of fragments produced by cleavage of EBV (B95-8) DNA with the same enzyme, whereas others were unique to each strain. Some Hsu I, EcoRI, and Sal I fragments of EBV (HR-1) DNA and Kpn I fragments of EBV (B95-8) DNA were present in half-molar abundance relative to the majority of the fragments. In these instances, the sum of the molecular weights of the fragments was in excess of 10(8), the known molecular weight of EBV (HR-1) and (B95-8) DNA. The simplest interpretation of this finding is that each EBV (HR-1), and possibly also (B95-8), DNA preparation contains two populations of DNA molecules that differ in the arrangement of DNA sequences about a single point, such as has been described for herpes simplex virus DNA. Minor fragments could also be observed if there were more than one difference in primary structure of the DNAs. The data do not exclude more extensive heterogeneity in primary structure of the DNA of the HR-1 strain. However, the observation that the relative molar abundance of major and minor fragments of EBV (HR-1) DNA did not vary between preparations from cultures that had been maintained separately for several years favors the former hypothesis over the latter.

A quantitative theory of transfection in B. subtilis.

A theory of transfection is developed which describes three different types of experiments in transfection: The concentration dependence of transfection, transfection with marker rescue, and the mapping function in transfection crosses. The theory is applicable to transfection systems exhibiting quadratic or higher order dose response (APP1, AP50, SP82). It pictures the essential process in transfection as follows: transfecting DNA molecules, having suffered inactivating events during uptake or intracellularly, have to recombine prior to replication under elimination of these lesions. The probability for recombination, successful in this sense, is calculated as a function of the number of DNA molecules within the competent cell, the mean number of inactivating events per DNA molecule, and the crossover probability per nucleotide. Under the assumption of random distribution of inactivating events over the population of DNA molecules and homogeneous crossover probabilities the theory explains on a quantitative basis a number of experimental observations in transfection, as for instance the relative efficiencies of different helper phage in transfection with marker rescue, the third order concentration dependence in SP50 transfection, and the high recombination frequencies observed in transfection crosses.

The evolution of new species of viral DNA during serial passage of simian virus 40 at high multiplicity

Viral DNA isolated during serial undiluted passage of SV40 was analysed by electrophoresis of fragments produced by H. influenzae restriction endonuclease and by DNA-DNA hybridization. Compared to parental SV40 DNA, early passage DNA showed alterations at many sites in the population of DNA molecules, although certain sites were more frequently altered than others as evidenced by the reduced yield of some original DNA fragments and the appearance of new fragments. With continued passage, new dominant species of viral DNA evolved which yielded simple restriction endonuclease digest patterns. DNA molecules from all passages retained SV40 DNA sequences; however, the percentage of total DNA sequences which corresponded to SV40 sequences diminished to about 30% or less in late passage DNA. Repetitive host cell DNA sequences were detected in a high percentage of late passage DNA molecules, but accounted for 20% or less of total sequences present. Therefore the new species of viral DNA which have evolved probably contain predominantly nonrepetitive host DNA sequences.

Genetic Homologies and Drift within Population of DNA Molecules

Genetic Homologies and Drift within Population of DNA Molecules

THE RELATIVE HOMOGENEITY OF MICROBIAL DNA

Density-gradient centrifugation has revealed that the population of DNA molecules from E. coli is relatively homogeneous with respect to buoyant density in a solution of cesium chloride.(1) Because of the notably small apparent atomic volume of nitrogen in aqueous solutions of organic compounds(2) and because the guanine-cytosine base pair is more rich in nitrogen than is the adenine-thymine pair, it was considered that the density homogeneity among E. coli DNA molecules might reflect a high degree of homogeneity with respect to base composition. To investigate this possibility, an examination has been made of the relationship between buoyant density in cesium chloride solution and base composition of DNA from various sources.