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,135 records · Page 63Linked to original sources

Making sense of olfaction through predictions of the 3-D structure and function of olfactory receptors.

We used the MembStruk first principles computational technique to predict the three-dimensional (3-D) structure of six mouse olfactory receptors (S6, S18, S19, S25, S46 and S50) for which experimental odorant recognition profiles are available for a set of 24 odorants (4-9 carbons aliphatic alcohols, acids, bromo-acids and diacids). We used the HierDock method to scan each predicted OR structure for potential odorant binding site(s) and to calculate binding energies of each odorant in these binding sites. The calculated binding affinity profiles are in good agreement with experimental activation profiles, validating the predicted 3-D structures and the predicted binding sites. For each of the six ORs, the binding site is located between trans-membrane domains (TMs) 3-6, with contributions from extracellular loops 2 and 3. In particular, we find six residue positions in TM3 and TM6 to be consistently involved in the binding modes of the odorants. Indeed, the differences in the experimental recognition profiles can be explained on the basis of these critical residues alone. These predictions are also consistent with mutation data on ligand binding for catecholamine receptors and sequence hypervariability studies for ORs. Based on this analysis, we defined amino acid patterns associated with the recognition of short aliphatic alcohols and mono-acids. Using these two sequence fingerprints to probe the alignment of 869 OR sequences from the mouse genome, we identified 34 OR sequences matching the fingerprint for aliphatic mono-acids and 36 corresponding to the recognition pattern for aliphatic alcohols. We suggest that these two sets of ORs might function as basic arrays for uniquely recognizing aliphatic alcohols and acids. We screened a library of 89 additional molecules against the six ORs and found that this set of ORs is likely to respond to aldehydes and esters with longer carbon chains than their currently known agonists. We also find that compounds associated with the flavor in foods are often among the best calculated binding affinities. This suggests that physiologic ligands for these ORs may be found among aldehydes and esters associated with flavor.

Acids↗

Transposon mutagenesis of the mouse germline.

Sleeping Beauty is a synthetic "cut-and-paste" transposon of the Tc1/mariner class. The Sleeping Beauty transposase (SB) was constructed on the basis of a consensus sequence obtained from an alignment of 12 remnant elements cloned from the genomes of eight different fish species. Transposition of Sleeping Beauty elements has been observed in cultured cells, hepatocytes of adult mice, one-cell mouse embryos, and the germline of mice. SB has potential as a random germline insertional mutagen useful for in vivo gene trapping in mice. Previous work in our lab has demonstrated transposition in the male germline of mice and transmission of novel inserted transposons in offspring. To determine sequence preferences and mutagenicity of SB-mediated transposition, we cloned and analyzed 44 gene-trap transposon insertion sites from a panel of 30 mice. The distribution and sequence content flanking these cloned insertion sites was compared to 44 mock insertion sites randomly selected from the genome. We find that germline SB transposon insertion sites are AT-rich and the sequence ANNTANNT is favored compared to other TA dinucleotides. Local transposition occurs with insertions closely linked to the donor site roughly one-third of the time. We find that approximately 27% of the transposon insertions are in transcription units. Finally, we characterize an embryonic lethal mutation caused by endogenous splicing disruption in mice carrying a particular intron-inserted gene-trap transposon.

Animals↗

trome, trEST and trGEN: databases of predicted protein sequences.

We previously introduced two new protein databases (trEST and trGEN) of hypothetical protein sequences predicted from EST and HTG sequences, respectively. Here, we present the updates made on these two databases plus a new database (trome), which uses alignments of EST data to HTG or full genomes to generate virtual transcripts and coding sequences. This new database is of higher quality and since it contains the information in a much denser format it is of much smaller size. These new databases are in a Swiss-Prot-like format and are updated on a weekly basis (trEST and trGEN) or every 3 months (trome). They can be downloaded by anonymous ftp from ftp://ftp.isrec.isb-sib.ch/pub/databases.

Animals↗

Draft versus finished sequence data for DNA and protein diagnostic signature development.

Sequencing pathogen genomes is costly, demanding careful allocation of limited sequencing resources. We built a computational Sequencing Analysis Pipeline (SAP) to guide decisions regarding the amount of genomic sequencing necessary to develop high-quality diagnostic DNA and protein signatures. SAP uses simulations to estimate the number of target genomes and close phylogenetic relatives (near neighbors or NNs) to sequence. We use SAP to assess whether draft data are sufficient or finished sequencing is required using Marburg and variola virus sequences. Simulations indicate that intermediate to high-quality draft with error rates of 10(-3)-10(-5) (approximately 8x coverage) of target organisms is suitable for DNA signature prediction. Low-quality draft with error rates of approximately 1% (3x to 6x coverage) of target isolates is inadequate for DNA signature prediction, although low-quality draft of NNs is sufficient, as long as the target genomes are of high quality. For protein signature prediction, sequencing errors in target genomes substantially reduce the detection of amino acid sequence conservation, even if the draft is of high quality. In summary, high-quality draft of target and low-quality draft of NNs appears to be a cost-effective investment for DNA signature prediction, but may lead to underestimation of predicted protein signatures.

Computational Biology↗

Detecting non-orthology in the COGs database and other approaches grouping orthologs using genome-specific best hits.

Correct orthology assignment is a critical prerequisite of numerous comparative genomics procedures, such as function prediction, construction of phylogenetic species trees and genome rearrangement analysis. We present an algorithm for the detection of non-orthologs that arise by mistake in current orthology classification methods based on genome-specific best hits, such as the COGs database. The algorithm works with pairwise distance estimates, rather than computationally expensive and error-prone tree-building methods. The accuracy of the algorithm is evaluated through verification of the distribution of predicted cases, case-by-case phylogenetic analysis and comparisons with predictions from other projects using independent methods. Our results show that a very significant fraction of the COG groups include non-orthologs: using conservative parameters, the algorithm detects non-orthology in a third of all COG groups. Consequently, sequence analysis sensitive to correct orthology assignments will greatly benefit from these findings.

Algorithms↗

Pristionchus.org: a genome-centric database of the nematode satellite species Pristionchus pacificus.

Comparative studies have been of invaluable importance to the understanding of evolutionary biology. The evolution of developmental programs can be studied in nematodes at a single cell resolution given their fixed cell lineage. We have established Pristionchus pacificus as a major satellite organism for evolutionary developmental biology relative to Caenorhabditis elegans, the model nematode. Online genomic information to support studies in this satellite system can be accessed at http://www.pristionchus.org. Our web resource offers diverse content covering genome browsing, genetic and physical maps, similarity searches, a community platform and assembly details. Content will be continuously improved as we annotate the P.pacificus genome, and will be an indispensable resource for P.pacificus genomics.

Animals↗

Class-level relationships in the phylum Cnidaria: molecular and morphological evidence.

The evolutionary history of cnidarian life cycles has been debated since the 1880s, with different hypotheses favored even by current textbooks. Contributing to the disagreement is the fact that the systematic relationships of the four cnidarian classes have received relatively little examination using modern systematic methods. Here we present analyses of class-level relationships based on 18S ribosomal DNA (rDNA) sequence, mitochondrial 16S rDNA sequence, mitochondrial genome structure, and morphological characters. DNA sequences were aligned using a repeatable parsimony-based approach incorporating a range of alignment parameters. Analyses of individual data sets and of all data combined are unanimous in grouping the classes possessing a medusa stage, leaving the holobenthic Anthozoa basal within the phylum.

Animals↗

Phylogenetic analysis of gag genes from 70 international HIV-1 isolates provides evidence for multiple genotypes.

OBJECTIVE: To determine the extent of genetic variation among internationally collected HIV-1 isolates, to analyse phylogenetic relationships and the geographic distribution of different variants. DESIGN: Phylogenetic comparison of 70 HIV-1 isolates collected in 15 countries on four continents. METHODS: To sequence the complete gag genome of HIV-1 isolates, build multiple sequence alignments and construct phylogenetic trees using distance matrix methods and maximum parsimony algorithms. RESULTS: Phylogenetic tree analysis identified seven distinct genotypes. The seven genotypes were evident by both distance matrix methods and maximum parsimony analysis, and were strongly supported by bootstrap resampling of the data. The intra-genotypic gag distances averaged 7%, whereas the inter-genotypic distances averaged 14%. The geographic distribution of variants was complex. Some genotypes have apparently migrated to several continents and many areas harbor a mixture of genotypes. Related variants may cluster in certain areas, particularly isolates from a single city collected over a short time. CONCLUSIONS: The genetic variation among HIV-1 isolates is more extensive than previously appreciated. At least seven distinct HIV-1 genotypes can be identified. Diversification, migration and establishment of local, temporal 'blooms' of particular variants may all occur concomitantly.

Africa↗

Comparative genomics provides evidence for an ancient genome duplication event in fish.

There are approximately 25 000 species in the division Teleostei and most are believed to have arisen during a relatively short period of time ca. 200 Myr ago. The discovery of 'extra' Hox gene clusters in zebrafish (Danio rerio), medaka (Oryzias latipes), and pufferfish (Fugu rubripes), has led to the hypothesis that genome duplication provided the genetic raw material necessary for the teleost radiation. We identified 27 groups of orthologous genes which included one gene from man, mouse and chicken, one or two genes from tetraploid Xenopus and two genes from zebrafish. A genome duplication in the ancestor of teleost fishes is the most parsimonious explanation for the observations that for 15 of these genes, the two zebrafish orthologues are sister sequences in phylogenies that otherwise match the expected organismal tree, the zebrafish gene pairs appear to have been formed at approximately the same time, and are unlinked. Phylogenies of nine genes differ a little from the tree predicted by the fish-specific genome duplication hypothesis: one tree shows a sister sequence relationship for the zebrafish genes but differs slightly from the expected organismal tree and in eight trees, one zebrafish gene is the sister sequence to a clade which includes the second zebrafish gene and orthologues from Xenopus, chicken, mouse and man. For these nine gene trees, deviations from the predictions of the fish-specific genome duplication hypothesis are poorly supported. The two zebrafish orthologues for each of the three remaining genes are tightly linked and are, therefore, unlikely to have been formed during a genome duplication event. We estimated that the unlinked duplicated zebrafish genes are between 300 and 450 Myr. Thus, genome duplication could have provided the genetic raw material for teleost radiation. Alternatively, the loss of different duplicates in different populations (i.e. 'divergent resolution') may have promoted speciation in ancient teleost populations.

Animals↗

A phylogenomic approach to bacterial phylogeny: evidence of a core of genes sharing a common history.

It has been claimed that complete genome sequences would clarify phylogenetic relationships between organisms, but up to now, no satisfying approach has been proposed to use efficiently these data. For instance, if the coding of presence or absence of genes in complete genomes gives interesting results, it does not take into account the phylogenetic information contained in sequences and ignores hidden paralogies by using a BLAST reciprocal best hit definition of orthology. In addition, concatenation of sequences of different genes as well as building of consensus trees only consider the few genes that are shared among all organisms. Here we present an attempt to use a supertree method to build the phylogenetic tree of 45 organisms, with special focus on bacterial phylogeny. This led us to perform a phylogenetic study of congruence of tree topologies, which allows the identification of a core of genes supporting similar species phylogeny. We then used this core of genes to infer a tree. This phylogeny presents several differences with the rRNA phylogeny, notably for the position of hyperthermophilic bacteria.

Computational Biology↗

Construction and characterization of bacterial artificial chromosome library of black-handed spider monkey (Ateles geoffroyi).

The large-insert genomic DNA library is a critical resource for genome-wide genetic dissection of target species. We constructed a high-redundancy bacterial artificial chromosome (BAC) library of a New World monkey species, the black-handed spider monkey (Ateles geoffroyi). A total of 193 152 BAC clones were generated in this library. The average insert size of the BAC clones was estimated to be 184.6 kb with the small inserts (50-100 kb) accounting for less than 3% and the non-recombinant clones only 1.2%. Assuming a similar genome size with humans, the spider monkey BAC library has about 11x genome coverage. In addition, by end sequencing of randomly selected BAC clones, we generated 367 sequence tags for the library. When blasted against human genome, they showed a good correlation between the number of hit clones and the size of the chromosomes, an indication of unbiased chromosomal distribution of the library. This black-handed spider monkey BAC library would serve as a valuable resource in comparative genomic study and large-scale genome sequencing of nonhuman primates.

Animals↗

Wheat EST sequence assembly facilitates comparison of gene contents among plant species and discovery of novel genes.

Using a strategy requiring only modest computational resources, wheat expressed sequence tag (EST) sequences from various sources were assembled into contigs and compared with a nonredundant barley sequence assembly, with ESTs, with complete draft genome sequences of rice and Arabidopsis thaliana, and with ESTs from other plant species. These comparisons indicate that (i) wheat sequences available from public sources represent a substantial proportion of the diversity of wheat coding sequences, (ii) prediction of open reading frames in the whole genome sequence improves when supplemented with EST information from other species, (iii) a substantial number of candidates for novel genes that are unique to wheat or related species can be identified, and (iv) a smaller number of genes can be identified that are common to monocots and dicots but absent from Arabidopsis. The sequences in the last group may have been lost from Arabidopsis after descendance from a common ancestor. Examples of potential novel wheat genes and Triticeae-specific genes are presented.

Arabidopsis↗

Semantic similarity measures as tools for exploring the gene ontology.

Many bioinformatics resources hold data in the form of sequences. Often this sequence data is associated with a large amount of annotation. In many cases this data has been hard to model, and has been represented as scientific natural language, which is not readily computationally amenable. The development of the Gene Ontology provides us with a more accessible representation of some of this data. However it is not clear how this data can best be searched, or queried. Recently we have adapted information content based measures for use with the Gene Ontology (GO). In this paper we present detailed investigation of the properties of these measures, and examine various properties of GO, which may have implications for its future design.

Classification↗

Similarities and differences of polyadenylation signals in human and fly.

BACKGROUND: Cleavage of messenger RNA (mRNA) precursors is an essential step in mRNA maturation. The signal recognized by the cleavage enzyme complex has been characterized as an A rich region upstream of the cleavage site containing a motif with consensus AAUAAA, followed by a U or UG rich region downstream of the cleavage site. RESULTS: We studied these signals using exhaustive databases of cleavage sites obtained from aligning raw expressed sequence tags (EST) sequences to genomic sequences in Homo sapiens and Drosophila melanogaster. These data show that the polyadenylation signal is highly conserved in human and fly. In addition, de novo motif searches generated a refined description of the U-rich downstream sequence (DSE) element, which shows more divergence between the two species. These refined motifs are applied, within a Hidden Markov Model (HMM) framework, to predict mRNA cleavage sites. CONCLUSION: We demonstrate that the DSE is a specific motif in both human and Drosophila. These findings shed light on the sequence correlates of a highly conserved biological process, and improve in silico prediction of 3' mRNA cleavage and polyadenylation sites.

3' Untranslated Regions↗

Methods for identifying and mapping recent segmental and gene duplications in eukaryotic genomes.

The aim of this chapter is to provide instruction for analyzing and mapping recent segmental and gene duplications in eukaryotic genomes. We describe a bioinformatics-based approach utilizing computational tools to manage eukaryotic genome sequences to characterize and understand the evolutionary fates and trajectories of duplicated genes. An introduction to bioinformatics tools and programs such as BLAST, Perl, BioPerl, and the GFF specification provides the necessary background to complete this analysis for any eukaryotic genome of interest.

Animals↗

Group 3 chromosome bin maps of wheat and their relationship to rice chromosome 1.

The focus of this study was to analyze the content, distribution, and comparative genome relationships of 996 chromosome bin-mapped expressed sequence tags (ESTs) accounting for 2266 restriction fragments (loci) on the homoeologous group 3 chromosomes of hexaploid wheat (Triticum aestivum L.). Of these loci, 634, 884, and 748 were mapped on chromosomes 3A, 3B, and 3D, respectively. The individual chromosome bin maps revealed bins with a high density of mapped ESTs in the distal region and bins of low density in the proximal region of the chromosome arms, with the exception of 3DS and 3DL. These distributions were more localized on the higher-resolution group 3 consensus map with intermediate regions of high-mapped-EST density on both chromosome arms. Gene ontology (GO) classification of mapped ESTs was not significantly different for homoeologous group 3 chromosomes compared to the other groups. A combined analysis of the individual bin maps using 537 of the mapped ESTs revealed rearrangements between the group 3 chromosomes. Approximately 232 (44%) of the consensus mapped ESTs matched sequences on rice chromosome 1 and revealed large- and small-scale differences in gene order. Of the group 3 mapped EST unigenes approximately 21 and 32% matched the Arabidopsis coding regions and proteins, respectively, but no chromosome-level gene order conservation was detected.

Chromosome Mapping↗

A 2500-locus bin map of wheat homoeologous group 5 provides insights on gene distribution and colinearity with rice.

We constructed high-density deletion bin maps of wheat chromosomes 5A, 5B, and 5D, including 2338 loci mapped with 1052 EST probes and 217 previously mapped loci (total 2555 loci). This information was combined to construct a consensus chromosome bin map of group 5 including 24 bins. A relatively higher number of loci were mapped on chromosome 5B (38%) compared to 5A (34%) and 5D (28%). Differences in the levels of polymorphism among the three chromosomes were partially responsible for these differences. A higher number of duplicated loci was found on chromosome 5B (42%). Three times more loci were mapped on the long arms than on the short arms, and a significantly higher number of probes, loci, and duplicated loci were mapped on the distal halves than on the proximal halves of the chromosome arms. Good overall colinearity was observed among the three homoeologous group 5 chromosomes, except for the previously known 5AL/4AL translocation and a putative small pericentric inversion in chromosome 5A. Statistically significant colinearity was observed between low-copy-number ESTs from wheat homoeologous group 5 and rice chromosomes 12 (88 ESTs), 9 (72 ESTs), and 3 (84 ESTs).

Chromosome Mapping↗