COMPAM :visualization of combining pairwise alignments for multiple genomes

COMPAM is a tool for visualizing relationships among multiple whole genomes by combining all pairwise genome alignments. It displays shared conserved regions (blocks) and where these blocks occur (edges) as block relation graphs which can be explored interactively. An unannotated genome, e.g. can then be explored using information from well-annotated genomes, COG-based genome annotation and genes. COMPAM can run either as a stand-alone application or through an applet that is provided as service to PLATCOM, a toolset for whole genome comparative analysis, where a wide variety of genomes can be easily selected. Features provided by COMPAM include the ability to export genome relationship information into file formats that can be used by other existing tools. AVAILABILITY: http://bio.informatics.indiana.edu/projects/compam/     

OWEN: aligning long collinear regions of genomes

OWEN is an interactive tool for aligning two long DNA sequences that represents similarity between them by a chain of collinear local similarities. OWEN employs several methods for constructing and editing local similarities and for resolving conflicts between them. Alignments of sequences of lengths over 10(6) can often be produced in minutes. OWEN requires memory below 20 L, where L is the sum of lengths of the compared sequences.     

Analysis of collinear regions of Oryza AA and CC genomes

Comparative analyses of genome structure and sequence of closely related species have yielded insights into the evolution and function of plant genomes. A total of 103,844 BAC end sequences delegated -73.8 Mb of O. officinalis that belongs to the CC genome type of the rice genus Oryza were obtained and compared with the genome sequences office cultivar, O. sativa ssp.japonica cv. Nipponbare. We found that more than 45% of O. officinalis genome consists of repeat sequences, which is higher than that of Nipponbare cultivar. To further investigate the evolutionary divergence of AA and CC genomes, two BAC-contigs of O. officinalis were compared with the collinear genomic regions of Nipponbare. Of 57 genes predicted in the AA genome orthologous regions, 39 had orthologs in the regions of the CC genome. Alignment of the orthologous regions indicated that the CC genome has undergone expansion in both genic and intergenic regions through primarily retroelement insertion. Particularly, the density of RNA transposable elements was 17.95% and 1.78% in O. officinalis and O. sativa, respectively. This explains why the orthologous region is about 100 kb longer in the CC genome in comparison to the AA genome.     

A hierarchical approach to aligning collinear regions of genomes

MOTIVATION: As a first approximation, similarity between two long orthologous regions of genomes can be represented by a chain of local similarities. Within such a chain, pairs of successive similarities are collinear (non-conflicting), i.e. segments involved in the nth similarity precede in both sequences segments involved in the (n+1)th similarity. However, when all similarities between two long sequences are considered, usually there are many conflicts between them. Although some conflicts can be avoided by masking transposons or low-complexity sequences, selecting only those similarities that reflect orthology and, thus, belong to the evolutionarily true chain is not trivial. RESULTS: We propose a simple, hierarchical algorithm of finding the true chain of local similarities. Starting from similarities with low P-values, we resolve each pairwise conflict by deleting a similarity with a higher P-value. This greedy approach constructs a chain of similarities faster than when a chain optimal with respect to some global criterion is sought, and makes more sense biologically.     

Comparative analysis of GS and BK virus genomes.

Analysis of heteroduplexing between the genomes of GS virus, a BK-group virus, and the prototype BK virus revealed one region of nonhomology. Further analysis by cleavage of viral DNA with the restriction endonucleases EcoRI, HindIII, and HaeIII revealed that base changes in the GS virus genome spanned 0.6 to 0.7 map unit from the EcoRI site. Large T and small t antigens of GS virus appear to be similar in size to the BK virus antigens.     

Integrated functional visualization of eukaryotic genomes

Background  

Increasing amounts of data from large scale whole genome analysis efforts demands convenient tools for manipulation, visualization and investigation. Whole genome plots offer an intuitive window to the analysis. We describe two applications that enable users to easily plot and explore whole genome data from their own or other researchers' experiments.     

YOC,A new strategy for pairwise alignment of collinear genomes

Background

Comparing and aligning genomes is a key step in analyzing closely related genomes. Despite the development of many genome aligners in the last 15 years, the problem is not yet fully resolved, even when aligning closely related bacterial genomes of the same species. In addition, no procedures are available to assess the quality of genome alignments or to compare genome aligners.

Results

We designed an original method for pairwise genome alignment, named YOC, which employs a highly sensitive similarity detection method together with a recent collinear chaining strategy that allows overlaps. YOC improves the reliability of collinear genome alignments, while preserving or even improving sensitivity. We also propose an original qualitative evaluation criterion for measuring the relevance of genome alignments. We used this criterion to compare and benchmark YOC with five recent genome aligners on large bacterial genome datasets, and showed it is suitable for identifying the specificities and the potential flaws of their underlying strategies.

Conclusions

The YOC prototype is available at https://github.com/ruricaru/YOC. It has several advantages over existing genome aligners: (1) it is based on a simplified two phase alignment strategy, (2) it is easy to parameterize, (3) it produces reliable genome alignments, which are easier to analyze and to use.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0530-3) contains supplementary material, which is available to authorized users.     

Dynamic gene copy number variation in collinear regions of grass genomes

A salient feature of genomes of higher organisms is the birth and death of gene copies. An example is the alpha prolamin genes, which encode seed storage proteins in grasses (Poaceae) and represent a medium-size gene family. To better understand the mechanism, extent, and pace of gene amplification, we compared prolamin gene copies in the genomes of two different tribes in the Panicoideae, the Paniceae and the Andropogoneae. We identified alpha prolamin (setarin) gene copies in the diploid foxtail millet (Paniceae) genome (490 Mb) and compared them with orthologous regions in diploid sorghum (730 Mb) and ancient allotetraploid maize (2,300 Mb) (Andropogoneae). Because sequenced genomes of other subfamilies of Poaceae like rice (389 Mb) (Ehrhartoideae) and Brachypodium (272 Mb) (Pooideae) do not have alpha prolamin genes, their collinear regions can serve as "empty" reference sites. A pattern emerged, where genes were copied and inserted into other chromosomal locations followed by additional tandem duplications (clusters). We observed both recent (species-specific) insertion events and older ones that are shared by these tribes. Many older copies were deleted by unequal crossing over of flanking sequences or damaged by truncations. However, some remain intact with active and inactive alleles. These results indicate that genomes reflect only a snapshot of the gene content of a species and are far less static than conventional genetics has suggested. Nucleotide substitution rates for active alpha prolamins genes were twice as high as for low copy number beta, gamma, and delta prolamin genes, suggesting that gene amplification accelerates the pace of divergence.     

MuGeN: simultaneous exploration of multiple genomes and computer analysis results

MOTIVATION: The availability of increasing amounts of sequence data about completely sequenced genomes spurs the development of new methods in the fields of automated annotation, and of comparative genomics. Tools allowing the visualization of results produced by analysis methods, superimposed on possibly annotated sequence data, and enabling synchronized navigation in multiple genomes, provide new means for interactive genome exploration. This kind of visual inspection can be used as a basis to assess the quality of new analysis algorithms, or to discover genome portions to be subjected to in-depth studies. RESULTS: We propose a software package, MuGeN, built for navigating through multiple annotated genomes. It is capable of retrieving annotated sequences in several formats, stored in local files, or available in databases over the network. From these, it then generates an interactive display, or an image file, in most common formats suitable for printing, further editing or integrating in Web pages. Genome maps may be mixed with computer analysis results loaded from XML files, whose format is generic enough to be adapted to a majority of sequence oriented analysis methods. AVAILABILITY: MuGeN is available at http://www-mig.jouy.inra.fr/bdsi/MuGeN.     

Genome analysis: A new approach for visualization of sequence organization in genomes

In this article we describe and demonstrate the versatility of a computer program, GENOME MAPPING, that uses interactive graphics and runs on an IRIS workstation. The program helps to visualize as well as analyse global and local patterns of genomic DNA sequences. It was developed keeping in mind the requirements of the human genome sequencing programme, which requires rapid analysis of the data. Using GENOME MAPPING one can discern signature patterns of different kinds of sequences and analyse such patterns for repetitive as well as rare sequence strings. Further, one can visualize the extent of global homology between different genomic sequences. An application of our method to the published yeast mitochondrial genome data shows similar sequence organizations in the entire sequence and in smaller subsequences     

Comparative analysis of processed pseudogenes in the mouse and human genomes

Pseudogenes are important resources in evolutionary and comparative genomics because they provide molecular records of the ancient genes that existed in the genome millions of years ago. We have systematically identified approximately 5000 processed pseudogenes in the mouse genome, and estimated that approximately 60% are lineage specific, created after the mouse and human diverged. In both mouse and human genomes, similar types of genes give rise to many processed pseudogenes. These tend to be housekeeping genes, which are highly expressed in the germ line. Ribosomal-protein genes, in particular, form the largest sub-group. The processed pseudogenes in the mouse occur with a distinctly different chromosomal distribution than LINEs or SINEs - preferentially in GC-poor regions. Finally, the age distribution of mouse-processed pseudogenes closely resembles that of LINEs, in contrast to human, where the age distribution closely follows Alus (SINEs).     

Comparative analysis of RNA genomes of mouse hepatitis viruses.

The RNA genomes of various murine hepatitis virus (MHV) strains were studied by T1-oligonucleotide fingerprinting analysis with regard to their structure and sequence relationship. It was found that the MHV particles contained only positive-stranded 60S RNA which had a "cap" structure at its 5' end. No negative-stranded RNA was found. It was also shown that most of the MHV strains had diverged quite extensively in their genetic sequences. However, MHV-3, a hepatotropic strain, and A59, a nonpathogenic strain, were found to have very similar oligonucleotide fingerprinting patterns. Yet, each of them contained two to four specific oligonucleotides. The MHV-3-specific oligonucleotides were conserved in almost all of the hepatotropic MHV strains studied. In contrast, two of the A59-specific oligonucleotides were absent from the genomes of all hepatotropic strains. These findings suggest that these unique oligonucleotides might be localized at the genetic region(s) associated with viral pathogenicity or other biological properties of the virus. Comparison of viral structural proteins also suggests that MHV-3 and A59 are more closely related than other MHV strains. The significance of these findings is discussed.     

Comparative analysis of the genomes of feline leukemia viruses.

The genomes of several strains of feline leukemia virus (FeLV) were compared by two-dimensional polyacrylamide gel electrophoresis of the large RNase T1-resistant oligonucleotides of the 70S RNA. Differences between each strain of FeLV tested were detected by this method. We estimate that the degree of sequence identity between the viruses is: FeLV A (Glasgow-1) to FeLV B (Snyder-Theilen), 52%; FeLV A (Glasgow-1) to FeLV C(Sarma), 66%; FeLV B(Snyder-Theilen) to FeLV C (Sarma), 37%. The fingerprints of two independent isolates of FeLV strains of subgroup A (Glasgow-1 and Rickard) were detectably different. We conclude that the RNase T1 oligonucleotide fingerprint pattern provides a useful tool for identification of FeLV strains.     

Comparative sequence and methylation analysis of chloroplast and amyloplast genomes from rice

Plant Molecular Biology - Grain amyloplast and leaf chloroplast DNA sequences are identical in rice plants but are differentially methylated. The leaf chloroplast DNA becomes more methylated as the...     

Comparative analysis of Acinetobacters: three genomes for three lifestyles

Acinetobacter baumannii is the source of numerous nosocomial infections in humans and therefore deserves close attention as multidrug or even pandrug resistant strains are increasingly being identified worldwide. Here we report the comparison of two newly sequenced genomes of A. baumannii. The human isolate A. baumannii AYE is multidrug resistant whereas strain SDF, which was isolated from body lice, is antibiotic susceptible. As reference for comparison in this analysis, the genome of the soil-living bacterium A. baylyi strain ADP1 was used. The most interesting dissimilarities we observed were that i) whereas strain AYE and A. baylyi genomes harbored very few Insertion Sequence elements which could promote expression of downstream genes, strain SDF sequence contains several hundred of them that have played a crucial role in its genome reduction (gene disruptions and simple DNA loss); ii) strain SDF has low catabolic capacities compared to strain AYE. Interestingly, the latter has even higher catabolic capacities than A. baylyi which has already been reported as a very nutritionally versatile organism. This metabolic performance could explain the persistence of A. baumannii nosocomial strains in environments where nutrients are scarce; iii) several processes known to play a key role during host infection (biofilm formation, iron uptake, quorum sensing, virulence factors) were either different or absent, the best example of which is iron uptake. Indeed, strain AYE and A. baylyi use siderophore-based systems to scavenge iron from the environment whereas strain SDF uses an alternate system similar to the Haem Acquisition System (HAS). Taken together, all these observations suggest that the genome contents of the 3 Acinetobacters compared are partly shaped by life in distinct ecological niches: human (and more largely hospital environment), louse, soil.     

Comparative analysis of variable regions in the genomes of variola virus

Nucleotide sequences of two extended segments of the terminal variable regions in variola virus genome were determined. The size of the left segment was 13.5 kbp and of the right, 10.5 kbp. Totally, over 540 kbp were sequenced for 22 variola virus strains. The conducted phylogenetic analysis and the data published earlier allowed us to find the interrelations between 70 variola virus isolates, the character of their clustering, and the degree of intergroup and intragroup variations of the clusters of variola virus strains. The most polymorphic loci of the genome segments studied were determined. It was demonstrated that that these loci are localized to either noncoding genome regions or to the regions of destroyed open reading frames, characteristic of the ancestor virus. These loci are promising for development of the strategy for genotyping variola virus strains. Analysis of recombination using various methods demonstrated that, with the only exception, no statistically significant recombinational events in the genomes of variola virus strains studied were detectable.