Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Genome alignment”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,531 records · Page 85Linked to original sources

Replication-associated strand asymmetries in mammalian genomes: toward detection of replication origins.

In the course of evolution, mutations do not affect both strands of genomic DNA equally. This imbalance mainly results from asymmetric DNA mutation and repair processes associated with replication and transcription. In prokaryotes, prevalence of G over C and T over A is frequently observed in the leading strand. The sign of the resulting TA and GC skews changes abruptly when crossing replication-origin and termination sites, producing characteristic step-like transitions. In mammals, transcription-coupled skews have been detected, but so far, no bias has been associated with replication. Here, analysis of intergenic and transcribed regions flanking experimentally identified human replication origins and the corresponding mouse and dog homologous regions demonstrates the existence of compositional strand asymmetries associated with replication. Multiscale analysis of human genome skew profiles reveals numerous transitions that allow us to identify a set of 1,000 putative replication initiation zones. Around these putative origins, the skew profile displays a characteristic jagged pattern also observed in mouse and dog genomes. We therefore propose that in mammalian cells, replication termination sites are randomly distributed between adjacent origins. Taken together, these analyses constitute a step toward genome-wide studies of replication mechanisms.

Animals↗

Archaebacterial DNA-dependent RNA polymerases testify to the evolution of the eukaryotic nuclear genome.

Genes for DNA-dependent RNA polymerase components B, A, and C from the archaebacterium Sulfolobus acidocaldarius and for components B", B', A, and C from the archaebacterium Halobacterium halobium were cloned and sequenced. They are organized in gene clusters in the order above, which corresponds to the order of the homologous rpoB and rpoC genes in the corresponding operon of the Escherichia coli genome. Derived amino acid sequences of archaebacterial components A and C were aligned with each other and with the sequences of corresponding (largest) subunits from the archaebacterium Methanobacterium thermoautotrophicum, with sequences of various eukaryotic nuclear RNA polymerases I, II, and III, and with the sequence of the beta' component from E. coli polymerase. The archaebacterial genes for component A are homologous to about the first two-thirds of genes for the eukaryotic component A and the eubacterial component beta', and the archaebacterial genes for component C are homologous to the last third of the genes for the eukaryotic component A and the eubacterial component beta'. Unrooted phylogenetic dendrograms derived from both distance matrix and parsimony analyses show the archaebacteria are a coherent group closely related to the eukaryotic nuclear RNA polymerase II and/or III lineages. The eukaryotic polymerase I lineage appears to arise separately from a bifurcation with the eubacterial beta' component lineage.

Amino Acid Sequence↗

Conservation of synteny between the genome of the pufferfish (Fugu rubripes) and the region on human chromosome 14 (14q24.3) associated with familial Alzheimer disease (AD3 locus)

The genome of the pufferfish (Fugu rubripes) (400 Mb) is approximately 7.5 times smaller than the human genome, but it has a similar gene repertoire to that of man. If regions of the two genomes exhibited conservation of gene order (i.e., were syntenic), it should be possible to reduce dramatically the effort required for identification of candidate genes in human disease loci by sequencing syntenic regions of the compact Fugu genome. We have demonstrated that three genes (dihydrolipoamide succinyltransferase, S31iii125, and S20i15), which are linked to FOS in the familial Alzheimer disease focus (AD3) on human chromosome 14, have homologues in the Fugu genome adjacent to Fugu cFOS. The relative gene order of cFOS, S31iii125, and S20i15 was the same in both genomes, but in Fugu these three genes lay within a 12.4-kb region, compared to >600 kb in the human AD3 locus. These results demonstrate the conservation of synteny between the genomes of Fugu and man and highlight the utility of this approach for sequence-based identification of genes in human disease loci.

Acyltransferases↗

Transposon diversity in Arabidopsis thaliana.

Recent availability of extensive genome sequence information offers new opportunities to analyze genome organization, including transposon diversity and accumulation, at a level of resolution that was previously unattainable. In this report, we used sequence similarity search and analysis protocols to perform a fine-scale analysis of a large sample ( approximately 17.2 Mb) of the Arabidopsis thaliana (Columbia) genome for transposons. Consistent with previous studies, we report that the A. thaliana genome harbors diverse representatives of most known superfamilies of transposons. However, our survey reveals a higher density of transposons of which over one-fourth could be classified into a single novel transposon family designated as Basho, which appears unrelated to any previously known superfamily. We have also identified putative transposase-coding ORFs for miniature inverted-repeat transposable elements (MITEs), providing clues into the mechanism of mobility and origins of the most abundant transposons associated with plant genes. In addition, we provide evidence that most mined transposons have a clear distribution preference for A + T-rich sequences and show that structural variation for many mined transposons is partly due to interelement recombination. Taken together, these findings further underscore the complexity of transposons within the compact genome of A. thaliana.

Amino Acid Sequence↗

DNA binding and bending by the human papillomavirus type 16 E2 protein. Recognition of an extended binding site.

The human papillomavirus (HPV) 16 E2 protein (hE2) binds to four sites present upstream of the P97 promoter and regulates transcription of the viral E6 and E7 oncogenes. We have determined the relative binding constants for the interaction of the full-length hE2 protein with these sites. Our results show that hE2 binds tightly to site 4, less tightly to sites 1 and 2, and weakly to site 3. Similar results have previously been obtained using a C-terminal fragment of the hE2 protein suggesting that the C-terminal domain is the sole determinant of DNA binding affinity and specificity. Using circular permutation assays we show that binding of the hE2 protein induces the formation of a significant DNA bend and that the hE2-induced DNA bend angle is the same at both tight and weak hE2-binding sites. An alignment of the four hE2-binding sites from the HPV 16 genome suggests that this protein recognizes an extended binding site when compared with the bovine papillomavirus E2 protein. Here we show that the hE2 protein binds tightly to sites containing an A:T or a G:C base pair at position 7 of its binding site but weakly to sites containing either C:G or T:A at this position. Using site-directed mutagenesis we demonstrate that an arginine at position 304 of the hE2 protein is responsible for the recognition of specific base pairs at this position.

Binding Sites↗

From Jabberwocky to genome: Lewis Carroll and computational biology.

In addition to his literary output, Lewis Carroll created a vast range of games and puzzles that depend upon wordplay of various kinds, especially the manipulation of alphabetic symbols in diverse contexts. Such wordplay reveals a turn of mind well suited to methodologies used in modern computational biology.

Computational Biology↗

Identifying conserved gene clusters in the presence of homology families.

The study of conserved gene clusters is important for understanding the forces behind genome organization and evolution, as well as the function of individual genes or gene groups. In this paper, we present a new model and algorithm for identifying conserved gene clusters from pairwise genome comparison. This generalizes a recent model called "gene teams." A gene team is a set of genes that appear homologously in two or more species, possibly in a different order yet with the distance of adjacent genes in the team for each chromosome always no more than a certain threshold. We remove the constraint in the original model that each gene must have a unique occurrence in each chromosome and thus allow the analysis on complex prokaryotic or eukaryotic genomes with extensive paralogs. Our algorithm analyzes a pair of chromosomes in O(mn) time and uses O(m+n) space, where m and n are the number of genes in the respective chromosomes. We demonstrate the utility of our methods by studying two bacterial genomes, E. coli K-12 and B. subtilis. Many of the teams identified by our algorithm correlate with documented E. coli operons, while several others match predicted operons, previously suggested by computational techniques. Our implementation and data are publicly available at euler.slu.edu/ approximately goldwasser/homologyteams/.

Bacillus subtilis↗

DNAskew: statistical analysis of base compositional asymmetry and prediction of replication boundaries in the genome sequences.

Sueoka and Lobry declared respectively that, in the absence of bias between the two DNA strands for mutation and selection, the base composition within each strand should be A=T and C=G (this state is called Parity Rule type 2, PR2). However, the genome sequences of many bacteria, vertebrates and viruses showed asymmetries in base composition and gene direction. To determine the relationship of base composition skews with replication orientation, gene function, codon usage biases and phylogenetic evolution, in this paper a program called DNAskew was developed for the statistical analysis of strand asymmetry and codon composition bias in the DNA sequence. In addition, the program can also be used to predict the replication boundaries of genome sequences. The method builds on the fact that there are compositional asymmetries between the leading and the lagging strand for replication. DNAskew was written in Perl script language and implemented on the LINUX operating system. It works quickly with annotated or unannotated sequences in GBFF (GenBank flatfile) or fasta format. The source code is freely available for academic use at http://www.epizooty.com/pub/stat/DNAskew.

Algorithms↗

SAMBA: hardware accelerator for biological sequence comparison.

MOTIVATION: SAMBA (Systolic Accelerator for Molecular Biological Applications) is a 128 processor hardware accelerator for speeding up the sequence comparison process. The short-term objective is to provide a low-cost board to boost PC or workstation performance on this class of applications. This paper places SAMBA amongst other existing systems and highlights the original features. RESULTS: Real performance obtained from the prototype is demonstrated. For example, a sequence of 300 amino acids is scanned against SWISS-PROT-34 (21 210 389 residues) in 30 s using the Smith and Waterman algorithm. More time-consuming applications, like the bank-to-bank comparison, are computed in a few hours instead of days on standard workstations. Technology allows the prototype to fit onto a single PCI board for plugging into any PC or workstation. AVAILABILITY: SAMBA can be tested on the WEB server at URL http://www.irisa.fr/SAMBA/.

Algorithms↗

bioWidgets: data interaction components for genomics.

MOTIVATION: The presentation of genomics data in a perspicuous visual format is critical for its rapid interpretation and validation. Relatively few public database developers have the resources to implement sophisticated front-end user interfaces themselves. Accordingly, these developers would benefit from a reusable toolkit of user interface and data visualization components. RESULTS: We have designed the bioWidget toolkit as a set of JavaBean components. It includes a wide array of user interface components and defines an architecture for assembling applications. The toolkit is founded on established software engineering design patterns and principles, including componentry, Model-View-Controller, factored models and schema neutrality. As a proof of concept, we have used the bioWidget toolkit to create three extendible applications: AnnotView, BlastView and AlignView.

Amino Acid Sequence↗

Pro-Frame: similarity-based gene recognition in eukaryotic DNA sequences with errors.

Performance of existing algorithms for similarity-based gene recognition in eukaryotes drops when the genomic DNA has been sequenced with errors. A modification of the spliced alignment algorithm allows for gene recognition in sequences with errors, in particular frameshifts. It tolerates up to 5% of sequencing errors without considerable drop of prediction reliability when a sufficiently close homologous protein is available (normalized evolutionary distance similarity score 50% or higher).

Algorithms↗

PRIMEGENS: robust and efficient design of gene-specific probes for microarray analysis.

MOTIVATION: DNA microarray is a powerful high-throughput tool for studying gene function and regulatory networks. Due to the problem of potential cross hybridization, using full-length genes for microarray construction is not appropriate in some situations. A bioinformatic tool, PRIMEGENS, has recently been developed for the automatic design of PCR primers using DNA fragments that are specific to individual open reading frames (ORFs). RESULTS: PRIMEGENS first carries out a BLAST search for each target ORF against all other ORFs of the genome to quickly identify possible homologous sequences. Then it performs optimal sequence alignment between the target ORF and each of its homologous ORFs using dynamic programming. PRIMEGENS uses the sequence alignments to select gene- specific fragments, and then feeds the fragments to the Primer3 program to design primer pairs for PCR amplification. PRIMEGENS can be run from the command line on Unix/Linux platforms as a stand-alone package or it can be used from a Web interface. The program runs efficiently, and it takes a few seconds per sequence on a typical workstation. PCR primers specific to individual ORFs from Shewanella oneidensis MR-1 and Deinococcus radiodurans R1 have been designed. The PCR amplification results indicate that this method is very efficient and reliable for designing specific probes for microarray analysis.

Algorithms↗

Detecting genomic features under weak selective pressure: the example of codon usage in animals and plants.

Large scale experiments of gene inactivation in yeast have shown that 50% of genes have no detectable impact on the phenotype, and similar observations have been made in other model organisms. This apparent paradox is probably due to the fact that many genes only have a marginal contribution to the fitness of organisms. Because of the size of populations and the number of generations that can be studied in laboratories, experimental approaches only permit to detect functional elements that have a strong phenotypic impact. Comparative sequence analysis can help to solve this problem: the analysis of sequences evolution permits to detect the action of selection, and hence to reveal functional features of genomes. This approach will be illustrated by the study of synonymous codon usage in animals and plants.

Animals↗

The fuzzy polynucleotide space: basic properties.

MOTIVATION: Any triplet codon may be regarded as a 12-dimensional fuzzy code. Sufficient information about a particular sequence may not be available in certain situations. The investigator will be confronted with imprecise sequences, yet want to make comparisons of sequences. Fuzzy polynucleotides can be compared by using geometrical interpretation of fuzzy sets as points in a hypercube. RESULTS: We introduce the space of fuzzy polynucleotides and a means of measuring dissimilitudes between them. We establish mathematical principles to measure dissimilarities between fuzzy polynucleotides and present several examples in this metric space. We calculate the frequencies of the nucleotides at the three base sites of a codon in the coding sequences of Escherichia coli K-12 and Mycobacterium tuberculosis H37Rv, and consider them as points in that fuzzy space. We compute the distance between the genomes of E.coli and M.tuberculosis.

Algorithms↗

MREPATT: detection and analysis of exact consecutive repeats in genomic sequences.

UNLABELLED: We have developed a program to determine the number, length and position of exact consecutive repeats of short sequences in DNA fragments or whole genomes. The program also gives the statistical significance of results by comparing them with those expected for a random sequence generated according to a Markovian model. AVAILABILITY: MREPATT can be accessed on line at http://www.lsi.upc.es/~alggen under the RESEARCH and SEARCH links.

Algorithms↗

MUTAGEN: multi-user tool for annotating genomes.

SUMMARY: MUTAGEN is a free prokaryotic annotation system. It offers the advantages of genome comparison, graphical sequence browsers, search facilities and open-source for user-specific adjustments. The web-interface allows several users to access the system from standard desktop computers. The Sulfolobus acidocaldarius genome, and several plasmids and viruses have so far been analysed and annotated using MUTAGEN. AVAILABILITY: MUTAGEN is released as open-source software under GPL. The code is available for download and/or contribution at http://dac.molbio.ku.dk/bioinformatics/MUTAGEN/

Database Management Systems↗

A systematic method to identify genomic islands and its applications in analyzing the genomes of Corynebacterium glutamicum and Vibrio vulnificus CMCP6 chromosome I.

MOTIVATION: Some genomic islands contain horizontally transferred genes, which play critical roles in altering the genotypes and phenotypes of organisms, and horizontal gene transfer has been recognized as a universal event throughout bacterial evolution. A windowless method to display the distribution of genomic GC content, the cumulative GC profile, is proposed to identify genomic islands in genomes whose complete genome sequences are available. Two new indices are proposed to assess the codon usage bias and amino acid usage bias in genomic islands. RESULTS: A 211 kb genomic island (CGGI-1) has been identified in the genome of Corynebacterium glutamicum, and three genomic islands VVGI-1, VVGI-2 and VVGI-3, with lengths 167, 40 and 33 kb, respectively, have been identified in the genome of Vibrio vulnificus CMCP6 chromosome I. The CGGI-1 is flanked by two approximately 500 bp direct repeats, and utilizes a Val-tRNA as the integration site. For the VVGI-1 and VVGI-2, each has an integrase gene at 5' junction. All the identified genomic islands show unusual GC content, codon usage and amino acid usage, compared with the rest of the genomes. In addition, it is found that genomic islands are fairly homogenous in terms of GC content variation. An index, h, to quantify the homogeneity of GC content for genomic islands is proposed, and it is shown that h is less than 0.1 for all the genomic islands analyzed. The cumulative GC profile, as well as various indices to assess the codon usage bias, amino acid usage bias and homogeneity of the genomic islands, will be useful in the analysis of other genomes. AVAILABILITY: Programs used in this work and numerical results are available upon request.

Algorithms↗

In silico fine-mapping: narrowing disease-associated loci by intergenomics.

UNLABELLED: Genetic linkage and association studies define quantitative trait loci (QTLs) and susceptibility loci (SLs) that influence the phenotype of polygenic traits. A web-accessible application was created to identify intergenomic consensuses to fine map QTLs and SLs in silico and select particularly promising candidate genes for such traits. Furthermore, this approach offers an empirical evaluation of animal models for their applicability to the study of human traits. AVAILABILITY: http://qtl.pzr.uni-rostock.de/qtlmix.php CONTACT: serrano@pzr.uni-rostock.de.

Animals↗