The RDP-II (Ribosomal Database Project)

The Ribosomal Database Project (RDP-II), previously described by Maidak et al. [Nucleic Acids Res. (2000), 28, 173–174], continued during the past year to add new rRNA sequences to the aligned data and to improve the analysis commands. Release 8.0 (June 1, 2000) consisted of 16 277 aligned prokaryotic small subunit (SSU) rRNA sequences while the number of eukaryotic and mitochondrial SSU rRNA sequences in aligned form remained at 2055 and 1503, respectively. The number of prokaryotic SSU rRNA sequences more than doubled from the previous release 14 months earlier, and ~75% are longer than 899 bp. An RDP-II mirror site in Japan is now available (http://wdcm.nig.ac.jp/RDP/html/index.html). RDP-II provides aligned and annotated rRNA sequences, derived phylogenetic trees and taxonomic hierarchies, and analysis services through its WWW server (http://rdp.cme.msu.edu/). Analysis services include rRNA probe checking, approximate phylogenetic placement of user sequences, screening user sequences for possible chimeric rRNA sequences, automated alignment, production of similarity matrices and services to plan and analyze terminal restriction fragment polymorphism experiments. The RDP-II email address for questions and comments has been changed from curator@cme.msu.edu to rdpstaff@msu.edu.

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.

Cleavage motifs of the yeast 20S proteasome β subunits deduced from digests of enolase 1

The 436-amino acid protein enolase 1 from yeast was degraded in vitro by purified wild-type and mutant yeast 20S proteasome particles. Analysis of the cleavage products at different times revealed a processive degradation mechanism and a length distribution of fragments ranging from 3 to 25 amino acids with an average length of 7 to 8 amino acids. Surprisingly, the average fragment length was very similar between wild-type and mutant 20S proteasomes with reduced numbers of active sites. This implies that the fragment length is not influenced by the distance between the active sites, as previously postulated. A detailed analysis of the cleavages also allowed the identification of certain amino acid characteristics in positions flanking the cleavage site that guide the selection of the P1 residues by the three active β subunits. Because yeast and mammalian proteasomes are highly homologous, similar cleavage motifs might be used by mammalian proteasomes. Therefore, our data provide a basis for predicting proteasomal degradation products from which peptides are sampled by major histocompatibility complex class I molecules for presentation to cytotoxic T cells.

Formation of stable adducts and absence of depurinating DNA adducts in cells and DNA treated with the potent carcinogen dibenzo[a,l]pyrene or its diol epoxides

Polycyclic aromatic hydrocarbons (PAH) are widespread environmental contaminants, and some are potent carcinogens in rodents. Carcinogenic PAH are activated in cells to metabolites that react with DNA to form stable covalent DNA adducts. It has been proposed [Cavalieri, E. L. & Roger, E. G. (1995) Xenobiotica 25, 677–688] that unstable DNA adducts are also formed and that apurinic sites in the DNA resulting from unstable PAH adducts play a key role in the initiation of cancer. The potent carcinogen dibenzo[a,l]pyrene (DB[a,l]P) is activated in cells to (+)-syn- and (−)-anti-DB[a,l]P-11,12-diol-13,14-epoxide (DB[a,l]PDE), which have been shown to form stable adducts with DNA. To evaluate the importance of unstable PAH adducts, we compared stable adduct formation to apurinic site formation. Stable DB[a,l]PDE adducts were determined by 33P-postlabeling and HPLC. To measure apurinic sites they were converted to strand breaks, and these were monitored by examining the integrity of a particular restriction fragment of the dihydrofolate reductase gene. The method easily detected apurinic sites resulting from methylation by treatment of cells or DNA with dimethyl sulfate or from reaction of DNA with DB[a,l]P in the presence of horseradish peroxidase. We estimate the method could detect 0.1 apurinic site in the 14-kb fragment examined. However, apurinic sites were below our limit of detection in DNA treated directly with (+)-syn- or (−)-anti-DB[a,l]PDE or in DNA from Chinese hamster ovary B11 cells so treated, although in these samples the frequency of stable adducts ranged from 3 to 10 per 14 kb. We also treated the human mammary carcinoma cell line MCF-7 with DB[a,l]P and again could not detect significant amounts of unstable adducts. These results indicate that the proportion of stable adducts formed by DB[a,l]P activated in cells and its diol epoxides is greater than 99% and suggest a predominant role for stable DNA adducts in the carcinogenic activity of DB[a,l]P.

Adaptive evolution via a major gene effect: Paedomorphosis in the Mexican axolotl

Although adaptive evolution is thought to depend primarily on mutations of small effect, major gene effects may underlie many of the important differences observed among species in nature. The Mexican axolotl (Ambystoma mexicanum) has a derived mode of development that is characterized by metamorphic failure (paedomorphosis), an adaptation for an entirely aquatic life cycle. By using an interspecific crossing design and genetic linkage analysis, a major quantitative trait locus for expression of metamorphosis was identified in a local map of amplified fragment length polymorphisms. These data are consistent with a major gene hypothesis for the evolution of paedomorphosis in A. mexicanum.

GenMapDB: a database of mapped human BAC clones

GenMapDB (http://genomics.med.upenn.edu/genmapdb) is a repository of human bacterial artificial chromosome (BAC) clones mapped by our laboratory to sequence-tagged site markers. Currently, GenMapDB contains over 3000 mapped clones that span 19 chromosomes, chromosomes 2, 4, 5, 9–22, X and Y. This database provides positional information about human BAC clones from the RPCI-11 human male BAC library. It also contains restriction fragment analysis data and end sequences of the clones. GenMapDB is freely available to the public. The main purpose of GenMapDB is to organize the mapping data and to allow the research community to search for mapped BAC clones that can be used in gene mapping studies and chromosomal mutation analysis projects.

Heat effects on DNA repair after ionising radiation: hyperthermia commonly increases the number of non-repaired double-strand breaks and structural rearrangements

After ionising radiation double-strand breaks (dsb) are lethal if not repaired or misrepaired. Cell killing is greatly enhanced by hyperthermia and it is questioned here whether heat not only affects dsb repair capacity but also fidelity in a chromosomal context. dsb repair experiments were designed so as to mainly score non-homologous end joining, while homologous recombination was largely precluded. Human male G0 fibroblasts were either preheated (45°C, 20 min) or not before X-irradiation. dsb induction and repair were measured by conventional gel electrophoresis and an assay combining restriction digestion using a rare cutting enzyme (NotI) and Southern hybridisation, which detects large chromosomal rearrangements (>100 kb). dsb induction rate in an X-chromosomal NotI fragment was 4.8 × 10–3 dsb/Gy/Mb. Similar values were found for the genome overall and also when cells were preheated. After 50 Gy, fibroblasts were competent to largely restore the original restriction fragment size. Five per cent of dsb remained non-rejoined and 14% were misrejoined. Correct restitution of restriction fragments occurred preferably during the first hour but continued at a slow rate for 12–16 h. In addition, dsb appeared to misrejoin throughout the entire repair period. After hyperthermia the fractions of non-rejoined and misrejoined dsb were similarly increased to 13 and 51%, respectively. It is suggested that heat increases the probability of dsb being incorrectly rejoined but it is not likely to interfere with one dsb repair pathway in particular.

Two genes abrogate the inhibition of murine hepatocarcinogenesis by ovarian hormones.

Hormonal and genetic factors strongly influence the susceptibility of inbred mice to hepatocarcinogenesis. Female C57BR/cdJ (BR) mice are extremely susceptible to liver tumor induction relative to other strains because they are genetically insensitive to the inhibition of hepatocarcinogenesis by ovarian hormones. To determine the genetic basis for the sensitivity of BR mice relative to resistant C57BL/6J (B6) mice, we treated 12-day-old B6BRF1 x B6 and B6BRF1 x B6BRF1 (F2) animals with N,N-diethylnitrosamine (0.1 micromol/g of body weight) and enumerated liver tumors at 32 weeks of age in males and at 50 weeks in females. Genomic DNA samples from backcross and F2 mice were analyzed for 70 informative simple sequence length polymorphism markers. Genetic markers on chromosome 17 (D17Mit21) and chromosome 1 (D1Mit33) cosegregated with high tumor multiplicity in both sexes. Together, these loci [designated Hcf1 and Hcf2 (Hepatocarcinogenesis in females), respectively] account for virtually all of the difference in sensitivity between BR and B6 mice. The Hcf1 locus accounts for a majority of the higher susceptibility of BR mice of both sexes. Backcross female mice heterozygous at both loci (33 +/- 23 tumors per mouse) and at Hcf1 only (17 +/- 18) were 15- and 8-fold more sensitive, respectively, than mice homozygous for the B6 alleles at Hcf1 and Hcf2 (2.2 +/- 3.9). In backcross male mice, the double heterozygotes (35 +/- 22) and Hcf1 heterozygotes (28 +/- 12) were 5.4- and 4.3-fold more sensitive than mice homozygous for B6 alleles at both loci (6.5 +/- 5.4).

Phylogenetic assessment of length variation at a microsatellite locus

Sixty-six haplotypes at a locus containing a simple dinucleotide (CA)n microsatellite repeat were isolated by PCR–single-strand conformational polymorphism from populations of the horseshoe crab Limulus polyphemus. These haplotypes were sequenced to assess nucleotide variation directly. Thirty-four distinct sequences (alleles) were identified in a region 570 bp long that included the microsatellite motif. In the repeat region itself, CA-number varied in integer values from 5 to 11 across alleles, except that a (CA)8 class was not observed. Differences among alleles were due also to polymorphisms at 22 sites in regions immediately flanking the microsatellite repeats. Nucleotide substitutions in these regions were used to estimate phylogenetic relationships among alleles, and the gene phylogeny was used to trace the evolution of length variation and CA repeat numbers. A low correlation between size variation and genealogical relationships among alleles suggests that absolute fragment size (as normally scored in microsatellite assays) is an unreliable indicator of historical affinities among alleles. This finding on the molecular fine structure of microsatellite variation suggests the need for caution in the use of repeat counts at microsatellite loci as secure indicators of allelic relationships.

Plasmon analyses of Triticum (wheat) and Aegilops: PCR–single-strand conformational polymorphism (PCR-SSCP) analyses of organellar DNAs

To investigate phylogenetic relationships among plasmons in Triticum and Aegilops, PCR–single-strand conformational polymorphism (PCR-SSCP) analyses were made of 14.0-kb chloroplast (ct) and 13.7-kb mitochondrial (mt)DNA regions that were isolated from 46 alloplasmic wheat lines and one euplasmic line. These plasmons represent 31 species of the two genera. The ct and mtDNA regions included 10 and 9 structural genes, respectively. A total of 177 bands were detected, of which 40.6% were variable. The proportion of variable bands in ctDNA (51.1%) was higher than that of mtDNA (28.9%). The phylogenetic trees of plasmons, derived by two different models, indicate a common picture of plasmon divergence in the two genera and suggest three major groups of plasmons (Einkorn, Triticum, and Aegilops). Because of uniparental plasmon transmission, the maternal parents of all but one polyploid species were identified. Only one Aegilops species, Ae. speltoides, was included in the Triticum group, suggesting that this species is the plasmon and B and G genome donor of all polyploid wheats. ctDNA variations were more intimately correlated with vegetative characters, whereas mtDNA variations were more closely correlated with reproductive characters. Plasmon divergence among the diploids of the two genera largely paralleled genome divergence. The relative times of origin of the polyploid species were inferred from genetic distances from their putative maternal parents.

PCR candidate region mismatch scanning: adaptation to quantitative, high-throughput genotyping

Linkage and association analyses were performed to identify loci affecting disease susceptibility by scoring previously characterized sequence variations such as microsatellites and single nucleotide polymorphisms. Lack of markers in regions of interest, as well as difficulty in adapting various methods to high-throughput settings, often limits the effectiveness of the analyses. We have adapted the Escherichia coli mismatch detection system, employing the factors MutS, MutL and MutH, for use in PCR-based, automated, high-throughput genotyping and mutation detection of genomic DNA. Optimal sensitivity and signal-to-noise ratios were obtained in a straightforward fashion because the detection reaction proved to be principally dependent upon monovalent cation concentration and MutL concentration. Quantitative relationships of the optimal values of these parameters with length of the DNA test fragment were demonstrated, in support of the translocation model for the mechanism of action of these enzymes, rather than the molecular switch model. Thus, rapid, sequence-independent optimization was possible for each new genomic target region. Other factors potentially limiting the flexibility of mismatch scanning, such as positioning of dam recognition sites within the target fragment, have also been investigated. We developed several strategies, which can be easily adapted to automation, for limiting the analysis to intersample heteroduplexes. Thus, the principal barriers to the use of this methodology, which we have designated PCR candidate region mismatch scanning, in cost-effective, high-throughput settings have been removed.

Enhancing the catalytic repertoire of nucleic acids: a systematic study of linker length and rigidity

The incorporation of potentially catalytic groups in DNA is of interest for the in vitro selection of novel deoxyribozymes. A series of 10 C5-modified analogues of 2′-deoxyuridine triphosphate have been synthesised that possess side chains of differing flexibility and bearing a primary amino or imidazole functionality. For each series of nucleotide analogues differing degrees of flexibility of the C5 side chain was achieved through the use of alkynyl, alkenyl and alkyl moieties. The imidazole function was conjugated to these C5-amino-modified nucleotides using either imidazole 4-acetic acid or imidazole 4-acrylic acid (urocanic acid). The substrate properties of the nucleotides (fully replacing dTTP) with Taq polymerase during PCR have been investigated in order to evaluate their potential applications for in vitro selection experiments. 5-(3-Aminopropynyl)dUTP and 5-(E-3-aminopropenyl)dUTP and their imidazole 4-acetic acid- and urocanic acid-modified conjugates were found to be substrates. In contrast, C5-amino-modified dUTPs with alkane or Z-alkene linkers and their corresponding conjugates were not substrates. The incorporation of these analogues during PCR has been confirmed by inhibition of restriction enzyme digestion using XbaI and by mass spectrometry of the PCR products.

Amplification of the full-length hepatitis A virus genome by long reverse transcription-PCR and transcription of infectious RNA directly from the amplicon.

The genetic study of RNA viruses is greatly facilitated by the availability of infectious cDNA clones. However, their construction has often been difficult. While exploring ways to simplify the construction of infectious clones, we have successfully modified and applied the newly described technique of "long PCR" to the synthesis of a full-length DNA amplicon from the RNA of a cytopathogenic mutant (HM 175/24a) of the hepatitis A virus (HAV). Primers were synthesized to match the two extremities of the HAV genome. The antisense primer, homologous to the 3' end, was used in both the reverse transcription (RT) and the PCR steps. With these primers we reproducibly obtained a full-length amplicon of approximately 7.5 kb. Further, since we engineered a T7 promoter in the sense primer, RNA could be transcribed directly from the amplicon with T7 RNA polymerase. Following transfection of cultured fetal rhesus kidney cells with the transcription mixture containing both the HAV cDNA and the transcribed RNA, replicating HAV was detected by immunofluorescence microscopy and, following passage to other cell cultures, by focus formation. The recovered virus displayed the cytopathic effect and large plaque phenotype typical of the original virus; this result highlights the fidelity of the modified long reverse transcription-PCR procedure and demonstrates the potential of this method for providing cDNAs of viral genomes and simplifying the construction of infectious clones.

Heterogeneity of primer extension products in asymmetric PCR is due both to cleavage by a structure-specific exo/endonuclease activity of DNA polymerases and to premature stops.

In PCR, DNA polymerases from thermophilic bacteria catalyze the extension of primers annealed to templates as well as the structure-specific cleavage of the products of primer extension. Here we show that cleavage by Thermus aquaticus and Thermus thermophilus DNA polymerases can be precise and substantial: it occurs at the base of the stem-loop structure assumed by the single strand products of primer extension using as template a common genetic element, the promoter-operator of the Escherichia coli lactose operon, and may involve up to 30% of the products. The cleavage is independent of primer, template, and triphosphates, is dependent on substrate length and temperature, requires free ends and Mg2+, and is absent in DNA polymerases lacking the 5'-->3' exonuclease, such as the Stoffel fragment and the T7 DNA polymerase. Heterogeneity of the extension products results also from premature detachment of the enzyme approaching the 5' end of the template.

Nonsense mutation in the phosphofructokinase muscle subunit gene associated with retention of intron 10 in one of the isolated transcripts in Ashkenazi Jewish patients with Tarui disease.

Mutations in the human phosphofructokinase muscle subunit gene (PFKM) are known to cause myopathy classified as glycogenosis type VII (Tarui disease). Previously described molecular defects include base substitutions altering encoded amino acids or resulting in abnormal splicing. We report a mutation resulting in phosphofructokinase deficiency in three patients from an Ashkenazi Jewish family. Using a reverse transcription PCR assay, PFKM subunit transcripts differing by length were detected in skeletal muscle tissue of all three affected subjects. In the longer transcript, an insertion of 252 nucleotides totally homologous to the structure of the 10th intron of the PFKM gene was found separating exon 10 from exon 11. In addition, two single base transitions were identified by direct sequencing: [exon 6; codon 95; CGA (Arg) to TGA (stop)] and [exon 7; codon 172; ACC (Thr) to ACT (Thr)] in either transcript. Single-stranded conformational polymorphism and restriction enzyme analyses confirmed the presence of these point substitutions in genomic DNA and strongly suggested homozygosity for the pathogenic allele. The nonsense mutation at codon 95 appeared solely responsible for the phenotype in these patients, further expanding genetic heterogeneity of Tarui disease. Transcripts with and without intron 10 arising from identical mutant alleles probably resulted from differential pre-mRNA processing and may represent a novel message from the PFKM gene.