Optimization of the Method for Determining the Motility Characteristics of Fish spermatozoa Using ImageJ Software and Excel Macros

Journal of Ichthyology - A method for computer-based determination of the motility characteristics of fish spermatozoa using ImageJ software and a number of original Excel macros that help...     

Molecular beacon sequence design algorithm

A method based on Web-based tools is presented to design optimally functioning molecular beacons. Molecular beacons, fluorogenic hybridization probes, are a powerful tool for the rapid and specific detection of a particular nucleic acid sequence. However, their synthesis costs can be considerable. Since molecular beacon performance is based on its sequence, it is imperative to rationally design an optimal sequence before synthesis. The algorithm presented here uses simple Microsoft Excel formulas and macros to rank candidate sequences. This analysis is carried out using mfold structural predictions along with other free Web-based tools. For smaller laboratories where molecular beacons are not the focus of research, the public domain algorithm described here may be usefully employed to aid in molecular beacon design.     

MAD: a suite of tools for microarray data management and processing

SUMMARY: Microarray data management and processing (MAD) is a set of Windows integrated software for microarray analysis. It consists of a relational database for data storage with many user-interfaces for data manipulation, several text file parsers and Microsoft Excel macros for automation of data processing, and a generator to produce text files that are ready for cluster analysis. AVAILABILITY: Executable is available free of charge on http://pompous.swmed.edu. The source code is also available upon request.     

MASIA: recognition of common patterns and properties in multiple aligned protein sequences

SUMMARY: MASIA is a software tool for pattern recognition in multiple aligned protein sequences. MASIA converts a sequence to a properties matrix that can be scanned in both vertical and horizontal steps. Consistent patterns are recognized based on the statistical significance of their occurrence. Preset macros can be altered on-line to seek any combination of amino acid properties or sequence characteristics. MASIA output can be used directly by our programs to predict the 3D structure of proteins. AVAILABILITY: Access MASIA at http://www.scsb.utmb.edu/masia/ma sia.html.     

DNA array analysis in a Microsoft Windows environment.

Microsoft Windows-based computers have evolved to the point that they provide sufficient computational and visualization power for robust analysis of DNA array data. In fact, smaller laboratories might prefer to carry out some or all of their analyses and visualization in a Windows environment, rather than alternative platforms such as UNIX. We have developed a series of manually executed macros written in Visual Basic for Microsoft Excel spreadsheets, that allows for rapid and comprehensive gene expression data analysis. The first macro assigns gene names to spots on the DNA array and normalizes individual hybridizations by expressing the signal intensity for each gene as a percentage of the sum of all gene intensities. The second macro streamlines statistical consideration of the confidence in individual gene measurements for sets of experimental replicates by calculating probability values with the Student's t test. The third macro introduces a threshold value, calculates expression ratios between experimental conditions, and calculates the standard deviation of the mean of the log ratio values. Selected columns of data are copied by a fourth macro to create a processed data set suitable for entry into a Microsoft Access database. An Access database structure is described that allows simple queries across multiple experiments and export of data into third-party data visualization software packages. These analysis tools can be used in their present form by others working with commercial E. coli membrane arrays, or they may be adapted for use with other systems. The Excel spreadsheets with embedded Visual Basic macros and detailed instructions for their use are available at http://www.ou.edu/microarray.     

MarC-V: a spreadsheet-based tool for analysis, normalization, and visualization of single cDNA microarray experiments

The comprehensive analysis and visualization of data extracted from cDNA microarrays can be a time-consuming and error-prone process that becomes increasingly tedious with increased number of gene elements on a particular microarray. With the increasingly large number of gene elements on today's microarrays, analysis tools must be developed to meet this challenge. Here, we present MarC-V, a Microsoft Excel spreadsheet tool with Visual Basic macros to automate much of the visualization and calculation involved in the analysis process while providing the familiarity and flexibility of Excel. Automated features of this tool include (i) lower-bound thresholding, (ii) data normalization, (iii) generation of ratio frequency distribution plots, (iv) generation of scatter plots color-coded by expression level, (v) ratio scoring based on intensity measurements, (vi) filtering of data based on expression level or specific gene interests, and (vii) exporting data for subsequent multi-array analysis. MarC-V also has an importing function included for GenePix results (GPR) raw data files.     

Useful Microsoft Word Macros for molecular biologists and protein chemists

Biologists today make extensive use of word processing programs for the production of research reports, literature reviews and grant proposals. Frequently, such programs become the default platform for viewing and the later publication of protein and nucleic acid sequence data. Thus, researchers often switch between their word processor and more specialized programs designed to analyze protein and nucleic acid sequences. It would be more convenient to perform these simple sequence analyses using the word processor without switching to another program. The focus here is on the use of the Visual Basic programming language, which is built into all recent versions of Microsoft Word to generate surprisingly complex and useful macros that can conveniently analyze several important features of protein and nucleic acid sequences. The standard Word interface can also be easily modified to display and run these macros from a pull-down menu. Several examples of this approach are provided.     

FiRe and microarrays: a fast answer to burning questions

FiRe is a user-friendly Excel macro designed to survey microarray data rapidly. This software interactively assembles data from different experiments and produces lists of candidate genes according to patterns of gene expression. Furthermore, macros bundled with FiRe can compare lists of genes, merge information from different spreadsheets, link candidates to information available from web-based databases, and produce heat-maps for easy visualization of microarray data. FiRe is freely available at http://www.unifr.ch/plantbio/FiRe/main.html .     

Genotype transposer: automated genotype manipulation for linkage disequilibrium analysis

SUMMARY: The purpose of this work is to provide the modern molecular geneticist with tools to perform more efficient and more accurate analysis of the genotype data they produce. By using Microsoft Excel macros written in Visual Basic, we can translate genotype data into a form readable by the versatile software 'Arlequin', read the Arlequin output, calculate statistics of linkage disequilibrium, and put the results in a format for viewing with the software 'GOLD'. AVAILABILITY: The software is available by FTP at: ftp://xcsg.iarc.fr/cox/Genotype_Transposer/. SUPPLEMENTARY INFORMATION: Detailed instruction and examples are available at: ftp://xcsg.iarc.fr/cox/Genotype&_Transposer/. Arlequin is available at: http://lgb.unige.ch/arlequin/. GOLD is available at: http://www.well.ox.ac.uk/asthma/GOLD/.     

CEDIT: a C interface and macro facility for protein sequence alignment editing in colour with Microsoft Word 5.0 for PCs

CEDIT, a C interface and macro facility that provides for thecolour editing of protein sequence alignments (up to 2000 sequences,5000 residues each) using Microsoft Word 5.0 for PCs is presented.CEDIT uses the ability of MS-Word 5.0 to display letters withthe desired colour to easily identify conservative homologiesacross the sequences. A glossary file with useful macros forthe sequence editing is provided, along with several utilitiesprograms for error checking, estimating sequence similaritiesand homology significance. CEDIT has a menu interface and acontext sensitive help.     

Sequence ordinations: a multivariate analysis approach to analysing large sequence data sets

Ordination is a powerful method for analysing complex data setsbut has been largely ignored in sequence analysis. This papershows how to use principal coordinates analysis to find low–dimensionalrepresentations of distance matrices derived from aligned setsof sequences. The method takes a matrix of Euclidean distancesbetween all pairs of sequence and finds a coordinate space wherethe distances are exactly preserved The main problem is to finda measure of distance between aligned sequences that is Euclidean.The simplest distance function is the square root of the percentagedifference (as measured by identities) between two sequences,where one ignores any positions in the alignment where thereis a gap in any sequence. If one does not ignore positions witha gap, the distances cannot be guaranteed to be Euclidean butthe deleterious effects are trivial. Two examples of using themethod are shown. A set of 226 aligned globins were analysedand the resulting ordination very successfully represents theknown patterns of relationship between the sequences. In theother example, a set of 610 aligned 5S rRNA sequences were analysed.Sequence ordinations complement phylogenetic analyses. Theyshould not be viewed as a complete alternative.     

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.     

SNUFER: A software for localization and presentation of single nucleotide polymorphisms using a Clustal multiple sequence alignment output file

SNUFER is a software for the automatic localization and generation of tables used for the presentation of single nucleotide polymorphisms (SNPs). After input of a fasta file containing the sequences to be analyzed, a multiple sequence alignment is generated using ClustalW ran inside SNUFER. The ClustalW output file is then used to generate a table which displays the SNPs detected in the aligned sequences and their degree of similarity. This table can be exported to Microsoft Word, Microsoft Excel or as a single text file, permitting further editing for publication. The software was written using Delphi 7 for programming and FireBird 2.0 for sequence database management. It is freely available for noncommercial use and can be downloaded from http://www.heranza.com.br/bioinformatica2.htm.     

Multiple sequence alignment with hierarchical clustering.

An algorithm is presented for the multiple alignment of sequences, either proteins or nucleic acids, that is both accurate and easy to use on microcomputers. The approach is based on the conventional dynamic-programming method of pairwise alignment. Initially, a hierarchical clustering of the sequences is performed using the matrix of the pairwise alignment scores. The closest sequences are aligned creating groups of aligned sequences. Then close groups are aligned until all sequences are aligned in one group. The pairwise alignments included in the multiple alignment form a new matrix that is used to produce a hierarchical clustering. If it is different from the first one, iteration of the process can be performed. The method is illustrated by an example: a global alignment of 39 sequences of cytochrome c.     

Automated quantification of nuclear immunohistochemical markers with different complexity

Manual quantification of immunohistochemically stained nuclear markers is still laborious and subjective and the use of computerized systems for digital image analysis have not yet resolved the problems of nuclear clustering. In this study, we designed a new automatic procedure for quantifying various immunohistochemical nuclear markers with variable clustering complexity. This procedure consisted of two combined macros. The first, developed with a commercial software, enabled the analysis of the digital images using color and morphological segmentation including a masking process. All information extracted with this first macro was automatically exported to an Excel datasheet, where a second macro composed of four different algorithms analyzed all the information and calculated the definitive number of positive nuclei for each image. One hundred and eighteen images with different levels of clustering complexity was analyzed and compared with the manual quantification obtained by a trained observer. Statistical analysis indicated a great reliability (intra-class correlation coefficient > 0.950) and no significant differences between the two methods. Bland–Altman plot and Kaplan–Meier curves indicated that the results of both methods were concordant around 90% of analyzed images. In conclusion, this new automated procedure is an objective, faster and reproducible method that has an excellent level of accuracy, even with digital images with a high complexity.     

Was the ANITA rooting of the angiosperm phylogeny affected by long-branch attraction? Amborella, Nymphaeales, Illiciales, Trimeniaceae, and Austrobaileya

Five groups of basal angiosperms, Amborella, Nymphaeales, Illiciales, Trimeniaceae, and Austrobaileya (ANITA), were identified in several recent studies as representing a series of the earliest-diverging lineages of the angiosperm phylogeny. All of these studies except one employed a multigene analysis approach and used gymnosperms as the outgroup to determine the ingroup topology. The high level of divergence between gymnosperms and angiosperms, however, has long been implicated in the difficulty of reconstructing relationships at the base of angiosperm phylogeny using DNA sequences, for fear of long-branch attraction (LBA). In this study, we replaced the gymnosperm sequences from the five-gene matrix (mitochondrial atp1 and matR, plastid atpB and rbcL, and nuclear 18S rDNA) used in our earlier study with four categories of divergent sequences--random sequences with equal base frequencies or equally AT- and GC-rich contents, homopolymers and heteropolymers, misaligned gymnosperm sequences, and aligned lycopod and bryophyte sequences--to evaluate whether the gymnosperms were an appropriate outgroup to angiosperms in our earlier study that identified the ANITA rooting. All 24 analyses performed rooted the angiosperm phylogeny at either Acorus or Alisma (or Alisma-Triglochin-Potamogeton in one case due to use of a slightly different alignment) and placed the monocots as a basal grade, producing genuine LBA results. These analyses demonstrate that the identification of ANITA as the basalmost extant angiosperms was based on historical signals preserved in the gymnosperm sequences and that the gymnosperms were an appropriate outgroup with which to root the angiosperm phylogeny in the multigene sequence analysis. This strategy of evaluating the appropriateness of an outgroup using artificial sequences and a series of outgroups with increments of divergence levels can be applied to investigations of phylogenetic patterns at the bases of other major clades, such as land plants, animals, and eukaryotes.     

Overestimation of nonsynonymous/synonymous rate ratio by reverse-translation of aligned amino acid sequences

In the analysis of protein-coding nucleotide sequences, the ratio of the number of nonsynonymous substitutions to that of synonymous substitutions (d(N)/d(S)) is used as an indicator for the direction and magnitude of natural selection operating at the amino acid sequence level. The d(S) and d(N) values are estimated based on the comparison of homologous codons, which are often identified by converting (reverse-translating) aligned amino acid sequences into codon sequences. In this method, however, homologous codons may be mis-identified when frame-shifts occurred or amino acid sequences were mis-aligned, which may lead to overestimation of the d(N)/d(S) ratio. Here the effect of reverse-translating aligned amino acid sequences on the estimation of d(N)/d(S) ratio was examined through a large-scale analysis of protein-coding nucleotide sequences from vertebrate species. Apparently, 1-9% of codon sites that were identified as homologous with reverse-translation contained non-homologous codons, where the d(N)/d(S) ratio was unduly high. By correcting the d(N)/d(S) ratio for these codon sites, it was inferred that the ratio was 5-43% overestimated with reverse-translation. These results suggest that caution should be exerted in the study of natural selection using the d(N)/d(S) ratio by reverse-translating aligned amino acid sequences.     

Phylogenetics of notothenioid fishes (Teleostei: Acanthomorpha): inferences from mitochondrial and nuclear gene sequences

Notothenioids represent an adaptive radiation of teleost fishes in the frigid and ice-laden waters of the Southern Ocean surrounding Antarctica. Phylogenetic hypotheses for this clade have resulted primarily from analyses of mtDNA gene sequences, and studies utilizing nuclear gene DNA sequence data have focused on particular sub-clades of notothenioid fishes. In this study, we provide the first phylogenetic analysis of notothenioids using both mtDNA and nuclear gene sequences for a comprehensive sampling of all major lineages in the clade. Maximum parsimony and Bayesian analyses of aligned mtDNA genes, an aligned nuclear gene (S7 ribosomal protein intron 1), and combined dataset containing the mtDNA and nuclear genes resulted in phylogenies that contained the previously identified Antarctic and High Antarctic Clades. There were areas of agreement and disagreement between different datasets and methods of phylogenetic analysis, and the phylogenies resulting from the nuclear encoded S7 ribosomal protein intron 1 sequences were considerably less resolved than those inferred from mtDNA gene sequences. However, we anticipate increased resolution of the notothenioid phylogeny from future analyses that sample DNA sequences from several nuclear genes.     

Profile scanning for three-dimensional structural patterns in protein sequences

Profile analysis measures the similarity between a target sequenceand a group of aligned sequences (the probe). The probe sequencesare used to produce a position-specific scoring table (the profile)that can be aligned with any sequence (the target) using standarddynamic programming methods. We are developing a library ofprofiles, each describing a different structural motif. Thisallows any target sequence to be rapidly scanned for the presenceof structural motifs. Levels of significance for the comparisonof target sequences with the profile are determined in advance,permitting an objective decision to be made as to whether aprotein is likely to possess a structural motif. Received on July 17, 1987; accepted on January 4, 1988     

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 approximately 75% are longer than 899 bp. An RDP-II mirror site in Japan is now available (http://wdcm.nig.ac.jp/RDP/html/index.h tml). 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.