Optimized homology searches of the gene and protein sequence data banks

Identification of a cytochrome P450 gene by reverse transcription--PCR using degenerate primers containing inosine

A cytochrome P450-like gene, tentatively named P450CMEF, was amplified by a mixed oligonucleotide-primed amplification of cDNA from C3H mouse embryo fibroblast cells, designated 10T1/2, that had been treated with 7,12-dimethylbenz[a]anthracene (DMBA) or benz[a]anthracene (BA). A set of inosine-containing degenerate primers that were targeted to two conserved regions of known cytochrome P450 cDNAs were used. One primer was coded for the well-described and conserved heme-binding region of P450 enzymes, and the second was designed based upon other considerations of homology among P450 molecules. One of the four PCR-amplified cDNA products hybridized to two major RNA bands, 4.2 and 5.3 kb, that were induced by DMBA or BA. The amino acid sequence of the fragment deduced from the base-sequence data indicate that the amplified cDNA has a 50-55% identity with the cytochrome P450 subfamily 1A. The induction of P450CMEF mRNA preceded the induction of aryl hydrocarbon hydroxylase activity after DMBA or BA treatment, suggesting that the product of P450CMEF is involved in the metabolism of these polycyclic aromatic hydrocarbons in 10T1/2 cells. From the partial sequence of the cDNA identified by this procedure, we propose that P450CMEF is a member of the P450 superfamily, possibly in a subfamily of family 1, that is induced in 10T1/2 cells by DMBA and BA. This method should be useful in identifying additional P450 genes and genes in other gene families.

A mouse Y chromosome pseudogene is related to human ubiquitin activating enzyme E1

A 2041 bp DNA fragment isolated from the Sxr (sex reversed) region of the mouse Y Chromosome (Chr) was sequenced and characterized. The sequence, pY8/b, contains four exons that are highly similar to 525 contiguous bases from the cDNA of human ubiquitin activating enzyme E1. Two of the exons contain stop codons, indicating that pY8/b is not part of a functional gene. Sequences related to pY8/b were amplified from the Y Chr of the inbred mouse strain, C57BL/6J. These sequences may be portions of the recently discovered functional equivalent of pY8/b. Despite a high degree of similarity with the human E1 gene, the functional equivalent of pY8/b is not the mouse E1 gene, because unlike E1, the functional equivalent of pY8/b is expressed in a tissue-specific manner. These data are discussed with respect to theory on the evolution of the mammalian Y Chr, and in particular, to the prediction that functional genes on the Y Chr have a male-specific function.

Sequence and analysis of a portion of the genomes of Shope fibroma virus and malignant rabbit fibroma virus that is important for viral replication in lymphocytes

The 10.7-kb BamHI "C" restriction fragment of malignant rabbit fibroma virus (MV) contains genes that are important for its immunosuppressive activity. When this fragment is transferred to a related avirulent leporipoxvirus, Shope fibroma virus (SFV), recombinant viruses show clinical features characteristic of MV: they replicate in lymphocytes and alter immune function in vitro, induce disseminated tumors in recipient rabbits, and are immunosuppressive in vivo. The 10.7-kb BamHI "C" restriction fragment of MV was sequenced in its entirety. Its DNA sequence and the 14 ORF's derived from analyzing this sequence are discussed. Analysis of known open reading frames to which the ORF's from MV's Bam "C" fragment show homology permits us to identify some MV ORF's showing high degrees of similarity to known and postulated proteins produced by vaccinia virus. Functions for some of these vaccinia proteins are known, while functions for others are hypothetical or unknown. Further analysis of genetic determinants of MV's virulence has indicated that two overlapping restriction subfragments of the BamHI "C" fragment can transfer MV's virulent behavior to SFV. The 0.7-kb region in which these two subfragments overlap includes the C-terminus of MV orf C-7 and the N terminus of MV orf C-8. These correspond to the C- and N-termini, respectively, of SFV orf's D-9 and D-10 and to vaccinia orf's D-6 (early transcription factor) and D-7 (subunit of RNA polymerase). We sequenced the region of SFV's BamHI "D" fragment in this area and illustrate here the comparative sequences of this portion of SFV's genome and orf's. On the basis of comparisons between MV, SFV, and vaccinia in this area we discuss the potential significance of these observations.

RNA secondary structure analysis of the packaging signal for Moloney murine leukemia virus

We have examined the cis-acting RNA packaging signal (psi) from Moloney murine leukemia virus using a combination of chemical and primary sequence analysis techniques. For our chemical analyses, we used dimethyl sulfate, kethoxal, and 1-cyclohexyl-3-(2-morpholinoethyl)-carbodiimide metho-p-toluene sulfonate as probes for RNA secondary structure. The structural information obtained from these studies was used to constrain computer algorithms for prediction of RNA secondary structure. In addition, we generated and analyzed a phylogenetic comparison of homologous sequences from related retroviruses. From these data, we have developed two models for the RNA secondary structure of the packaging signal psi. Both of these models suggest the presence of secondary structure elements in a region of the psi RNA known to be required for function.

The use of primers from highly conserved pol regions to identify uncharacterized retroviruses by the polymerase chain reaction

Two degenerate oligonucleotide primers derived from regions of pol conserved among retroviruses have been synthesized. Polymerase chain reactions utilizing these primers amplify a 135-bp pol fragment in every retrovirus DNA tested to date. The polymerase chain reaction has been linked to a reverse transcriptase step so that a pol-specific DNA fragment can be obtained from a moderate amount of a purified retrovirus or viral RNA. The identity of an unknown retrovirus can be determined by sequencing of the amplified fragment following molecular cloning. This procedure was tested on an unidentified (non-HIV) retrovirus expressed by a B-cell lymphoma line obtained from an AIDS patient. Our PCR assay identified the retrovirus as being highly similar to Mason-Pfizer monkey virus (MPMV) and simian retrovirus 1, which are closely related immunosuppressive type D viruses that cause simian AIDS.

Use of homology domains in sequence similarity detection

We have found the detection of homology domains using a nonlinear similarity score and the DD algorithm to be a useful approach for identifying similarity between sequences and evaluating potential homology. There are several reasons for the success of the method. (i) Homology domains are identified by a rigorous method that guarantees they will be locally optimal. (ii) The relative significance of different homology domains can be directly compared using the nonlinear similarity score. (iii) Different cost matrices can be used in the calculation of the similarity score. (iv) Relatively long sequences can be compared in a single pass as storage requirements are proportional to the shorter of the two sequences being compared. (v) The method has proved to be very sensitive in practice. (vi) The boundaries of authentic regions of homology are accurately identified. (vii) The information required to define a homology domain (its location, size, similarity score, etc.) can be stored in a compact data structure, facilitating the sharing of homology domain data among different software tools. (viii) The method can be applied to similarity searches of the nucleotide and protein sequence data banks. These properties make the identification of homology domains for studying sequence similarity a useful companion to other accepted methods, such as dynamic programming based analyses.

Multicopy suppression of the cdc24 budding defect in yeast by CDC42 and three newly identified genes including the ras-related gene RSR1

Genes CDC24, CDC42, and CDC43 are required for the establishment of cell polarity and the localization of secretion in Saccharomyces cerevisiae; mutants defective in these genes fail to form buds and display isotropic expansion of the cell surface. To identify other genes that may be involved in these processes, we screened yeast genomic DNA libraries for heterologous genes that, when overexpressed from a plasmid, can suppress a temperature-sensitive cdc24 mutation. We identified four such genes. One of these proved to be CDC42, which has previously been shown to be a member of the rho (ras-homologous) family of genes, and a second is a newly identified ras-related gene that we named RSR1. RSR1 maps between CDC62 and ADE3 on the right arm of chromosome VII; its predicted product is approximately 50% identical to other proteins in the ras family. Deletion of RSR1 is nonlethal but disrupts the normal pattern of bud site selection. Although both CDC42 and RSR1 can suppress cdc24 and both appear to encode GTP-binding proteins, these genes do not themselves appear to be functionally interchangeable. However, one of the other genes that was isolated by virtue of its ability to suppress cdc24 can also suppress cdc42. This gene, named MSB1, maps between ADE9 and HIS3 on the right arm of chromosome XV.

Two new members of the OmpR superfamily detected by homology to a sensor-binding core domain

The OmpR superfamily includes proteins that act as transcriptional regulators of operons that respond to environmental stimuli. A homologous domain near the N-terminus, termed a sensor-binding core domain, is thought to play a role in recognition of a signal transduction protein. We have identified two previously unrecognized members of this regulator family of proteins: a 23.8-kd protein transcribed from the uvrC transcription unit and the PgtA gene product, which is a phosphoglycerate transport regulatory protein. The sensor-binding core domain is also present in four proteins that regulate bacterial sporulation and chemotaxis. The 23.8-kd protein also has sequence similarity to elongation factor Tu and two regulatory proteins: HtpR, the heat-shock regulatory protein, and TraJ, a regulator of expression of genes involved in conjugation. There is a 77-amino acid region near the C-terminus of the 23.8-kd protein that has 30% similarity with a 28.1-kd protein coded for by an open reading frame 5' to the reading frame of the 23.8-kd protein in the uvrC transcription unit. Genetic distance analysis of amino acid sequences of proteins with a sensor-binding core domain suggests that the 23.8-kd protein and the chemotaxis regulatory proteins are distantly related to the other regulatory proteins in the OmpR superfamily.

DNA----DNA, and DNA----RNA----protein: orchestration by a single complex operon

In Escherichia coli, the workhorse of molecular biology, a single operon is involved in the replication, transcription and translation of genetic information. This operon is controlled in a complex manner involving multiple cis-acting regulatory sequences and trans-acting regulatory proteins. It interacts with global regulatory networks by mechanisms which are presently being dissected.

Identification of a conserved sequence in the non-coding regions of many human genes

We have analyzed a sequence of approximately 70 base pairs (bp) that shows a high degree of similarity to sequences present in the non-coding regions of a number of human and other mammalian genes. The sequence was discovered in a fragment of human genomic DNA adjacent to an integrated hepatitis B virus genome in cells derived from human hepatocellular carcinoma tissue. When one of the viral flanking sequences was compared to nucleotide sequences in GenBank, more than thirty human genes were identified that contained a similar sequence in their non-coding regions. The sequence element was usually found once or twice in a gene, either in an intron or in the 5' or 3' flanking regions. It did not share any similarities with known short interspersed nucleotide elements (SINEs) or presently known gene regulatory elements. This element was highly conserved at the same position within the corresponding human and mouse genes for myoglobin and N-myc, indicating evolutionary conservation and possible functional importance. Preliminary DNase I footprinting data suggested that the element or its adjacent sequences may bind nuclear factors to generate specific DNase I hypersensitive sites. The size, structure, and evolutionary conservation of this sequence indicates that it is distinct from other types of short interspersed repetitive elements. It is possible that the element may have a cis-acting functional role in the genome.

The vitamin D receptor: a primitive steroid receptor related to thyroid hormone receptor

We have previously reported the cloning and sequencing of both the chicken and human vitamin D3 receptor cDNAs. A comparison of their deduced amino acid sequence with that of the other classic steroid hormone receptors and the receptor for thyroid hormone indicates that there are two regions of conservation between these molecules. The first is a 70 amino acid, cysteine-rich sequence (C1), the second region (C2) is a 62 amino acid region located towards the carboxyl terminus of the proteins. In other systems the former has been identified as a region responsible for DNA binding activity, whereas the latter represents the NH2-terminal boundary of the hormone binding domain. We present here evidence utilizing eucaryotic expression of cDNA encoding the hVDR C1 domain, followed by a DNA cellulose chromatography assay, which confirms that the DNA binding activity resides in this region of the receptor for vitamin D3. Additionally, the vitamin D3 receptor contains a 60 amino acid portion at its carboxyl terminus (C3) which exhibits homology with the receptor for thyroid hormone. Conservation in this region of the molecule is found only between homologous or closely related receptors. This indicates a relationship between the vitamin D3 receptor and the receptor for thyroid hormone and may suggest that they evolved from a single primordial gene.

Genetic relatedness of human DNA polymerase beta and terminal deoxynucleotidyltransferase

The Protein Identification Resource (PIR) protein sequence data bank was searched for sequence similarity between known proteins and human DNA polymerase beta (Pol beta) or human terminal deoxynucleotidyltransferase (TdT). Pol beta and TdT were found to exhibit amino acid sequence similarity only with each other and not with any other of the 4750 entries in release 12.0 of the PIR data bank. Optimal amino acid sequence alignment of the entire 39-kDa Pol beta polypeptide with the C-terminal two thirds of TdT revealed 24% identical aa residues and 21% conservative aa substitutions. The Monte Carlo score of 12.6 for the entire aligned sequences indicates highly significant aa sequence homology. The hydropathicity profiles of the aligned aa sequences were remarkably similar throughout, suggesting structural similarity of the polypeptides. The most significant regions of homology are aa residues 39-224 and 311-333 of Pol beta vs. aa residues 191-374 and 484-506 of TdT. In addition, weaker homology was seen between a large portion of the 'nonessential' N-terminal end of TdT (aa residues 33-130) and the first region of strong homology between the two proteins (aa residues 31-128 of Pol beta and aa residues 183-280 of TdT), suggestive of genetic duplication within the ancestral gene. On the basis of nucleotide differences between conserved regions of Pol beta and TdT genes (aligned according to optimally aligned aa sequences) it was estimated that Pol beta and TdT diverged on the order of 250 million years ago, corresponding roughly to a time before radiation of mammals and birds.