Alignment editing and identification of consensus secondary structures for nucleic acid sequences: interactive use of dot matrix representations.

We present a computer-aided approach for identifying and aligning consensus secondary structure within a set of functionally related oligonucleotide sequences aligned by sequence. The method relies on visualization of secondary structure using a generalization of the dot matrix representation appropriate for consensus sequence data sets. An interactive computer program implementing such a visualization of consensus structure has been developed. The program allows for alignment editing, data and display filtering and various modes of base pair representation, including co-variation. The utility of this approach is demonstrated with four sample data sets derived from in vitro selection experiments and one data set comprising tRNA sequences.     

A high speed, high capacity homology matrix: zooming through SV40 and polyoma.

We present a new homology matrix program which owes its basic conception to the two-dimensional dot matrices previously described (1,2), but has important improvements and new features. It scores sequence homology over an adjustable range and plots the scores which are above an operator-determined filtration level. Its powerful noise-filtration system, capacity for compression without much loss of information, and speed of execution make this program a valuable tool in the analysis of homologies, internal direct repeats and reverse repeats, including palindromic sequences. The properties of the program are exemplified by analysis of SV40 and polyoma DNA sequences.     

Matrix program to analyze primary structure homology

A FORTRAN program to analyze homology of letter strings (nucleotide or amino acid sequences) and to display the result in the form of a dot matrix is presented. The program is generally usable, user-friendly and has a number of options (filtering, "fudging," i.e., consideration of groups of homologous residues, and screening, i.e., display of only particular groups of residues) which greatly potentiate its analytical power.     

A simple image analysis program for quantitative dot assays

Summary Dot assays are versatile, and are widely used for determining antigens and proteins. Because they require expensive equipment to be quantitative, often only qualitative dot results are reported. However, because the dot pattern is so regular, a simple image analysis program can determine mean dot grey levels, while handling irregular dot outlines, radial variations in colour within a dot, and varying or noisy sheet backgrounds. We describe a dot analysis program that runs on PCs under Windows, and permits quantitative dot assays to be run with inexpensive grayscale scanner input. The program is available from us as source and executable files. We present demonstration results for antigen and protein determinations.     

A Macintosh program for the versatile generation of random nucleic acid sequences and their structural analysis

The program ‘MacStAn’ for the Apple Macintosh generatesrandom sequences and can analyze their tendency to form secondarystructure or translation products as well as their mono-, di-and trinucleotide composition. Generation of random sequencesis versatile in that one can (i) predefine the G + C content,maximal base repetitions and constant regions; (ii) preset theentire dinucleotide composition; or (iii) shuffle an existingsequence. The program constitutes an integrated package witha graphical user interface, fill-featured editing, saving, printing,text import and export, dot plot and sequence alignment.     

Screen for stage-specific expression genes between tail bud stage and heartbeat beginning stage in embryogenesis of gynogenetic silver crucian carp

A systemic study was initiated to identify stage-specific expression genes in fish embryogenesis by using suppression subtractive hybridization (SSH) technique. In this study, we presented a preliminary result on screen for stage-specific expression genes between tail bud stage (TBS) and heartbeat beginning stage (HBS) in gynogenetic silver crucian carp (Carassius auratus gibelio). Two SSH plasmid libraries specific for TBS embryos and HBS embryos were constructed, and stage-specific expression genes were screened between the two stages. 1963 TBS positive clones and 2466 HBS positive clones were sampled to PCR amplification, and 1373 TBS and 1809 HBS PCR positive clones were selected to carry out dot blots. 169 TBS dot blot positive clones and 272 HBS dot blot positive clones were sequenced. Searching GenBank by using these nucleotide sequences indicated that most of the TBS dot blot positive clones could not be found homologous sequences in the database, while known genes were mainly detected from HBS d     

SIMS: computation of a smooth invariant molecular surface.

SIMS, a new method of calculating a smooth invariant molecular dot surface, is presented. The SIMS method generates the smooth molecular surface by rolling two probe spheres. A solvent probe sphere is rolled over the molecule and produces a Richards-Connolly molecular surface (MS), which envelops the solvent-excluded volume of the molecule. In deep crevices, Connolly's method of calculating the MS has two deficiencies. First, it produces self-intersecting parts of the molecular surface, which must be removed to obtain the correct MS. Second, the correct MS is not smooth, i.e., the direction of the normal vector of the MS is not continuous, and some points of the MS are singular. We present an exact method for removing self-intersecting parts and smoothing the singular regions of the MS. The singular MS is smoothed by rolling a smoothing probe sphere over the inward side of the singular MS. The MS in the vicinity of singularities is replaced with the reentrant surface of the smoothing probe sphere. The smoothing method does not disturb the topology of a singular MS, and the smooth MS is a better approximation of the dielectric border between high dielectric solvent and the low dielectric molecular interior. The SIMS method generates a smooth molecular dot surface, which has a quasi-uniform dot distribution in two orthogonal directions on the molecular surface, which is invariant with molecular rotation and stable under changes in the molecular conformation, and which can be used in a variety of implicit methods of modeling solvent effects. The SIMS program is faster than the Connolly MS program, and in a matter of seconds generates a smooth dot MS of a 200-residue protein. The program is available from the authors on request (see http:@femto.med.unc.edu/SIMS).     

反向斑点杂交法检测HBV基本核心启动子区A1762T／G1764A突变的实验研究

目的建立一种基于尼龙膜的反向斑点杂交法,用于检测乙型肝炎病毒（HBV）基本核心启动子区（BCP）A1762T／G1764A突变。方法根据我国HBV主要流行的基因型为B和C,从GenBank上查出4种HBVBCP序列。利用在线工具ClustalW进行比对,针对该突变位点设计引物和检测探针。探针经合成和修饰后点在带正电的尼龙膜上。将反向斑点杂交法结合地高辛检测试剂盒用于检测A1762T／G1764A突变,以测序法确定该区域序列的标本为检测对象。结果反向斑点杂交法分别检测5例A1762／G1764病毒株、2例T1762／G1764病毒株、5例A1762／A1764病毒株和4例T1762／A1764病毒株,结果与测序完全相同。结论应用本方法可以快速、准确地HBV相关的热点突变。     

Consensus sequence L/PKSSLL mimics crucial epitope on Loop III of Taiwan cobra cardiotoxin

Phage display is effective in screening peptides that mimic venom’s neutralizing epitopes. A phage display cyclized heptapeptide library (C7C library) was panned with purified divalent antivenin IgG, which neutralizes Naja naja atra venom (NAV) and Bungarus multicinctus venom (BMV). The selected heptapeptide sequences were aligned with known protein sequences of NAV and BMV in GenBank. One of the four consensus sequences, L/PKSSLL, mimicked the crucial epitope on Loop III of Taiwan cobra cardiotoxin that is associated with the venom’s lethal potency. In dot blot analysis, several clones showed varying reactivities for NAV monovalent antivenin and lesser cross-reactions with BMV monovalent antivenin. The KSSLLRN-carrying phage occurred four times in selected clones and showed the strongest reactivity to NAV monovalent antivenin. Furthermore, the QDSLLPS-carrying phage also presented significant dot blot signal, indicating that the SLL sequence shared by these two clones may be a crucial antibody-binding site.     

MATLAB 7.X生物信息工具箱的应用——序列比对(二)

序列比对是基因序列分析中的一项重要工作.本文以人和鼠的基因为对象,介绍MATLAB 7.X生物信息工具箱中的序列比对方法,内容包括从数据库获取序列信息,查找序列的开放阅读框,将核苷酸序列转换为氨基酸序列,绘制比较两氨基酸序列的散点图,用Needleman-Wunsch算法和Smith-Waterman算法进行比对,以及计算两序列的同一性.     

Adplot: detection and visualization of repetitive patterns in complete genomes

MOTIVATION: Repetitive DNA sequences are abundant in genomes and efficient mining of significant repeats is important as the first step of repetitive sequence research. Although many computational tools for the purpose, either automatic or visualization ones, have been developed, detection and analysis of approximate repeats are still non-trivial task. RESULTS: Auto Dot PLOT (Adplot), a dotplot-like repetitive pattern visualization program with a window filtering based on iid Bernoulli trials, is developed and applied to yeast chromosomes and human T cell receptor locus sequence. Typical examples found in yeast chromosomes 1 and 10 and a tandem repeat of periods longer than 10,000 bp in human T cell receptor locus are presented. A complex structure composed of both direct and palindromic repeats found in yeast chromosome 10 is also visualized as specific dot pattern. Computational time measured by a Pentium 3 PC for each yeast auto chromosome with a standard parameter setting is linearly scaled and below 10 s per one chromosome, indicating efficiency of the program. From the examples, it is shown that Adplot can visualize approximate local repeat structures and give us a diagnosis power for inferring a duplicational history of repeats. AVAILABILITY: Adplot can be obtained by an e-mail request.     

A fast word search algorithm for the representation of sequence similarity in genomic DNA.

Representation of sequence similarity by dot matrix plots is a method widely used for comparing biological sequences. The user is presented with an overall view of similarity between two sequences. Computation of this plot has been reconsidered here. An improvement is proposed through the preprocessing of the data into an automation recognizing the word structure of a sequence. The main advantage of this approach is to systematically eliminate the repetitions during word comparison. Simple heuristics are also considered to greatly speed up pattern matching. As a result, large sequences are handled very efficiently. This is illustrated by a comparison of large genomic DNA. The algorithm has been implemented in an interactive application on a microcomputer.     

Analyses of the splenic mRNA expressed by rabbits of different immunoglobulin kappa-light chain allotypes: conserved sequences in the 3' untranslated region and allotype-specific probes

The allelism of the structural genes for the complex rabbit b allotypes of immunoglobulin kappa-light chains has been questioned because of observations of unexpected phenotypic expression of "latent" allotypes. We find that the coding sequences of the b4 and b5 "alleles" are only 80% homologous for the last 60 nucleotides but there is a high degree of homology (96%) in the 3' untranslated region (3'DT). The high conservation of 3' DT region sequences enabled us to detect kappa-light chain mRNAs from rabbits of different genetic types (b4, b5, b9 and bbas) on northern blots and dot blots. We can distinguish mRNA encoding b9 and b5 allotypes on dot blots with b5 fragment-probes of known sequence and detect mRNA produced by unstimulated cultured splenic lymphocytes. Analyses of mRNA from cultured cells manipulated to enhance mRNA synthesis and production of unexpected or "latent" b allotypes can now be conducted.     

Identification and Analysis of Genes Present in Leptospira interrogans serovar lai but Absent in L. biflexa serovar monvalerio

Leptospirosis is a globally important zoonotic diseasecaused by the pathogenic species of the spirochete genus,Leptospira including L. interrogans, L. kirschneri, L.noguchii, L. borgpetersenii, L. santarosai, L. weilii andetc. [1]. Pathogenic leptospires …     

Identification and analysis of genes present in Leptospira interrogans serovar lai but absent in L. biflexa serovar monvalerio

Genes present in virulent bacterial strains but absent in avirulent close relatives can be of great biologic and clinical interest. This project aimed to identify strain specific DNA sequences of Leptospira interrogens serovar lai, which is absent in the saprophytic L. biflexa serovar monvalerio, via suppression subtractive hybridization with the former as the tester while the latter as the driver. The mixture of PCR amplified DNA fragments from two subtractive hybridization experiments were cloned into pMD 18-T vector and the positive clones were identified by dot blotting against the chromosome DNA of the two strains individually. After DNA sequencing and analysis, the distribution of these genomic fragment sequences in a panel of pathogenic and nonpathogenic leptospires was investigated employing dot blot analysis. Among the 188 positive clones randomly chosen, 24 contained the tester strain specific genomic regions, of which, 5 were non-coding fragments while the others contained 23 distinct protein coding sequences. Besides 9 genes encoding functional proteins, 12 genes encode unknown proteins and the rest two genes encode proteins with recognizable domain structures, one for a putative leucine-rich repeats (LRR) family protein while the other as an outer-membrane protein. Our experiment results indicated that suppression subtractive hybridization is effective for screening specific DNA sequences between two leptospiral strains, and some of these sequences might be responsible for virulence determination. Further analysis of these DNA sequences will provide important information on the pathogenesis of Leptospira.     

Recognition of 3' -processing sites of human mRNA precursors

We have developed a computer program POLYAH and an algorithmfor the identification of 3'-processing sites of human mRNAprecursors. The algorithm is based on a linear discriminantfunction (LDF) trained to discriminate real poly(A) signal regionsfrom the other regions of human genes possessing the AATAAAsequence which is most likely nonfunctional. As the parametersof LDF, various significant contextual characteristics of sequencessurrounding AATAAA signals were used. An accuracy of methodhas been estimated on a set of 131 poly(A) regions and 1466regions of human genes having the AATAAA sequence. When thethreshold was set to predict 86% of poly(A) regions correctly,specificity of 51% and correlation coefficient of 0.62 had beenachieved. The precision of this approach is better than forthe other methods and has been tested on a larger data set.POLYAH can be used through World Wide Web (at Gene-Finder Homepage: URL http: //dot.imgen.bcm.tmc.edu: 9331/gene-finder/ gf.html)or by sending files with uncharacterized human sequences tothe University of Houston or Weizmann Institute of Science e-mailservers.     

The Mirror-Type Internal Symmetry in the Primary Structure of Proteins: Detection and Functional Role (Review)

This review summarized the data obtained by the author in studies on internal symmetry of the mirror type in primary structures of proteins. The methods for detection of symmetric segments in amino acid sequences are analyzed: (1) the method based on analysis of sequences of roots of amino acid codons; (2) the dot matrix method; (3) the method of internal symmetry scanning. The results of studies of internal symmetry in enzymes and signaling proteins are presented. The probable role of the internal symmetry in the structural-functional organization of proteins is discussed.     

Cue-triggered activity replay in human early visual cortex

The recall of learned temporal sequences by a visual cue is an important form of experience-based neural plasticity. Here we observed such reactivation in awake human visual cortex using intracranial recording. After repeated exposure to a moving dot, a flash of the dot was able to trigger neural reactivation in the downstream receptive field along the motion path. This effect was observed only when the cue appeared near the receptive field. The estimated traveling speed was faster compared to the activation induced by the real motion. We suggest a range-limited, time-compressed reactivation as a result of repeated visual exposure in awake human visual cortex.     

LINKER: a program to generate linker sequences for fusion proteins

The construction of functional fusion proteins often requires a linker sequence that adopts an extended conformation to allow for maximal flexibility. Linker sequences are generally selected based on intuition. Without a reliable selection criterion, the design of such linkers is often difficult, particularly in situations where longer linker sequences are required. Here we describe a program called LINKER which can automatically generate a set of linker sequences that are known to adopt extended conformations as determined by X-ray crystallography and NMR. The only required input to the program is the desired linker sequence length. The program is specifically designed to assist in fusion protein construction. A number of optional input parameters have been incorporated so that users are able to enhance sequence selection based on specific applications. The program output simply contains a set of sequences with a specified length. This program should be a useful tool in both the biotechnology industry and biomedical research. It can be accessed through the Web page http://www.fccc. edu/research/labs/feng/linker.html.