Mouse mammary tumor virus/v-Ha-ras transgene-induced mammary tumors exhibit strain-specific allelic loss on mouse chromosome 4

Hybrid mice carrying oncogenic transgenes afford powerful systems for investigating loss of heterozygosity (LOH) in tumors. Here, we apply this approach to a neoplasm of key importance in human medicine: mammary carcinoma. We performed a whole genome search for LOH using the mouse mammary tumor virus/v-Ha-ras mammary carcinoma model in female (FVB/N × Mus musculus castaneus)F1 mice. Mammary tumors developed as expected, as well as a few tumors of a second type (uterine leiomyosarcoma) not previously associated with this transgene. Genotyping of 94 anatomically independent tumors revealed high-frequency LOH (≈38%) for markers on chromosome 4. A marked allelic bias was observed, with M. musculus castaneus alleles almost exclusively being lost. No evidence of genomic imprinting effects was noted. These data point to the presence of a tumor suppressor gene(s) on mouse chromosome 4 involved in mammary carcinogenesis induced by mutant H-ras expression, and for which a significant functional difference may exist between the M. musculus castaneus and FVB/N alleles. Provisional subchromosomal localization of this gene, designated Loh-3, can be made to a distal segment having syntenic correspondence to human chromosome 1p; LOH in this latter region is observed in several human malignancies, including breast cancers. Evidence was also obtained for a possible second locus associated with LOH with less marked allele bias on proximal chromosome 4.

CD68+ cells of monocyte/macrophage lineage in the environment of AIDS-associated and classic-sporadic Kaposi sarcoma are singly or doubly infected with human herpesviruses 7 and 6B

Earlier studies have shown that Kaposi sarcomas contain cells infected with human herpesvirus (HHV) 6B, and in current studies we report that both AIDS-associated and classic-sporadic Kaposi sarcoma contain HHV-7 genome sequences detectable by PCR. To determine the distribution of HHV-7-infected cells relative to those infected with HHV-6, sections from paraffin-embedded tissues were allowed to react with antibodies to HHV-7 virion tegument phosphoprotein pp85 and to HHV-6B protein p101. The antibodies are specific for HHV-7 and HHV-6B, respectively, and they retained reactivity for antigens contained in formalin-fixed, paraffin-embedded tissue samples. We report that (i) HHV-7 pp85 was present in 9 of 32 AIDS-associated Kaposi sarcomas, and in 1 of 7 classical-sporadic HIV-negative Kaposi sarcomas; (ii) HHV-7 pp85 was detected primarily in cells bearing the CD68 marker characteristic of the monocyte/macrophage lineage present in or surrounding the Kaposi sarcoma lesions; and (iii) in a number of Kaposi sarcoma specimens, tumor-associated CD68+ monocytes/macrophages expressed simultaneously antigens from both HHV-7 and HHV-6B, and therefore appeared to be doubly infected with the two viruses. CD68+ monocytes/macrophages infected with HHV-7 were readily detectable in Kaposi sarcoma, but virtually absent from other normal or pathological tissues that harbor macrophages. Because all of the available data indicate that HHV-7 infects CD4+ T lymphocytes, these results suggest that the environment of the Kaposi sarcoma (i) attracts circulating peripheral lymphocytes and monocytes, triggers the replication of latent viruses, and thereby increases the local concentration of viruses, (ii) renders CD68+ monocytes/macrophages susceptible to infection with HHV-7, and (iii) the combination of both events enables double infections of cells with both HHV-6B and HHV-7.

Determining divergence times with a protein clock: Update and reevaluation

A recent study of the divergence times of the major groups of organisms as gauged by amino acid sequence comparison has been expanded and the data have been reanalyzed with a distance measure that corrects for both constraints on amino acid interchange and variation in substitution rate at different sites. Beyond that, the availability of complete genome sequences for several eubacteria and an archaebacterium has had a great impact on the interpretation of certain aspects of the data. Thus, the majority of the archaebacterial sequences are not consistent with currently accepted views of the Tree of Life which cluster the archaebacteria with eukaryotes. Instead, they are either outliers or mixed in with eubacterial orthologs. The simplest resolution of the problem is to postulate that many of these sequences were carried into eukaryotes by early eubacterial endosymbionts about 2 billion years ago, only very shortly after or even coincident with the divergence of eukaryotes and archaebacteria. The strong resemblances of these same enzymes among the major eubacterial groups suggest that the cyanobacteria and Gram-positive and Gram-negative eubacteria also diverged at about this same time, whereas the much greater differences between archaebacterial and eubacterial sequences indicate these two groups may have diverged between 3 and 4 billion years ago.

Gene content phylogeny of herpesviruses

Clusters of orthologous groups [COGs; Tatusov, R. L., Koonin, E. V. & Lipman, D. J. (1997) Science 278, 631–637] were identified for a set of 13 completely sequenced herpesviruses. Each COG represented a family of gene products conserved across several herpes genomes. These families were defined without using an arbitrary threshold criterion based on sequence similarity. The COG technique was modified so that variable stringency in COG construction was possible. High stringencies identify a core set of highly conserved genes. Varying COG stringency reveals differences in the degree of conservation between functional classes of genes. The COG data were used to construct whole-genome phylogenetic trees based on gene content. These trees agree well with trees based on other methods and are robust when tested by bootstrap analysis. The COG data also were used to construct a reciprocal tree that clustered genes with similar phylogenetic profiles. This clustering may give clues to genes with related functions or with related histories of acquisition and loss during herpesvirus evolution.

Building a dictionary for genomes: Identification of presumptive regulatory sites by statistical analysis

The availability of complete genome sequences and mRNA expression data for all genes creates new opportunities and challenges for identifying DNA sequence motifs that control gene expression. An algorithm, “MobyDick,” is presented that decomposes a set of DNA sequences into the most probable dictionary of motifs or words. This method is applicable to any set of DNA sequences: for example, all upstream regions in a genome or all genes expressed under certain conditions. Identification of words is based on a probabilistic segmentation model in which the significance of longer words is deduced from the frequency of shorter ones of various lengths, eliminating the need for a separate set of reference data to define probabilities. We have built a dictionary with 1,200 words for the 6,000 upstream regulatory regions in the yeast genome; the 500 most significant words (some with as few as 10 copies in all of the upstream regions) match 114 of 443 experimentally determined sites (a significance level of 18 standard deviations). When analyzing all of the genes up-regulated during sporulation as a group, we find many motifs in addition to the few previously identified by analyzing the subclusters individually to the expression subclusters. Applying MobyDick to the genes derepressed when the general repressor Tup1 is deleted, we find known as well as putative binding sites for its regulatory partners.

The Comprehensive Microbial Resource

One challenge presented by large-scale genome sequencing efforts is effective display of uniform information to the scientific community. The Comprehensive Microbial Resource (CMR) contains robust annotation of all complete microbial genomes and allows for a wide variety of data retrievals. The bacterial information has been placed on the Web at http://www.tigr.org/CMR for retrieval using standard web browsing technology. Retrievals can be based on protein properties such as molecular weight or hydrophobicity, GC-content, functional role assignments and taxonomy. The CMR also has special web-based tools to allow data mining using pre-run homology searches, whole genome dot-plots, batch downloading and traversal across genomes using a variety of datatypes.

Does recombination improve selection on codon usage? Lessons from nematode and fly complete genomes

Understanding the factors responsible for variations in mutation patterns and selection efficacy along chromosomes is a prerequisite for deciphering genome sequences. Population genetics models predict a positive correlation between the efficacy of selection at a given locus and the local rate of recombination because of Hill–Robertson effects. Codon usage is considered one of the most striking examples that support this prediction at the molecular level. In a wide range of species including Caenorhabditis elegans and Drosophila melanogaster, codon usage is essentially shaped by selection acting for translational efficiency. Codon usage bias correlates positively with recombination rate in Drosophila, apparently supporting the hypothesis that selection on codon usage is improved by recombination. Here we present an exhaustive analysis of codon usage in C. elegans and D. melanogaster complete genomes. We show that in both genomes there is a positive correlation between recombination rate and the frequency of optimal codons. However, we demonstrate that in both species, this effect is due to a mutational bias toward G and C bases in regions of high recombination rate, possibly as a direct consequence of the recombination process. The correlation between codon usage bias and recombination rate in these species appears to be essentially determined by recombination-dependent mutational patterns, rather than selective effects. This result highlights that it is necessary to take into account the mutagenic effect of recombination to understand the evolutionary role and impact of recombination.

Functional genomics

Complete genome sequences are providing a framework to allow the investigation of biological processes by the use of comprehensive approaches. Genome analysis also is having a dramatic impact on medicine through its identification of genes and mutations involved in disease and the elucidation of entire microbial gene sets. Studies of the sequences of model organisms, such as that of the nematode worm Caenorhabditis elegans, are providing extraordinary insights into development and differentiation that aid the study of these processes in humans. The field of functional genomics seeks to devise and apply technologies that take advantage of the growing body of sequence information to analyze the full complement of genes and proteins encoded by an organism.

Testing ancient RNA–protein interactions

The past decade in molecular biology has seen remarkable advances in the study of the origin and early evolution of life. The mathematical tools for analyzing DNA and protein sequences, coupled with the availability of complete microbial genome sequences, provide insight almost as far back as the age of the nucleic acids themselves. Experimental evolution in the laboratory and especially in vitro evolution of RNA provide insight into a hypothetical world where RNA, or a close relative, may have debuted as a primary functional and informational molecule. The ability to isolate new functional RNAs from random sequences now ultimately makes the world of possible primitive chemical interactions accessible even when the molecules or reactions are no longer present in modern species. Thus we can at last form direct experimental tests of specific models for the origin of RNA–protein associations, such as those that influenced the genetic code. This marks a turning point for probing the origin and early history of life at the molecular level.

Comparative organization of cattle chromosome 5 revealed by comparative mapping by annotation and sequence similarity and radiation hybrid mapping

A whole genome cattle-hamster radiation hybrid cell panel was used to construct a map of 54 markers located on bovine chromosome 5 (BTA5). Of the 54 markers, 34 are microsatellites selected from the cattle linkage map and 20 are genes. Among the 20 mapped genes, 10 are new assignments that were made by using the comparative mapping by annotation and sequence similarity strategy. A LOD-3 radiation hybrid framework map consisting of 21 markers was constructed. The relatively low retention frequency of markers on this chromosome (19%) prevented unambiguous ordering of the other 33 markers. The length of the map is 398.7 cR, corresponding to a ratio of ≈2.8 cR5,000/cM. Type I genes were binned for comparison of gene order among cattle, humans, and mice. Multiple internal rearrangements within conserved syntenic groups were apparent upon comparison of gene order on BTA5 and HSA12 and HSA22. A similarly high number of rearrangements were observed between BTA5 and MMU6, MMU10, and MMU15. The detailed comparative map of BTA5 should facilitate identification of genes affecting economically important traits that have been mapped to this chromosome and should contribute to our understanding of mammalian chromosome evolution.

Regulation of the Mycobacterium tuberculosis hypoxic response gene encoding α-crystallin

Unlike many pathogens that are overtly toxic to their hosts, the primary virulence determinant of Mycobacterium tuberculosis appears to be its ability to persist for years or decades within humans in a clinically latent state. Since early in the 20th century latency has been linked to hypoxic conditions within the host, but the response of M. tuberculosis to a hypoxic signal remains poorly characterized. The M. tuberculosis α-crystallin (acr) gene is powerfully and rapidly induced at reduced oxygen tensions, providing us with a means to identify regulators of the hypoxic response. Using a whole genome microarray, we identified >100 genes whose expression is rapidly altered by defined hypoxic conditions. Numerous genes involved in biosynthesis and aerobic metabolism are repressed, whereas a high proportion of the induced genes have no known function. Among the induced genes is an apparent operon that includes the putative two-component response regulator pair Rv3133c/Rv3132c. When we interrupted expression of this operon by targeted disruption of the upstream gene Rv3134c, the hypoxic regulation of acr was eliminated. These results suggest a possible role for Rv3132c/3133c/3134c in mycobacterial latency.

Use of cDNA subtraction and RNA interference screens in combination reveals genes required for germ-line development in Caenorhabditis elegans

Caenorhabditis elegans is an ideal organism for the study of the molecular basis of fundamental biological processes such as germ-line development, especially because of availability of the whole genome sequence and applicability of the RNA interference (RNAi) technique. To identify genes involved in germ-line development, we produced subtracted cDNA pools either enriched for or deprived of the cDNAs from germ-line tissues. We then performed differential hybridization on the high-density cDNA grid, on which about 7,600 nonoverlapping expressed sequence tag (EST) clones were spotted, to identify a set of genes specifically expressed in the germ line. One hundred and sixty-eight clones were then tested with the RNAi technique. Of these, 15 clones showed sterility with a variety of defects in germ-line development. Seven of them led to the production of unfertilized eggs, because of defects in spermatogenesis (4 clones), or defects in the oocytes (3 clones). The other 8 clones led to failure of oogenesis. These failures were caused by germ-line proliferation defect (Glp phenotype), meiotic arrest, and defects in sperm–oocyte switch (Mog phenotype) among others. These results demonstrate the efficacy of the screening strategy using the EST library combined with the RNAi technique in C. elegans.

Intron insertion facilitates amplification of cloned virus cDNA in Escherichia coli while biological activity is reestablished after transcription in vivo.

Insertion of introns into cloned cDNA of two isolates of the plant potyvirus pea seedborne mosaic virus facilitated plasmid amplification in Escherichia coli. Multiple stop codons in the inserted introns interrupted the open reading frame of the virus cDNA, thereby terminating undesired translation of virus proteins in E. coli. Plasmids containing the full-length virus sequences, placed under control of the cauliflower mosaic virus 35S promoter and the nopaline synthase termination signal, were stable and easy to amplify in E. coli if one or more introns were inserted into the virus sequence. These plasmids were infectious when inoculated mechanically onto Pisum sativum leaves. Examination of the cDNA-derived viruses confirmed that intron splicing of in vivo transcribed pre-mRNA had occurred as predicted, reestablishing the virus genome sequences. Symptom development and virus accumulation of the cDNA derived viruses and parental viruses were identical. It is proposed that intron insertion can be used to facilitate manipulation and amplification of cloned DNA fragments that are unstable in, or toxic to, E. coli. When transcribed in vivo in eukaryotic cells, the introns will be eliminated from the sequence and will not interfere with further analysis of protein expression or virus infection.

Repair of x-ray-induced DNA double-strand breaks in specific Not I restriction fragments in human fibroblasts: joining of correct and incorrect ends.

An assay that allows measurement of absolute induction frequencies for DNA double-strand breaks (dsbs) in defined regions of the genome and that quantitates rejoining of correct DNA ends has been used to study repair of dsbs in normal human fibroblasts after x-irradiation. The approach involves hybridization of single-copy DNA probes to Not I restriction fragments separated according to size by pulsed-field gel electrophoresis. Induction of dsbs is quantitated from the decrease in the intensity of the hybridizing restriction fragment and an accumulation of a smear below the band. Rejoining of dsbs results in reconstitution of the intact restriction fragment only if correct DNA ends are joined. By comparing results from this technique with results from a conventional electrophoresis assay that detects all rejoining events, it is possible to quantitate the misrejoining frequency. Three Not I fragments on the long arm of chromosome 21 were investigated with regard to dsb induction, yielding an identical induction rate of 5.8 X 10(-3) break per megabase pair per Gy. Correct dsb rejoining was measured for two of these Not I fragments after initial doses of 80 and 160 Gy. The misrejoining frequency was about 25% for both fragments and was independent of dose. This result appears to be representative for the whole genome as shown by analysis of the entire Not I fragment distribution. The correct rejoining events primarily occurred within the first 2 h, while the misrejoining kinetics included a much slower component, with about half of the events occurring between 2 and 24 h. These misrejoining kinetics are similar to those previously reported for production of exchange aberrations in interphase chromosomes.

Ancient mtDNA sequences in the human nuclear genome: A potential source of errors in identifying pathogenic mutations

Nuclear-localized mtDNA pseudogenes might explain a recent report describing a heteroplasmic mtDNA molecule containing five linked missense mutations dispersed over the contiguous mtDNA CO1 and CO2 genes in Alzheimer’s disease (AD) patients. To test this hypothesis, we have used the PCR primers utilized in the original report to amplify CO1 and CO2 sequences from two independent ρ° (mtDNA-less) cell lines. CO1 and CO2 sequences amplified from both of the ρ° cells, demonstrating that these sequences are also present in the human nuclear DNA. The nuclear pseudogene CO1 and CO2 sequences were then tested for each of the five “AD” missense mutations by restriction endonuclease site variant assays. All five mutations were found in the nuclear CO1 and CO2 PCR products from ρ° cells, but none were found in the PCR products obtained from cells with normal mtDNA. Moreover, when the overlapping nuclear CO1 and CO2 PCR products were cloned and sequenced, all five missense mutations were found, as well as a linked synonymous mutation. Unlike the findings in the original report, an additional 32 base substitutions were found, including two in adjacent tRNAs and a two base pair deletion in the CO2 gene. Phylogenetic analysis of the nuclear CO1 and CO2 sequences revealed that they diverged from modern human mtDNAs early in hominid evolution about 770,000 years before present. These data would be consistent with the interpretation that the missense mutations proposed to cause AD may be the product of ancient mtDNA variants preserved as nuclear pseudogenes.