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,621 records · Page 90Linked to original sources
PubMed2002

A genomic schism in birds revealed by phylogenetic analysis of DNA strings.

The molecular systematics of vertebrates has been based entirely on alignments of primary structures of macromolecules; however, higher order features of DNA sequences not used in traditional studies also contain valuable phylogenetic information. Recent molecular data sets conflict over the phylogenetic placement of flightless birds (ratites - paleognaths), but placement of this clade critically influences interpretation of character change in birds. To help resolve this issue, we applied a new bioinformatics approach to the largest molecular data set currently available. We distilled nearly one megabase (1 million base pairs) of heterogeneous avian genomic DNA from 20 birds and an alligator into genomic signatures, defined as the complete set of frequencies of short sequence motifs (strings), thereby providing a way to directly compare higher order features of nonhomologous DNA sequences. Phylogenetic analysis and principal component analysis of the signatures strongly support the traditional hypothesis of basal ratites and monophyly of the nonratite birds (neognaths) and imply that ratite genomes are linguistically primitive within birds, despite their base compositional similarity to neognath genomes. Our analyses show further that the phylogenetic signal of genomic signatures are strongest among deep splits within vertebrates. Despite clear problems with phylogenetic analysis of genomic signatures, our study raises intriguing issues about the biological and genomic differences that fundamentally differentiate paleognaths and neognaths.

Animals↗
PubMed1997

Analysis of the sequence diversity of the P1, HC, P3, NIb and CP genomic regions of several yam mosaic potyvirus isolates: implications for the intraspecies molecular diversity of potyviruses.

Partial sequences from serologically characterized yam mosaic potyvirus (YMV) isolates were determined in conserved (helper-component proteinase, HC; nuclear inclusion b, NIb) and variable (first protein, P1; third protein, P3; and coat protein, CP) regions of the potyviral genome in order to investigate the intraspecies molecular diversity of YMV. Multiple sequence alignments and pairwise comparisons were used to quantify the sequence polymorphism in these regions. Two levels of diversity were observed among YMV isolates: above 90% nucleotide (nt) sequence identities were found between YMV isolates of the same group (intragroup) regardless of the region considered, whereas identities between isolates from different groups (intergroup) were lower and depended upon the protein chosen. For instance, the average intergroup nt sequence identity between YMV isolates was about 65% in the P1 protein and the N terminus of the CP while there was more than 80% nt identity in the HC, P3 and NIb proteins. Thus P3 appeared to be conserved between YMV isolates even though this region was variable between potyvirus species. Similar analysis of the intraspecies molecular diversity of other potyviruses (potato virus Y, zucchini yellow mosaic virus, plum pox virus, pea seed-borne mosaic virus) led to the same results: (i) two levels of intraspecies molecular diversity were found (intragroup and intergroup); (ii) intraspecies molecular diversity differed from interspecies molecular diversity in the P3, P1 and N-terminal regions.

Amino Acid Sequence↗
PubMed2003

Mouse SCO-spondin, a gene of the thrombospondin type 1 repeat (TSR) superfamily expressed in the brain.

SCO-spondin is specifically expressed in the subcommissural organ (SCO), a secretory ependymal differentiation lining the roof of the third ventricular cavity of the brain. When released into the cerebro-spinal fluid (CSF), SCO-spondin aggregates and forms Reissner's fiber (RF), a structure present in the central canal of the spinal cord. SCO-spondin belongs to the superfamily of proteins exhibiting conserved motifs called TSRs for 'thrombospondin type 1 repeats' and involved in axonal pathfinding during development. The mouse SCO-spondin coding sequence was searched by alignement of the coding bovine SCO-spondin sequence with the mouse whole genome shotgun (WGS) supercontig (NW 000250). Compared to the bovine, mouse SCO-spondin shows 66.8% identity of amino acids. This extracellular matrix glycoprotein has a modular arrangement of several conserved domains including 25 TSRs, 10 low-density lipoprotein receptor (LDLr) type A repeats and cystein-rich regions in the -NH2 and -COOH ends. The spatio-temporal expression of SCO-spondin was analyzed using specific antisera and an homospecific SCO-spondin riboprobe. In the adult, the patterns obtained by in situ hybridization (ISH) and immunohistochemistry correlated well in the SCO, while Reissner's fiber and the ampulla caudalis were immunoreactive only. In the fetus, both the immuno and ISH reactions appeared between 14 and 15 days post coïtum (dpc) in the SCO anlage. In addition, the mouse SCO-spondin gene was located at chromosome 6, between marker D6Mit352 and D6Mit119, in a conserved syntenic region.

Amino Acid Sequence↗
PubMed2002

Singular value decomposition analysis of protein sequence alignment score data.

One of the standard tools for the analysis of data arranged in matrix form is singular value decomposition (SVD). Few applications to genomic data have been reported to date mainly for the analysis of gene expression microarray data. We review SVD properties, examine mathematical terms and assumptions implicit in the SVD formalism, and show that SVD can be applied to the analysis of matrices representing pairwise alignment scores between large sets of protein sequences. In particular, we illustrate SVD capabilities for data dimension reduction and for clustering protein sequences. A comparison is performed between SVD-generated clusters of proteins and annotation reported in the SWISS-PROT Database for a set of protein sequences forming the calycin superfamily, entailing all entries corresponding to the lipocalin, cytosolic fatty acid-binding protein, and avidin-streptavidin Prosite patterns.

Amino Acid Sequence↗
PubMed1997

The cox1 gene from Euglena gracilis: a protist mitochondrial gene without introns and genetic code modifications.

We present the nucleotide sequence of the cox 1 gene encoding subunit 1 of cytochrome c oxidase in Euglena gracilis, the first report on a mitochondrial gene from this protist. Its study reveals that the Euglena mitochondrial genome does not appear as a compact and homogeneous structure and that its A+T content is high (about 76%) whereas this value is less than 50% in nuclear DNA. The Euglena cox1 gene does not exhibit any intron, and an amino-acid alignment of Euglena COX1 with homologous proteins shows that the universal genetic code is used. Comparisons of the genomic and cDNA sequences of Euglena cox1 indicate that the transcript does not undergo RNA editing as found in trypanosomes and in higher plants. The phylogeny obtained with COX1 protein sequences is in agreement with that obtained with nuclear rRNA sequences and places Euglena and Trypanosoma far apart from other eukaryotes. This result strengthens the hypothesis that these protists represent the earliest mitochondrion-containing organisms.

Amino Acid Sequence↗
PubMed2006

Multiple sequence alignment with user-defined anchor points.

BACKGROUND: Automated software tools for multiple alignment often fail to produce biologically meaningful results. In such situations, expert knowledge can help to improve the quality of alignments. RESULTS: Herein, we describe a semi-automatic version of the alignment program DIALIGN that can take pre-defined constraints into account. It is possible for the user to specify parts of the sequences that are assumed to be homologous and should therefore be aligned to each other. Our software program can use these sites as anchor points by creating a multiple alignment respecting these constraints. This way, our alignment method can produce alignments that are biologically more meaningful than alignments produced by fully automated procedures. As a demonstration of how our method works, we apply our approach to genomic sequences around the Hox gene cluster and to a set of DNA-binding proteins. As a by-product, we obtain insights about the performance of the greedy algorithm that our program uses for multiple alignment and about the underlying objective function. This information will be useful for the further development of DIALIGN. The described alignment approach has been integrated into the TRACKER software system.

Journal Article↗
PubMed2004

PhyloGenie: automated phylome generation and analysis.

Phylogenetic reconstruction is the method of choice to determine the homologous relationships between sequences. Difficulties in producing high-quality alignments, which are the basis of good trees, and in automating the analysis of trees have unfortunately limited the use of phylogenetic reconstruction methods to individual genes or gene families. Due to the large number of sequences involved, phylogenetic analyses of proteomes preclude manual steps and therefore require a high degree of automation in sequence selection, alignment, phylogenetic inference and analysis of the resulting set of trees. We present a set of programs that automates the steps from seed sequence to phylogeny and a utility to extract all phylogenies that match specific topological constraints from a database of trees. Two example applications that show the type of questions that can be answered by phylome analysis are provided. The generation and analysis of the Thermoplasma acidophilum phylome with regard to lateral gene transfer between Thermoplasmata and Sulfolobus, showed best BLAST hits to be far less reliable indicators of lateral transfer than the corresponding protein phylogenies. The generation and analysis of the Danio rerio phylome provided more than twice as many proteins as described previously, supporting the hypothesis of an additional round of genome duplication in the actinopterygian lineage.

Amino Acid Sequence↗
PubMed1997

Evolutionary shifts in three major structural features of the mitochondrial genome among iguanian lizards.

A phylogenetic tree for major lineages of iguanian lizards is estimated from 1,488 aligned base positions (858 informative) of newly reported mitochondrial DNA sequences representing coding regions for eight tRNAs, ND2, and portions of ND1 and COI. Two well-supported groups are defined, the Acrodonta and the Iguanidae (sensu lato). This phylogenetic hypothesis is used to investigate evolutionary shifts in mitochondrial gene order, origin for light-strand replication, and secondary structure of tRNACys. These three characters shift together on the branch leading to acrodont lizards. Plate tectonics and the fossil record indicate that these characters changed in the Jurassic. We propose that changes to the secondary structure of tRNACys may destroy function of the origin for light-strand replication which, in turn, may facilitate shifts in gene order.

Animals↗
PubMed2005

Hemiascomycetous yeasts at the forefront of comparative genomics.

With more than a dozen species fully sequenced, as many as this partially sequenced, and more in development, yeasts are now used to explore the frontlines of comparative genomics of eukaryotes. Innovative procedures have been developed to compare and annotate genomes at various evolutionary distances, to identify short cis-acting regulatory elements, to map duplications, or to align syntenic blocks. Human and plant pathogens, in addition to yeasts that show a variety of interesting physiological properties, are included in this multidimensional comparative survey, which encompasses a very broad evolutionary range. As major steps of the evolutionary history of hemiascomycetous genomes emerge, precise questions on the general mechanisms of their evolution can be addressed, using both experimental and in silico methods.

Adaptation, Physiological↗
PubMed1999

Cloning of an ORF with homology to Mycobacterium echA1, encoding the enoyl-CoA hydratase, in Rhodococcus fascians.

An open reading frame encoding a polypeptide of significant homology (55.7% identity) with the enoyl-CoA hydratase encoded by the gene echA1 from Mycobacterium tuberculosis has been found in the genome of the plant-pathogen bacteria Rhodococcus fascians strain NRRL-B-15096. Sequence alignments showed that it possesses several conserved blocks common to E. coli, M. tuberculosis and human mitochondria. One of such blocks includes a glutamate residue located at position 149, corresponding to the glutamate 139 of Escherichia coli. This glutamate was previously shown to be the catalytic residue of enoyl-CoA hydratase in the multienzyme complex of fatty acid oxidation from E. coli. Our results provide additional information on the conserved domains of this enzyme. Significant homologies in other genome regions between R. fascians and M. tuberculosis confirm their phylogenetic relationship.

Amino Acid Sequence↗
PubMed1998

Gene structure and chromosomal localization of the human P2X7 receptor.

The genomic organization for the human P2X7 receptor gene was determined to comprise 13 exons. Alignment of the exon-intron junctions with those for the rat P2X2 gene demonstrated a precise conservation of the boundaries for the first 10 introns. The human P2X7 receptor gene was localized by in situ hybridization to chromosome 12q24. Radiation hybrid mapping indicated that this is within 130 kb of the gene for the homologous P2X4 receptor.

Animals↗
PubMed1999

The promise of comparative genomics in mammals.

Dense genetic maps of human, mouse, and rat genomes that are based on coding genes and on microsatellite and single-nucleotide polymorphism markers have been complemented by precise gene homolog alignment with moderate-resolution maps of livestock, companion animals, and additional mammal species. Comparative genetic assessment expands the utility of these maps in gene discovery, in functional genomics, and in tracking the evolutionary forces that sculpted the genome organization of modern mammalian species.

Animals↗
PubMed2002

Analysis of similarity within 142 pairs of orthologous intergenic regions of Caenorhabditis elegans and Caenorhabditis briggsae.

Patterns of similarity between genomes of related species reflect the distribution of selective constraint within DNA. We analyzed alignments of 142 orthologous intergenic regions of Caenorhabditis elegans and Caenorhabditis briggsae and found a mosaic pattern with regions of high similarity (phylogenetic footprints) interspersed with non-alignable sequences. Footprints cover approximately 20% of intergenic regions, often occur in clumps and are rare within 5' UTRs but common within 3' UTRs. The footprints have a higher ratio of transitions to transversions than expected at random and a higher GC content than the rest of the intergenic region. The number of footprints and the GC content of footprints within an intergenic region are higher when genes are oriented so that their 5' ends form the boundaries of the intergenic region. Overall, the patterns and characteristics identified here, along with other comparative and experimental studies, suggest that many footprints have a regulatory function, although other types of function are also possible. These conclusions may be quite general across eukaryotes, and the characteristics of conserved regulatory elements determined from genomic comparisons can be useful in prediction of regulation sites within individual DNA sequences.

Animals↗
PubMed2004

1001 optimal PDB structure alignments: integer programming methods for finding the maximum contact map overlap.

Protein structure comparison is a fundamental problem for structural genomics, with applications to drug design, fold prediction, protein clustering, and evolutionary studies. Despite its importance, there are very few rigorous methods and widely accepted similarity measures known for this problem. In this paper we describe the last few years of developments on the study of an emerging measure, the contact map overlap (CMO), for protein structure comparison. A contact map is a list of pairs of residues which lie in three-dimensional proximity in the protein's native fold. Although this measure is in principle computationally hard to optimize, we show how it can in fact be computed with great accuracy for related proteins by integer linear programming techniques. These methods have the advantage of providing certificates of near-optimality by means of upper bounds to the optimal alignment value. We also illustrate effective heuristics, such as local search and genetic algorithms. We were able to obtain for the first time optimal alignments for large similar proteins (about 1,000 residues and 2,000 contacts) and used the CMO measure to cluster proteins in families. The clusters obtained were compared to SCOP classification in order to validate the measure. Extensive computational experiments showed that alignments which are off by at most 10% from the optimal value can be computed in a short time. Further experiments showed how this measure reacts to the choice of the threshold defining a contact and how to choose this threshold in a sensible way.

Algorithms↗
PubMed2005

An Eulerian path approach to local multiple alignment for DNA sequences.

Expensive computation in handling a large number of sequences limits the application of local multiple sequence alignment. We present an Eulerian path approach to local multiple alignment for DNA sequences. The computational time and memory usage of this approach is approximately linear to the total size of sequences analyzed; hence, it can handle thousands of sequences or millions of letters simultaneously. By constructing a De Bruijn graph, most of the conserved segments are amplified as heavy Eulerian paths in the graph, and the original patterns distributed in sequences are recovered even if they do not exist in any single sequence. This approach can accurately detect unknown conserved regions, for both short and long, conserved and degenerate patterns. We further present a Poisson heuristic to estimate the significance of a local multiple alignment. The performance of our method is demonstrated by finding Alu repeats in the human genome. We compare the results with Alus marked by repeatmasker, where the two programs are in good agreement. Our method is robust under various conditions and superior to other methods in terms of efficiency and accuracy.

Base Sequence↗
PubMed2003

Identifying related L1 retrotransposons by analyzing 3' transduced sequences.

BACKGROUND: A large fraction of the human genome is attributable to L1 retrotransposon sequences. Not only do L1s themselves make up a significant portion of the genome, but L1-encoded proteins are thought to be responsible for the transposition of other repetitive elements and processed pseudogenes. In addition, L1s can mobilize non-L1, 3'-flanking DNA in a process called 3' transduction. Using computational methods, we collected DNA sequences from the human genome for which we have high confidence of their mobilization through L1-mediated 3' transduction. RESULTS: The precursors of L1s with transduced sequence can often be identified, allowing us to reconstruct L1 element families in which a single parent L1 element begot many progeny L1s. Of the L1s exhibiting a sequence structure consistent with 3' transduction (L1 with transduction-derived sequence, L1-TD), the vast majority were located in duplicated regions of the genome and thus did not necessarily represent unique insertion events. Of the remaining L1-TDs, some lack a clear polyadenylation signal, but the alignment between the parent-progeny sequences nevertheless ends in an A-rich tract of DNA. CONCLUSIONS: Sequence data suggest that during the integration into the genome of RNA representing an L1-TD, reverse transcription may be primed internally at A-rich sequences that lie downstream of the L1 3' untranslated region. The occurrence of L1-mediated transduction in the human genome may be less frequent than previously thought, and an accurate estimate is confounded by the frequent occurrence of segmental genomic duplications.

Base Sequence↗
PubMed2005

[Bioinformatic analysis of the 14-3-3 gene family in rice].

Using two-step HMM (hidden markov model) scan strategy,eight 14-3-3-like proteins were identified by searching the Oryza sativa L. ssp. japonica protein database. From them four genes were newly detected in this study. We confined the genes expressing in Nipponbare by EST search. Expression analysis also showed each gene expressed diversely within any individual,this tends suggested specific function of particular gene. Alignment of amino acid sequences suggested that there could be isoform function of the specific residues. The analyses of gene structure and chromosome location indicated that rice genome contains both epsilon and no-epsilon forms of 14-3-3 proteins. In addition,we analyzed the evolution of the rice 14-3-3 protein family.

14-3-3 Proteins↗
PubMed2004

Polymorphic segmental duplications at 8p23.1 challenge the determination of individual defensin gene repertoires and the assembly of a contiguous human reference sequence.

BACKGROUND: Defensins are important components of innate immunity to combat bacterial and viral infections, and can even elicit antitumor responses. Clusters of defensin (DEF) genes are located in a 2 Mb range of the human chromosome 8p23.1. This DEF locus, however, represents one of the regions in the euchromatic part of the final human genome sequence which contains segmental duplications, and recalcitrant gaps indicating high structural dynamics. RESULTS: We find that inter- and intraindividual genetic variations within this locus prevent a correct automatic assembly of the human reference genome (NCBI Build 34) which currently even contains misassemblies. Manual clone-by-clone alignment and gene annotation as well as repeat and SNP/haplotype analyses result in an alternative alignment significantly improving the DEF locus representation. Our assembly better reflects the experimentally verified variability of DEF gene and DEF cluster copy numbers. It contains an additional DEF cluster which we propose to reside between two already known clusters. Furthermore, manual annotation revealed a novel DEF gene and several pseudogenes expanding the hitherto known DEF repertoire. Analyses of BAC and working draft sequences of the chimpanzee indicates that its DEF region is also complex as in humans and DEF genes and a cluster are multiplied. Comparative analysis of human and chimpanzee DEF genes identified differences affecting the protein structure. Whether this might contribute to differences in disease susceptibility between man and ape remains to be solved. For the determination of individual DEF gene repertoires we provide a molecular approach based on DEF haplotypes. CONCLUSIONS: Complexity and variability seem to be essential genomic features of the human DEF locus at 8p23.1 and provides an ongoing challenge for the best possible representation in the human reference sequence. Dissection of paralogous sequence variations, duplicon SNPs ans multisite variations as well as haplotypes by sequencing based methods is the way for future studies of interindividual DEF locus variability and its disease association.

Chromosome Mapping↗