Genomic evolution during a 10,000-generation experiment with bacteria

Molecular methods are used widely to measure genetic diversity within populations and determine relationships among species. However, it is difficult to observe genomic evolution in action because these dynamics are too slow in most organisms. To overcome this limitation, we sampled genomes from populations of Escherichia coli evolving in the laboratory for 10,000 generations. We analyzed the genomes for restriction fragment length polymorphisms (RFLP) using seven insertion sequences (IS) as probes; most polymorphisms detected by this approach reflect rearrangements (including transpositions) rather than point mutations. The evolving genomes became increasingly different from their ancestor over time. Moreover, tremendous diversity accumulated within each population, such that almost every individual had a different genetic fingerprint after 10,000 generations. As has been often suggested, but not previously shown by experiment, the rates of phenotypic and genomic change were discordant, both across replicate populations and over time within a population. Certain pivotal mutations were shared by all descendants in a population, and these are candidates for beneficial mutations, which are rare and difficult to find. More generally, these data show that the genome is highly dynamic even over a time scale that is, from an evolutionary perspective, very brief.

Quantifying single gene copy number by measuring fluorescent probe lengths on combed genomic DNA

An approach was developed for the quantification of subtle gains and losses of genomic DNA. The approach relies on a process called molecular combing. Molecular combing consists of the extension and alignment of purified molecules of genomic DNA on a glass coverslip. It has the advantage that a large number of genomes can be combed per coverslip, which allows for a statistically adequate number of measurements to be made on the combed DNA. Consequently, a high-resolution approach to mapping and quantifying genomic alterations is possible. The approach consists of applying fluorescence hybridization to the combed DNA by using probes to identify the amplified region. Measurements then are made on the linear hybridization signals to ascertain the region's exact size. The reliability of the approach first was tested for low copy number amplifications by determining the copy number of chromosome 21 in a normal and trisomy 21 cell line. It then was tested for high copy number amplifications by quantifying the copy number of an oncogene amplified in the tumor cell line GTL-16. These results demonstrate that a wide range of amplifications can be accurately and reliably quantified. The sensitivity and resolution of the approach likewise was assessed by determining the copy number of a single allele (160 kb) alteration.

Molecular cytogenetic delineation of a novel critical genomic region in chromosome bands 11q22.3-923.1 in lymphoproliferative disorders.

Aberrations of the long arm of chromosome 11 are among the most common chromosome abnormalities in lymphoproliferative disorders (LPD). Translocations involving BCL1 at 11q13 are strongly associated with mantle cell lymphoma. other nonrandom aberrations, especially deletions and, less frequently, translocations, involving bands 11q21-923 have been identified by chromosome banding analysis. To date, the critical genomic segment and candidate genes involved in these deletions have not been identified. In the present study, we have analyzed tumors from 43 patients with LPD (B-cell chronic lymphocytic leukemia, n = 40; mantle cell lymphoma, n = 3) showing aberrations of bands 11q21-923 by fluorescence in situ hybridization. As probes we used Alu-PCR products from 17 yeast artificial chromosome clones spanning chromosome bands 11q14.3-923.3, including a panel of yeast artificial chromosome clones recognizing a contiguous genomic DNA fragment of approximately 9-10 Mb in bands 11q22.3-923.3. In the 41 tumors exhibiting deletions, we identified a commonly deleted segment in band 11q22.3-923.1; this region is approximately 2-3 Mb in size and contains the genes coding for ATM (ataxia telangiectasia mutated), RDX (radixin), and FDX1 (ferredoxin 1). Furthermore, two translocation break-points were localized to a 1.8-Mb genomic fragment contained within the commonly deleted segment. Thus, we have identified a single critical region of 2-3 Mb in size in which 11q14-923 aberrations in LPD cluster. This provides the basis for the identification of the gene(s) at 11q22.3-923.1 that are involved in the pathogenesis of LPD.

The genomic and physical organization of Ty1-copia-like sequences as a component of large genomes in Pinus elliottii var. elliottii and other gymnosperms.

A DNA sequence, TPE1, representing the internal domain of a Ty1-copia retroelement, was isolated from genomic DNA of Pinus elliottii Engelm. var. elliottii (slash pine). Genomic Southern analysis showed that this sequence, carrying partial reverse transcriptase and integrase gene sequences, is highly amplified within the genome of slash pine and part of a dispersed element >4.8 kbp. Fluorescent in situ hybridization to metaphase chromosomes shows that the element is relatively uniformly dispersed over all 12 chromosome pairs and is highly abundant in the genome. It is largely excluded from centromeric regions and intercalary chromosomal sites representing the 18S-5.8S-25S rRNA genes. Southern hybridization with specific DNA probes for the reverse transcriptase gene shows that TPE1 represents a large subgroup of heterogeneous Ty1-copia retrotransposons in Pinus species. Because no TPE1 transcription could be detected, it is most likely an inactive element--at least in needle tissue. Further evidence for inactivity was found in recombinant reverse transcriptase and integrase sequences. The distribution of TPE1 within different gymnosperms that contain Ty1-copia group retrotransposons, as shown by a PCR assay, was investigated by Southern hybridization. The TPE1 family is highly amplified and conserved in all Pinus species analyzed, showing a similar genomic organization in the three- and five-needle pine species investigated. It is also present in spruce, bald cypress (swamp cypress), and in gingko but in fewer copies and a different genomic organization.

Transient excess of MYC activity can elicit genomic instability and tumorigenesis

Overexpression of the MYC protooncogene has been implicated in the genesis of diverse human tumors. Tumorigenesis induced by MYC has been attributed to sustained effects on proliferation and differentiation. Here we report that MYC may also contribute to tumorigenesis by destabilizing the cellular genome. A transient excess of MYC activity increased tumorigenicity of Rat1A cells by at least 50-fold. The increase persisted for >30 days after the return of MYC activity to normal levels. The brief surfeit of MYC activity was accompanied by evidence of genomic instability, including karyotypic abnormalities, gene amplification, and hypersensitivity to DNA-damaging agents. MYC also induced genomic destabilization in normal human fibroblasts, although these cells did not become tumorigenic. Stimulation of Rat1A cells with MYC accelerated their passage through G1/S. Moreover, MYC could force normal human fibroblasts to transit G1 and S after treatment with N-(phosphonoacetyl)-l-aspartate (PALA) at concentrations that normally lead to arrest in S phase by checkpoint mechanisms. Instead, the cells subsequently appeared to arrest in G2. We suggest that the accelerated passage through G1 was mutagenic but that the effect of MYC permitted a checkpoint response only after G2 had been reached. Thus, MYC may contribute to tumorigenesis through a dominant mutator effect.

Extent of genomic rearrangement after genome duplication in yeast

Whole-genome duplication approximately 108 years ago was proposed as an explanation for the many duplicated chromosomal regions in Saccharomyces cerevisiae. Here we have used computer simulations and analytic methods to estimate some parameters describing the evolution of the yeast genome after this duplication event. Computer simulation of a model in which 8% of the original genes were retained in duplicate after genome duplication, and 70–100 reciprocal translocations occurred between chromosomes, produced arrangements of duplicated chromosomal regions very similar to the map of real duplications in yeast. An analytical method produced an independent estimate of 84 map disruptions. These results imply that many smaller duplicated chromosomal regions exist in the yeast genome in addition to the 55 originally reported. We also examined the possibility of determining the original order of chromosomal blocks in the ancestral unduplicated genome, but this cannot be done without information from one or more additional species. If the genome sequence of one other species (such as Kluyveromyces lactis) were known it should be possible to identify 150–200 paired regions covering the whole yeast genome and to reconstruct approximately two-thirds of the original order of blocks of genes in yeast. Rates of interchromosome translocation in yeast and mammals appear similar despite their very different rates of homologous recombination per kilobase.

Genetic analysis using genomic representations

Analysis of the genetic changes in human tumors is often problematical because of the presence of normal stroma and the limited availability of pure tumor DNA. However, large amounts of highly reproducible “representations” of tumor and normal genomes can be made by PCR from nanogram amounts of restriction endonuclease cleaved DNA that has been ligated to oligonucleotide adaptors. We show here that representations are useful for many types of genetic analyses, including measuring relative gene copy number, loss of heterozygosity, and comparative genomic hybridization. Representations may be prepared even from sorted nuclei from fixed and archived tumor biopsies.

Discrimination of the closely related A and D genomes of the hexaploid oat Avena sativa L.

A satellite DNA sequence, As120a, specific to the A-genome chromosomes in the hexaploid oat, Avena sativa L., was isolated by subcloning a fragment with internal tandem repeats from a plasmid, pAs120, that had been obtained from an Avena strigosa (As genome) genomic library. Southern and in situ hybridization showed that sequences with homology to sequences within pAs120 were dispersed throughout the genome of diploid (A and C genomes), tetraploid (AC genomes), and hexaploid (ACD genomes) Avena species. In contrast, sequences homologous to As120a were found in two A-genome species (A. strigosa and Avena longiglumis) and in the hexaploid A. sativa whereas this sequence was little amplified in the tetraploid Avena murphyi and was absent in the remaining A- and C-genome diploid species. In situ hybridization of pAs120a to hexaploid oat species revealed the distribution of elements of the As120a repeated family over both arms of 14 of 42 chromosomes of this species. By using double in situ hybridization with pAs120a and a C genome-specific probe, three sets of 14 chromosomes were revealed corresponding to the A, C, and D genomes of the hexaploid species. Simultaneous in situ hybridizations with pAs120a and ribosomal probes were used to assign the SAT chromosomes of hexaploid species to their correct genomes. This work reports a sequence able to distinguish between the closely related A and D genomes of hexaploid oats. This sequence offers new opportunities to analyze the relationships of Avena species and to explore the possible evolution of various polyploid oat species.

Human ATP-binding cassette transporter 1 (ABC1): Genomic organization and identification of the genetic defect in the original Tangier disease kindred

Tangier disease is characterized by low serum high density lipoproteins and a biochemical defect in the cellular efflux of lipids to high density lipoproteins. ABC1, a member of the ATP-binding cassette family, recently has been identified as the defective gene in Tangier disease. We report here the organization of the human ABC1 gene and the identification of a mutation in the ABC1 gene from the original Tangier disease kindred. The organization of the human ABC1 gene is similar to that of the mouse ABC1 gene and other related ABC genes. The ABC1 gene contains 49 exons that range in size from 33 to 249 bp and is over 70 kb in length. Sequence analysis of the ABC1 gene revealed that the proband for Tangier disease was homozygous for a deletion of nucleotides 3283 and 3284 (TC) in exon 22. The deletion results in a frameshift mutation and a premature stop codon starting at nucleotide 3375. The product is predicted to encode a nonfunctional protein of 1,084 aa, which is approximately half the size of the full-length ABC1 protein. The loss of a Mnl1 restriction site, which results from the deletion, was used to establish the genotype of the rest of the kindred. In summary, we report on the genomic organization of the human ABC1 gene and identify a frameshift mutation in the ABC1 gene of the index case of Tangier disease. These results will be useful in the future characterization of the structure and function of the ABC1 gene and the analysis of additional ABC1 mutations in patients with Tangier disease.

Probes bound to myosin Cys-707 rotate during length transients in contraction

It is widely conjectured that muscle shortens because portions of myosin molecules (the “cross-bridges”) impel the actin filament to which they transiently attach and that the impulses result from rotation of the cross-bridges. Crystallography indicates that a cross-bridge is articulated–consisting of a globular catalytic/actin-binding domain and a long lever arm that may rotate. Conveniently, a rhodamine probe with detectable attitude can be attached between the globular domain and the lever arm, enabling the observer to tell whether the anchoring region rotates. Well-established signature effects observed in shortening are tension changes resulting from the sudden release or quick stretch of active muscle fibers. In this investigation we found that closely correlated with such tension changes are changes in the attitude of the rhodamine probes. This correlation strongly supports the conjecture about how shortening is achieved.

Panhandle PCR strategy to amplify MLL genomic breakpoints in treatment-related leukemias

Panhandle PCR amplifies genomic DNA with known 5′ and unknown 3′ sequences from a template with an intrastrand loop schematically shaped like a pan with a handle. We used panhandle PCR to clone MLL genomic breakpoints in two pediatric treatment-related leukemias. The karyotype in a case of treatment-related acute lymphoblastic leukemia showed the t(4;11)(q21;q23). Panhandle PCR amplified the translocation breakpoint at position 2158 in intron 6 in the 5′ MLL breakpoint cluster region (bcr). The karyotype in a case of treatment-related acute myeloid leukemia was normal, but Southern blot analysis showed a single MLL gene rearrangement. Panhandle PCR amplified the breakpoint at position 1493 in MLL intron 6. Screening of somatic cell hybrid and radiation hybrid DNAs by PCR and reverse transcriptase-PCR analysis of the leukemic cells indicated that panhandle PCR identified a fusion of MLL intron 6 with a previously uncharacterized sequence in MLL intron 1, consistent with a partial duplication. In both cases, the breakpoints in the MLL bcr were in Alu repeats, and there were Alu repeats in proximity to the breakpoints in the partner DNAs, suggesting that Alu sequences were relevant to these rearrangements. This study shows that panhandle PCR is an effective method for cloning MLL genomic breakpoints in treatment-related leukemias. Analysis of additional pediatric cases will determine whether breakpoint distribution deviates from the predilection for 3′ distribution in the bcr that has been found in adult cases.

Fluorescence correlation analysis of probe diffusion simplifies quantitative pathogen detection by PCR

A sensitive, labor-saving, and easily automatable nonradioactive procedure named APEX-FCS (amplified probe extension detected by fluorescence correlation spectroscopy) has been established to detect specific in vitro amplification of pathogen genomic sequences. As an example, Mycobacterium tuberculosis genomic DNA was subjected to PCR amplification with the Stoffel fragment of Thermus aquaticus DNA polymerase in the presence of nanomolar concentrations of a rhodamine-labeled probe (third primer), binding to the target in between the micromolar amplification primers. The probe becomes extended only when specific amplification occurs. Its low concentration avoids false-positives due to unspecific hybridization under PCR conditions. With increasing portion of extended probe molecules, the probe’s average translational diffusion properties gradually change over the course of the reaction, reflecting amplification kinetics. Following PCR, this change from a stage of high to a stage of low mobility can directly be monitored during a 30-s measurement using a fluorescence correlation spectroscopy device. Quantitation down to 10 target molecules in a background of 2.5 μg unspecific DNA without post-PCR probe manipulations could be achieved with different primer/probe combinations. The assay holds the promise to concurrently perform amplification, probe hybridization, and specific detection without opening the reaction chamber, if sealable foils are used.

Genomic analysis of human and mouse TCL1 loci reveals a complex of tightly clustered genes

TCL1 and TCL1b genes on human chromosome 14q23.1 are activated in T cell leukemias by translocations and inversions at 14q32.1, juxtaposing them to regulatory elements of T cell receptor genes. In this report we present the cloning, mapping, and expression analysis of the human and murine TCL1/Tcl1 locus. In addition to TCL1 and TCL1b, the human locus contains two additional genes, TCL1-neighboring genes (TNG) 1 and 2, encoding proteins of 141 and 110 aa, respectively. Both genes show no homology to any known genes, but their expression profiles are very similar to those of TCL1 and TCL1b. TNG1 and TNG2 also are activated in T cell leukemias with rearrangements at 14q32.1. To aid in the development of a mouse model we also have characterized the murine Tcl1 locus and found five genes homologous to human TCL1b. Tcl1b1–Tcl1b5 proteins range from 117 to 123 aa and are 65–80% similar, but they show only a 30–40% similarity to human TCL1b. All five mouse Tcl1b and murine Tcl1 mRNAs are abundant in mouse oocytes and two-cell embryos but rare in various adult tissues and lymphoid cell lines. These data suggest a similar or complementary function of these proteins in early embryogenesis.

A model system to study genomic imprinting of human genes

Somatic-cell hybrids have been shown to maintain the correct epigenetic chromatin states to study developmental globin gene expression as well as gene expression on the active and inactive X chromosomes. This suggests the potential use of somatic-cell hybrids containing either a maternal or a paternal human chromosome as a model system to study known imprinted genes and to identify as-yet-unknown imprinted genes. Testing gene expression by using reverse transcription followed by PCR, we show that functional imprints are maintained at four previously characterized 15q11–q13 loci in hybrids containing a single human chromosome 15 and at two chromosome 11p15 loci in hybrids containing a single chromosome 11. In contrast, three γ-aminobutyric acid type A receptor subunit genes in 15q12–q13 are nonimprinted. Furthermore, we have found that differential DNA methylation imprints at the SNRPN promoter and at a CpG island in 11p15 are also maintained in somatic-cell hybrids. Somatic-cell hybrids therefore are a valid and powerful system for studying known imprinted genes as well as for rapidly identifying new imprinted genes.

Genomic structure and ecdysone regulation of the prophenoloxidase 1 gene in the malaria vector Anopheles gambiae

Prophenoloxidase, a melanin-synthesizing enzyme, is considered to be an important arthropod immune protein. In mosquitoes, prophenoloxidase has been shown to be involved in refractory mechanisms against malaria parasites. In our study we used Anopheles gambiae, the most important human malaria vector, to characterize the first arthropod prophenoloxidase gene at the genomic level. The complete nucleotide sequence, including the immediate 5′ flanking sequence (−855 bp) of the prophenoloxidase 1 gene, was determined. The gene spans 10 kb and is composed of five exons and four introns coding for a 2.5-kb mRNA. In the 5′ flanking sequence, we found several putative regulatory motifs, two of which were identified as ecdysteroid regulatory elements. Electrophoretic mobility gel-shift assays and supershift assays demonstrated that the Aedes aegypti ecdysone receptor/Ultraspiracle nuclear receptor complex, and, seemingly, the endogenous Anopheles gambiae nuclear receptor complex, was able to bind one of the ecdysteroid response elements. Furthermore, 20-hydroxyecdysone stimulation was shown to up-regulate the transcription of the prophenoloxidase 1 gene in an A. gambiae cell line.