Recombination Analysis Tool (RAT): a program for the high-throughput detection of recombination

MOTIVATION: Recombination can be a prevailing drive in shaping genome evolution. RAT (Recombination Analysis Tool) is a Java-based tool for investigating recombination events in any number of aligned sequences (protein or DNA) of any length (short viral sequences to full genomes). It is an uncomplicated and intuitive application and allows the user to view only the regions of sequence alignments they are interested in. RESULTS: RAT was applied to viral sequences. Its utility was demonstrated through the detection of a known recombinant of HIV and a detailed analysis of Noroviruses, the most common cause of viral gastroenteritis in humans. AVAILABILITY: RAT, along with a user's guide, is freely available from http://jic-bioinfo.bbsrc.ac.uk/bioinformatics-research/staff/graham_etherington/RAT.htm.     

TreeMos: a high-throughput phylogenomic approach to find and visualize phylogenetic mosaicism

SUMMARY: TreeMos is a novel high-throughput graphical analysis application that allows the user to search for phylogenetic mosaicism among one or more DNA or protein sequence multiple alignments and additional unaligned sequences. TreeMos uses a sliding window and local alignment algorithm to identify the nearest neighbour of each sequence segment, and visualizes instances of sequence segments whose nearest neighbour is anomalous to that identified using the global alignment. Data sets can include whole genome sequences allowing phylogenomic analyses in which mosaicism may be attributed to recombination between any two points in the genome. TreeMos can be run from the command line, or within a web browser allowing the relationships between taxa to be explored by drill-through. AVAILABILITY: http://www2.warwick.ac.uk/fac/sci/whri/research/archaeobotany.     

Phylogenetic mapping of recombination hotspots in human immunodeficiency virus via spatially smoothed change-point processes

We present a Bayesian framework for inferring spatial preferences of recombination from multiple putative recombinant nucleotide sequences. Phylogenetic recombination detection has been an active area of research for the last 15 years. However, only recently attempts to summarize information from several instances of recombination have been made. We propose a hierarchical model that allows for simultaneous inference of recombination breakpoint locations and spatial variation in recombination frequency. The dual multiple change-point model for phylogenetic recombination detection resides at the lowest level of our hierarchy under the umbrella of a common prior on breakpoint locations. The hierarchical prior allows for information about spatial preferences of recombination to be shared among individual data sets. To overcome the sparseness of breakpoint data, dictated by the modest number of available recombinant sequences, we a priori impose a biologically relevant correlation structure on recombination location log odds via a Gaussian Markov random field hyperprior. To examine the capabilities of our model to recover spatial variation in recombination frequency, we simulate recombination from a predefined distribution of breakpoint locations. We then proceed with the analysis of 42 human immunodeficiency virus (HIV) intersubtype gag recombinants and identify a putative recombination hotspot.     

Chromosomes for molecular hybridization. Assignment of repetitive and single copy genes using a rapid filter-fixation method

Specific recombinant DNA sequences (5S rRNA, B1, albumin) were assigned to flow sorted chromosomes of the Chinese hamster cell line CHV79. For this purpose, a rapid protocol was developed using filterbound chromosomal DNA and probing with various nucleic acids, that allows sequence identification in chromosomes. A flow histogram and a flow karyogram of the CHV79 cell line were established by flow analysis in order to calculate the amount of DNA per CHV79 cell and their chromosomes. Subsequently, metaphase chromosomes or chromosomal groups were fractionated by electronic sorting and a defined number of chromosomes was directly bound to nitrocellulose filters for sequence homology analysis by a dot blot hybridization procedure. This procedure not only allows the assigning of specific DNA sequences to particular chromosomes, it is also applicable to studies of changes in karyotypes, for example translocations of given sequences.     

GARD: a genetic algorithm for recombination detection

MOTIVATION: Phylogenetic and evolutionary inference can be severely misled if recombination is not accounted for, hence screening for it should be an essential component of nearly every comparative study. The evolution of recombinant sequences can not be properly explained by a single phylogenetic tree, but several phylogenies may be used to correctly model the evolution of non-recombinant fragments. RESULTS: We developed a likelihood-based model selection procedure that uses a genetic algorithm to search multiple sequence alignments for evidence of recombination breakpoints and identify putative recombinant sequences. GARD is an extensible and intuitive method that can be run efficiently in parallel. Extensive simulation studies show that the method nearly always outperforms other available tools, both in terms of power and accuracy and that the use of GARD to screen sequences for recombination ensures good statistical properties for methods aimed at detecting positive selection. AVAILABILITY: Freely available http://www.datamonkey.org/GARD/     

EUGENE'HOM: A generic similarity-based gene finder using multiple homologous sequences

EUGENE'HOM is a gene prediction software for eukaryotic organisms based on comparative analysis. EUGENE'HOM is able to take into account multiple homologous sequences from more or less closely related organisms. It integrates the results of TBLASTX analysis, splice site and start codon prediction and a robust coding/non-coding probabilistic model which allows EUGENE'HOM to handle sequences from a variety of organisms. The current target of EUGENE'HOM is plant sequences. The EUGENE'HOM web site is available at http://genopole.toulouse.inra.fr/bioinfo/eugene/EuGeneHom/cgi-bin/EuGeneHom.pl.     

The Potential of Genomic Approaches to Rotifer Ecology

Rotifers are a key component of many freshwater ecosystems, but surveys of the composition of rotifer communities are limited by the labor-intensiveness of sample processing, particularly of non-planktonic taxa, and by the shortage of investigators qualified to identify a broad range of rotifer species. Additional problems are posed by species that must be identified from living specimens, and by members of cryptic species complexes. As DNA sequencing becomes easier and cheaper, it has become practical to obtain representative DNA sequences from identified rotifer species for use in genome-based surveys to determine which rotifers are present in a new sample, avoiding the difficulties of traditional surveys. Here we discuss two genome-based tools used in surveys of microbial communities: serial analysis of gene tags (SAGT) and microarray hybridization. SAGT is a method for inexpensively obtaining characteristic short DNA sequences from a sample that can both identify taxa for which the tag sequence is known and signal the presence of additional uncharacterized species. Microarray hybridization allows detection of DNA sequences in the sample that are identical or similar to sequences present on the microarray. We also report the construction and hybridization of a small microarray of rotifer sequences, demonstrating that this method can discriminate among bdelloid families, and is likely to make much finer discriminations if appropriate sequences are present on the microarray. These techniques are most powerful when combined with traditional systematics in collaborative efforts, which may be fostered through the data base of rotifer biology, WheelBase (http://jbpc.mbl.edu/wheelbase).     

PhyloPro: a web-based tool for the generation and visualization of phylogenetic profiles across Eukarya

SUMMARY: With increasing numbers of eukaryotic genome sequences, phylogenetic profiles of eukaryotic genes are becoming increasingly informative. Here, we introduce a new web-tool Phylopro (http://compsysbio.org/phylopro/), which uses the 120 available eukaryotic genome sequences to visualize the evolutionary trajectories of user-defined subsets of model organism genes. Applied to pathways or complexes, PhyloPro allows the user to rapidly identify core conserved elements of biological processes together with those that may represent lineage-specific innovations. PhyloPro thus provides a valuable resource for the evolutionary and comparative studies of biological systems.     

HUGE: a database for human large proteins identified in the Kazusa cDNA sequencing project

We have been developing a HUGE database to summarize results from the sequence analysis of human novel large (>4 kb) cDNAs identified in the Kazusa cDNA sequencing project, systematically designated KIAA plus a four-digit number. HUGE currently contains nearly 2000 gene/protein characteristic tables harboring the results of the computer-assisted analysis of the cDNA and the predicted protein sequences together with those of expression profiling and chromosomal mapping. In the updated version of HUGE, we made it possible to compare each KIAA cDNA sequence with the corresponding entry in the human draft genome sequence that was published recently. Approximately 90% of KIAA cDNAs in HUGE can be localized along the human genome for at least half or more of the cDNA’s length. Any nucleotide differences between the cDNA and the corresponding genomic sequences are also presented in detail. This new version of HUGE greatly helps us evaluate the completeness of cDNA clones and the accuracy of cDNA/genomic sequences. More interestingly, in some cases, the ability to compare cDNA with genomic sequences allows us to identify candidate sites of RNA editing. HUGE is available on the World Wide Web at http://www.kazusa.or.jp/huge.     

Chromosomes for molecular hybridization

Summary Specific recombinant DNA sequences (5S rRNA, B1, albumin) were assigned to flow sorted chromosomes of the Chinese hamster cell line CHV79. For this purpose, a rapid protocol was developed using filterbound chromosomal DNA and probing with various nucleic acids, that allows sequence identification in chromosomes. A flow histogram and a flow karyogram of the CHV79 cell line were established by flow analysis in order to calculate the amount of DNA per CHV79 cell and their chromosomes. Subsequently, metaphase chromosomes or chromosomal groups were fractionated by electronic sorting and a defined number of chromosomes was directly bound to nitrocellulose filters for sequence homology analysis by a dot blot hybridization procedure. This procedure not only allows the assigning of specific DNA sequences to particular chromosomes, it is also applicable to studies of changes in karyotypes, for example translocations of given sequences.Some results shown constituted a Diploma Thesis by G.L. submitted to and accepted by the Department of Biology, University of Kaiserslautern     

Construction, analysis, and application to 46,XY gonadal dysgenesis of a recombinant phage DNA library from flow-sorted human Y chromosomes

The analysis of a recombinant human Y-enriched Hind III total digest phage library prepared from the DNA of flow sorted human Y chromosomes is described. Out of 43 phage inserts from the library thus far mapped, 25 revealed hybridization with Y chromosomal DNA. These inserts may be divided into five groups according to their degree of Y specific hybridization: inserts that hybridize with one single copy or slightly repeated Y-specific DNA sequence, Y-specific repeated sequences of various restriction fragment lengths, Y-chromosomal DNA sequence(s) shared by a sequence on the X and/or on autosomes, Y-specific DNA sequences in addition to multiple X and/or autosomal sequences, or Y-specific repeated DNA in addition to multiple X and/or autosomal sequences. Application of probes from this library for diagnostic purposes is shown in two 46,XY patients with gonadal dysgenesis and small deletions of the Y short arm.     

JaDis: computing distances between nucleic acid sequences.

SUMMARY: JaDis is a Java application for computing evolutionary distances between nucleic acid sequences and G+C base frequencies. It allows specific comparison of coding sequences, of non-coding sequences or of a non-coding sequence with coding sequences. AVAILABILITY: http://pbil.univ-lyon1.fr/software/jadis.html     

An exact nonparametric method for inferring mosaic structure in sequence triplets

Statistical tests for detecting mosaic structure or recombination among nucleotide sequences usually rely on identifying a pattern or a signal that would be unlikely to appear under clonal reproduction. Dozens of such tests have been described, but many are hampered by long running times, confounding of selection and recombination, and/or inability to isolate the mosaic-producing event. We introduce a test that is exact, nonparametric, rapidly computable, free of the infinite-sites assumption, able to distinguish between recombination and variation in mutation/fixation rates, and able to identify the breakpoints and sequences involved in the mosaic-producing event. Our test considers three sequences at a time: two parent sequences that may have recombined, with one or two breakpoints, to form the third sequence (the child sequence). Excess similarity of the child sequence to a candidate recombinant of the parents is a sign of recombination; we take the maximum value of this excess similarity as our test statistic Delta(m,n,b). We present a method for rapidly calculating the distribution of Delta(m,n,b) and demonstrate that it has comparable power to and a much improved running time over previous methods, especially in detecting recombination in large data sets.     

ecoPrimers: inference of new DNA barcode markers from whole genome sequence analysis

Using non-conventional markers, DNA metabarcoding allows biodiversity assessment from complex substrates. In this article, we present ecoPrimers, a software for identifying new barcode markers and their associated PCR primers. ecoPrimers scans whole genomes to find such markers without a priori knowledge. ecoPrimers optimizes two quality indices measuring taxonomical range and discrimination to select the most efficient markers from a set of reference sequences, according to specific experimental constraints such as marker length or specifically targeted taxa. The key step of the algorithm is the identification of conserved regions among reference sequences for anchoring primers. We propose an efficient algorithm based on data mining, that allows the analysis of huge sets of sequences. We evaluate the efficiency of ecoPrimers by running it on three different sequence sets: mitochondrial, chloroplast and bacterial genomes. Identified barcode markers correspond either to barcode regions already in use for plants or animals, or to new potential barcodes. Results from empirical experiments carried out on a promising new barcode for analyzing vertebrate diversity fully agree with expectations based on bioinformatics analysis. These tests demonstrate the efficiency of ecoPrimers for inferring new barcodes fitting with diverse experimental contexts. ecoPrimers is available as an open source project at: http://www.grenoble.prabi.fr/trac/ecoPrimers.     

The use of R-looping for structural gene identification and mRNA purification.

A method is presented for the purification of mRNAs and the identification of structural gene sequences in recombinant DNA molecules. RNA is hybridized to double-stranded linear DNA such that R-loops are formed between most DNAs and their complementary RNA sequences. These R-loops are purified from unhybridized RNAs by gel filtration chromatography in the presence of a high concentration of salt. The complementary RNAs are released from the R-loops by heating, and are assayed by gel electrophoresis or cell free translation to determine their purity and to identify the proteins for which they code. We have demonstrated that recombinant DNAs containing sequences for abundant or moderately abundant mRNAs of Saccharomyces cerevisiae can be identified by this means.     

SciRoKo: a new tool for whole genome microsatellite search and investigation

SciRoKo is a user-friendly software tool for the identification of microsatellites in genomic sequences. The combination of an extremely fast search algorithm with a built-in summary statistic tool makes SciRoKo an excellent tool for full genome analysis. Compared to other already existing tools, SciRoKo also allows the analysis of compound microsatellites. AVAILABILITY: free for use: www.kofler.or.at/Bioinformatics. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.     

OligoMatcher

OligoMatcher is a web-based tool for analysis and selection of unique oligonucleotide sequences for gene silencing by antisense oligonucleotides (ASOs) or small interfering RNA (siRNA). A specific BLAST server was built for analysing sequences of ASOs that target pre-mRNA in the cell nucleus. Tissue- and cell-specific expression data of potential cross-reactive genes are integrated in the OligoMatcher program, which allows biologists to select unique oligonucleotide sequences for their target genes in specific experimental systems. AVAILABILITY: The OligoMatcher web server is available at http://shelob.cs.iupui.edu:18081/oligomatch.php. The source code is freely available for non-profit use on request to the authors. CONTACT: Mathew Palakal (mpalakal@cs.iupui.edu) or Shuyu Li (li_shuyu_dan@lilly.com).     

Signal-CF: a subsite-coupled and window-fusing approach for predicting signal peptides

We have developed an automated method for predicting signal peptide sequences and their cleavage sites in eukaryotic and bacterial protein sequences. It is a 2-layer predictor: the 1st-layer prediction engine is to identify a query protein as secretory or non-secretory; if it is secretory, the process will be automatically continued with the 2nd-layer prediction engine to further identify the cleavage site of its signal peptide. The new predictor is called Signal-CF, where C stands for "coupling" and F for "fusion", meaning that Signal-CF is formed by incorporating the subsite coupling effects along a protein sequence and by fusing the results derived from many width-different scaled windows through a voting system. Signal-CF is featured by high success prediction rates with short computational time, and hence is particularly useful for the analysis of large-scale datasets. Signal-CF is freely available as a web-server at http://chou.med.harvard.edu/bioinf/Signal-CF/ or http://202.120.37.186/bioinf/Signal-CF/.     

ACT: the Artemis Comparison Tool

The Artemis Comparison Tool (ACT) allows an interactive visualisation of comparisons between complete genome sequences and associated annotations. The comparison data can be generated with several different programs; BLASTN, TBLASTX or Mummer comparisons between genomic DNA sequences, or orthologue tables generated by reciprocal FASTA comparison between protein sets. It is possible to identify regions of similarity, insertions and rearrangements at any level from the whole genome to base-pair differences. ACT uses Artemis components to display the sequences and so inherits powerful searching and analysis tools. ACT is part of the Artemis distribution and is similarly open source, written in Java and can run on any Java enabled platform, including UNIX, Macintosh and Windows.