Computerized polymorphic marker identification: Experimental validation and a predicted human polymorphism catalog

A computational system for the prediction of polymorphic loci directly and efficiently from human genomic sequence was developed and verified. A suite of programs, collectively called pompous (polymorphic marker prediction of ubiquitous simple sequences) detects tandem repeats ranging from dinucleotides up to 250 mers, scores them according to predicted level of polymorphism, and designs appropriate flanking primers for PCR amplification. This approach was validated on an approximately 750-kilobase region of human chromosome 3p21.3, involved in lung and breast carcinoma homozygous deletions. Target DNA from 36 paired B lymphoblastoid and lung cancer lines was amplified and allelotyped for 33 loci predicted by pompous to be variable in repeat size. We found that among those 36 predominately Caucasian individuals 22 of the 33 (67%) predicted loci were polymorphic with an average heterozygosity of 0.42. Allele loss in this region was found in 27/36 (75%) of the tumor lines using these markers. pompous provides the genetic researcher with an additional tool for the rapid and efficient identification of polymorphic markers, and through a World Wide Web site, investigators can use pompous to identify polymorphic markers for their research. A catalog of 13,261 potential polymorphic markers and associated primer sets has been created from the analysis of 141,779,504 base pairs of human genomic sequence in GenBank. This data is available on our Web site (pompous.swmed.edu) and will be updated periodically as GenBank is expanded and algorithm accuracy is improved.

trEST, trGEN and Hits: access to databases of predicted protein sequences

High throughput genome (HTG) and expressed sequence tag (EST) sequences are currently the most abundant nucleotide sequence classes in the public database. The large volume, high degree of fragmentation and lack of gene structure annotations prevent efficient and effective searches of HTG and EST data for protein sequence homologies by standard search methods. Here, we briefly describe three newly developed resources that should make discovery of interesting genes in these sequence classes easier in the future, especially to biologists not having access to a powerful local bioinformatics environment. trEST and trGEN are regularly regenerated databases of hypothetical protein sequences predicted from EST and HTG sequences, respectively. Hits is a web-based data retrieval and analysis system providing access to precomputed matches between protein sequences (including sequences from trEST and trGEN) and patterns and profiles from Prosite and Pfam. The three resources can be accessed via the Hits home page (http://hits.isb-sib.ch).

Predicted highly expressed and putative alien genes of Deinococcus radiodurans and implications for resistance to ionizing radiation damage

Predicted highly expressed (PHX) and putative alien genes determined by codon usages are characterized in the genome of Deinococcus radiodurans (strain R1). Deinococcus radiodurans (DEIRA) can survive very high doses of ionizing radiation that are lethal to virtually all other organisms. It has been argued that DEIRA is endowed with enhanced repair systems that provide protection and stability. However, predicted expression levels of DNA repair proteins with the exception of RecA tend to be low and do not distinguish DEIRA from other prokaryotes. In this paper, the capability of DEIRA to resist extreme doses of ionizing and UV radiation is attributed to an unusually high number of PHX chaperone/degradation, protease, and detoxification genes. Explicitly, compared with all current complete prokaryotic genomes, DEIRA contains the greatest number of PHX detoxification and protease proteins. Other sources of environmental protection against severe conditions of UV radiation, desiccation, and thermal effects for DEIRA are the several S-layer (surface structure) PHX proteins. The top PHX gene of DEIRA is the multifunctional tricarboxylic acid (TCA) gene aconitase, which, apart from its role in respiration, also alerts the cell to oxidative damage.

A Determinant of Substrate Specificity Predicted from the Acyl-Acyl Carrier Protein Desaturase of Developing Cat's Claw Seed1

Cat's claw (Doxantha unguis-cati L.) vine accumulates nearly 80% palmitoleic acid (16:1Δ9) plus cis-vaccenic acid (18:1Δ11) in its seed oil. To characterize the biosynthetic origin of these unusual fatty acids, cDNAs for acyl-acyl carrier protein (acyl-ACP) desaturases were isolated from developing cat's claw seeds. The predominant acyl-ACP desaturase cDNA identified encoded a polypeptide that is closely related to the stearoyl (Δ9–18:0)-ACP desaturase from castor (Ricinis communis L.) and other species. Upon expression in Escherichia coli, the cat's claw polypeptide functioned as a Δ9 acyl-ACP desaturase but displayed a distinct substrate specificity for palmitate (16:0)-ACP rather than stearate (18:0)-ACP. Comparison of the predicted amino acid sequence of the cat's claw enzyme with that of the castor Δ9–18:0-ACP desaturase suggested that a single amino acid substitution (L118W) might account in large part for the differences in substrate specificity between the two desaturases. Consistent with this prediction, conversion of leucine-118 to tryptophan in the mature castor Δ9–18:0-ACP desaturase resulted in an 80-fold increase in the relative specificity of this enzyme for 16:0-ACP. The alteration in substrate specificity observed in the L118W mutant is in agreement with a crystallographic model of the proposed substrate-binding pocket of the castor Δ9–18:0-ACP desaturase.

Alternatively spliced mRNAs predicted to yield frame-shift proteins and stable intron 1 RNAs of the herpes simplex virus 1 regulatory gene alpha 0 accumulate in the cytoplasm of infected cells.

The infected cell protein no. 0 (ICP0), the product of the alpha 0 gene, and an important herpes simplex virus 1 regulatory protein is encoded by three exons. We report that intron 1 forms a family of four stable nonpolyadenylylated cytoplasmic RNAs sharing a common 5' end but differing in 3' ends. The 5' and 3' ends correspond to the accepted splice donor and four splice acceptor sites within the mapped intron domain. The most distant splice acceptor site yields the mRNA encoding the 775-aa protein known as ICP0. The mRNAs resulting from the use of alternative splice acceptor sites were also present in the cytoplasm of infected cells and would be predicted to encode proteins of 152 (ICP0-B), 87 (ICP0-C), and 90 (ICP0-D) amino acids, respectively. Both the stability of the alpha 0 mRNA and the utilization of at least one splice acceptor site was regulated by ICP22 and or US1.5 protein inasmuch as cells infected with a mutant from which these genes had been deleted accumulated smaller amounts of alpha 0 mRNA than would be predicted from the amounts of accumulated intron RNAs. In addition, one splice acceptor site was at best underutilized. These results indicate that both the splicing pattern and longevity of alpha 0 mRNA are regulated. These and other recent examples indicate that herpes simplex virus 1 regulates its own gene expression and that of the infected cells through control of mRNA splicing and longevity.

Human myosin VIIA responsible for the Usher 1B syndrome: a predicted membrane-associated motor protein expressed in developing sensory epithelia.

The gene encoding human myosin VIIA is responsible for Usher syndrome type III (USH1B), a disease which associates profound congenital sensorineural deafness, vestibular dysfunction, and retinitis pigmentosa. The reconstituted cDNA sequence presented here predicts a 2215 amino acid protein with a typical unconventional myosin structure. This protein is expected to dimerize into a two-headed molecule. The C terminus of its tail shares homology with the membrane-binding domain of the band 4.1 protein superfamily. The gene consists of 48 coding exons. It encodes several alternatively spliced forms. In situ hybridization analysis in human embryos demonstrates that the myosin VIIA gene is expressed in the pigment epithelium and the photoreceptor cells of the retina, thus indicating that both cell types may be involved in the USH1B retinal degenerative process. In addition, the gene is expressed in the human embryonic cochlear and vestibular neuroepithelia. We suggest that deafness and vestibular dysfunction in USH1B patients result from a defect in the morphogenesis of the inner ear sensory cell stereocilia.

Types of Listeria monocytogenes predicted by the positions of EcoRI cleavage sites relative to ribosomal RNA sequences.

By using taxonomic characters derived from EcoRI restriction endonuclease digestion of genomic DNA and hybridization with a labeled rRNA operon from Escherichia coli, a polymorphic structure of Listeria monocytogenes, characterized by fragments with different frequencies of occurrence, was observed. This structure was expanded by creating predicted patterns through a recursive process of observation, expectation, prediction, and assessment of completeness. This process was applied, in turn, to normalized strain patterns, fragment bands, and positions of EcoRI recognition sites relative to rRNA regions. Analysis of 1346 strains provided observed patterns, fragment sizes, and their frequencies of occurrence in the patterns. Fragment size statistics led to the creation of unobserved combinations of bands, predicted pattern types. The observed fragment bands revealed positions of EcoRI sites relative to rRNA sequences. Each EcoRI site had a frequency of occurrence, and unobserved fragment sizes were postulated on the basis of knowing the restriction site locations. The result of the recursion process applied to the components of the strain data was an extended classification with observed and predicted members.

Proprotein convertases in amphioxus: predicted structure and expression of proteases SPC2 and SPC3.

SPC2 and SPC3 are two members of a family of subtilisin-related proteases which play essential roles in the processing of prohormones into their mature forms in the pancreatic B cell and many other neuroendocrine cells. To investigate the phylogenetic origins and evolutionary functions of SPC2 and SPC3 we have identified and cloned cDNAs encoding these enzymes from amphioxus (Branchiostoma californiensis), a primitive chordate. The amino acid sequence of preproSPC2 contains 689 aa and is 71% identical to human SPC2. In contrast, amphioxus prproSPC3 consists of 774 aa and exhibits 55% identity to human SPC3. These results suggest that the primary structure of SPC2 has been more highly conserved during evolution than that of SPC3. To further investigate the function(s) of SPC2 and SPC3 in amphioxus, we have determined the regional expression of these genes by using a reverse transcriptase-linked polymerase chain reaction (RT-PCR) assay. Whole amphioxus was dissected longitudinally into four equal-length segments and RNA was extracted. Using RT-PCR to simultaneously amplify SPC2 and SPC3 DNA fragments, we found that the cranial region (section 1) expressed equal amounts of SPC2 and SPC3 mRNAs, whereas in the caudal region (section 4) the SPC2-to-SPC3 ratio was 5:1. In the mid-body sections 2 and 3 the SPC2-to-SPC3 ratio was 1:5. By RT-PCR we also determined that amphioxus ILP, a homologue of mammalian insulin/insulin-like growth factor, was expressed predominately in section 3. These results suggest that the relative levels of SPC2 and SPC3 mRNAs are specifically regulated in various amphioxus tissues. Furthermore, the ubiquitous expression of these mRNAs in the organism indicates that they are involved in the processing of other precursor proteins in addition to proILP.