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

Sequence alignments, variabilities, and vagaries.

It seems as if the algorithms and weighting matrices for multiple sequence alignments of the highly divergent members of the P450 gene superfamily have advanced to the point that unknown proteins can be aligned to structurally known members with reasonable accuracy. As stated earlier, the alignment tends to break down at gaps in the sequence alignments, but these regions can be improved manually. This type of alignment and analysis is especially useful for extracting and analyzing the various genome databases. Variations of the conservation analysis can be used to identify charged and uncharged residues that may be important in domain/domain interactions with redox partners or effector molecules (e.g., cytochrome b5). From these alignments and with comparative analysis within families and across P450 families, one can readily obtain an estimation of those residues that might be involved in substrate binding, in redox partner interaction, and in the catalytic mechanism.

Amino Acid Sequence↗
PubMed2001

Catarrhine phylogeny: noncoding DNA evidence for a diphyletic origin of the mangabeys and for a human-chimpanzee clade.

Maximum-parsimony and maximum-likelihood analyses of two of the serum albumin gene's intron sequences from 24 catarrhines (17 cercopithecid and 7 hominid) and 3 platyrrhines (an outgroup to the catarrhines) yielded results on catarrhine phylogeny that are congruent with those obtained with noncoding sequences of the gamma(1)-gamma(2) globin gene genomic region, using only those flanking and intergenic gamma sequences that in their history were not involved in gene conversion. A data set that combined in a tandem alignment these two sets of noncoding DNA orthologues from the two unlinked nuclear genomic loci yielded the following confirmatory results both on the course of cladistic branchings (the divisions in a cladistic classification of higher ranking taxa into subordinate taxa) and on the ages of the taxa (each taxon representing a clade). The cercopithecid branch of catarrhines, at approximately 14 Ma (mega annum) divided into Colobini (the leaf-eating Old World monkeys) and Cercopithecini (the cheek-pouched Old World monkeys). At approximately 10-9 Ma, Colobini divided into an African clade, Colobina, and an Asian clade, Presbytina; similarly at this time level, Cercopithecini divided into Cercopithecina (the guenons, patas, and green monkeys) and Papionina. At approximately 7 Ma, Papionina divided into Macaca, Cercocebus, and Papio. At approximately 5 Ma, Cercocebus divided subgenerically into C. (Cercocebus) for terrestrial mangabeys and C. (Mandrillus) for drills and mandrills, while at approximately 4 Ma Papio divided subgenerically into P. (Locophocebus) for arboreal mangabeys, P. (Theropithecus) for gelada baboons, and P. (Papio) for hamadryas baboons. In turn, the hominid branch of catarrhines at approximately 18 Ma divided into Hylobatini (gibbons and siamangs) and Hominini; at approximately 14 Ma, Hominini divided into Pongina (orangutans) and Hominina; at approximately 7 Ma, Hominina divided into Gorilla and Homo; and at approximately 6-5 Ma, Homo divided subgenerically into H. (Homo) for humans and H. (Pan) for common and bonobo chimpanzees. Rates of noncoding DNA evolution were assessed using a data set of noncoding gamma sequence orthologues that represented 18 catarrhines, 16 platyrrhines, 3 non-anthropoid primates (2 tarsiers and 1 strepsirhine), and rabbit (as outgroup to the primates). Results obtained with this data set revealed a faster rate of nucleotide substitutions in the early primate lineage to the anthropoid (platyrrhine/catarrhine) ancestor than from that ancestor to the present. Rates were slower in catarrhines than in platyrrhines, slower in the cheek-pouched than in the leaf-eating cercopithecids, and slower yet in the hominids. On relating these results to data on brain sizes and life spans, it was suggested that life-history strategies that favor intelligence and longer life spans also select for decreases in de novo mutation rates.

Albumins↗
PubMed1995

Interspecific reassortment of genomic segments in the evolution of cucumoviruses.

Segmented genomes of RNA viruses are thought to evolve and be maintained in analogy to sexual recombination and reassortment in eukaryotic systems. If reassortment among genomes is an important event in cucumoviral evolution, then such events should be detectable among extant viruses. In this study, phylogenetic analyses of cucumoviruses were performed using aligned amino acid sequences. The results reveal different relationships among species when the three genomic segments are compared, suggesting that reassortment events have given rise to extant forms. In addition, we describe a cucumovirus isolate that is composed of genomic segments from two distinct viral species. These results indicate that reassortment events may provide a mechanism for speciation in cucumoviruses.

Amino Acid Sequence↗
PubMed1999

Transproteomic evidence of a loop-deletion mechanism for enhancing protein thermostability.

Understanding the molecular determinants of protein thermostability is of theoretical and practical importance. While numerous determinants have been suggested, no molecular feature has been judged of paramount importance, with the possible exception of ion-pair networks. The difficulty in identifying the main determinants may have been the limited structural information available on the thermostable proteins. Recently the complete genomes for mesophilic, thermophilic and hyperthermophilic organisms have been sequenced, vastly improving the potential for uncovering general trends in sequence and structure evolution related to thermostability and, thus, for isolating the more important determinants. From a comparative analysis of 20 complete genomes, we find a trend towards shortened thermophilic proteins relative to their mesophilic homologs. Moreover, sequence alignments to proteins of known structure indicate that thermophilic sequences are more likely than their mesophilic homologs to have deletions in exposed loop regions. The new genomes offer enough comparable sequences to compute meaningful statistics that point to loop deletion as a general evolutionary strategy for increasing thermostability.

Archaeal Proteins↗
PubMed1999

The Yeast Proteome Database (YPD): a model for the organization and presentation of genome-wide functional data.

The Yeast Proteome Database (YPD) is a model for the organization and presentation of comprehensive protein information. Based on the detailed curation of the scientific literature for the yeast Saccharomyces cerevisiae, YPD contains more than 50 000 annotations lines derived from the review of 8500 research publications. The information concerning each of the approximately 6100 yeast proteins is structured around a convenient one-page format, the Yeast Protein Report, with additional information provided as pop-up windows. Protein classification schema have been revised this year, defining each protein's cellular role, function and pathway, and adding a Functional to the Yeast Protein Report. These changes provide the user with a succinct summary of the protein's function and its place in the biology of the cell, and they enhance the power of YPD Search functions. Precalculated sequence alignments have been added, to provide a crossover point for comparative genomics. The first transcript profiling data has been integrated into the YPD Protein Reports, providing the framework for the presentation of genome-wide functional data. The Yeast Proteome Database can be accessed on the Web at http://www.proteome.com/YPDhome.html

Computational Biology↗
PubMed1999

From global expression data to gene networks.

Allowing the parallel monitoring of the transcription of thousands of genes, microarrays constitute a powerful technique for functional genomics. In a recent paper, a clustering method and a local alignment software were combined to identify DNA motifs in sets of yeast genes endowed with similar transcription profiles throughout mitosis (1). Identifying various known transcriptional binding sites together with new putative ones, the authors made a significant step towards a systematic characterization of the regulatory structure of genomic networks. BioEssays 1999;21:895-899.

Biosensing Techniques↗
PubMed2006

Detecting disease-causing mutations in the human genome by haplotype matching.

Comparisons between haplotypes from affected patients and the human reference genome are frequently used to identify candidates for disease-causing mutations, even though these alignments are expected to reveal a high level of background neutral polymorphism. This limits the scope of genetic studies to relatively small genomic intervals, because current methods for distinguishing potential causal mutations from neutral variation are inefficient. Here we describe a new strategy for detecting mutations that is based on comparing affected haplotypes with closely matched control sequences from healthy individuals, rather than with the human reference genome. We use theory, simulation, and a real data set to show that this approach is expected to reduce the number of sequence variants that must be subjected to follow-up analysis by at least a factor of 20 when closely matched control sequences are selected from a reference panel with as few as 100 control genomes. We also define a reference data resource that would allow efficient application of this strategy to large critical intervals across the genome.

Alleles↗
PubMed1994

Isolation, characterization, and disruption of the yeast gene encoding cytosolic NADP-specific isocitrate dehydrogenase.

The cytosolic isozyme of NADP-specific isocitrate dehydrogenase (IDP2) was purified from a Saccharomyces cerevisiae mutant containing a chromosomal disruption in the gene encoding the mitochondrial isozyme (IDP1). IDP2 was shown to be a homodimer with a subunit molecular weight of approximately 45,000 and an isoelectric point of 5.5. Amino acid sequences were obtained for tryptic peptides of IDP2 and used to plan polymerase chain reactions. A resulting 400 bp DNA fragment was used as a hybridization probe to isolate the IDP2 gene from a yeast genomic DNA library. The complete nucleotide sequence of the IDP2 coding region was determined and translated into a 412-residue amino acid sequence. IDP2 and IDP1 were found to be identical in 71% of the aligned residue positions. The identity of the IDP2 gene was confirmed by genomic replacement with a disrupted IDP2 coding region. Haploid yeast strains lacking either or both IDP2 and IDP1 were constructed by genetic crosses of mutant strains containing disruptions in chromosomal IDP2 and IDP1 loci. No dramatic differences in growth rates with common carbon sources could be attributed to these disruptions.

Amino Acid Sequence↗
PubMed2005

Genome comparison of Pseudomonas aeruginosa large phages.

Pseudomonas aeruginosa phage EL is a dsDNA phage related to the giant phiKZ-like Myoviridae. The EL genome sequence comprises 211,215 bp and has 201 predicted open reading frames (ORFs). The EL genome does not share DNA sequence homology with other viruses and micro-organisms sequenced to date. However, one-third of the predicted EL gene products (gps) shares similarity (Blast alignments of 17-55% amino acid identity) with phiKZ proteins. Comparative EL and phiKZ genomics reveals that these giant phages are an example of substantially diverged genetic mosaics. Based on the position of similar EL and phiKZ predicted gene products, five genome regions can be delineated in EL, four of which are relatively conserved between EL and phiKZ. Region IV, a 17.7 kb genome region with 28 predicted ORFs, is unique to EL. Fourteen EL ORFs have been assigned a putative function based on protein similarity. Assigned proteins are involved in DNA replication and nucleotide metabolism (NAD+-dependent DNA ligase, ribonuclease HI, helicase, thymidylate kinase), host lysis and particle structure. EL-gp146 is the first chaperonin GroEL sequence identified in a viral genome. Besides a putative transposase, EL harbours predicted mobile endonucleases related to H-N-H and LAGLIDADG homing endonucleases associated with group I intron and intein intervening sequences.

Amino Acid Sequence↗
PubMed2002

Purification and partial characterization of an entomopoxvirus (DLEPV) from a parasitic wasp of tephritid fruit flies.

An insect poxvirus [entomopoxvirus (EPV)] occurs in the poison gland apparatus of female Diachasmimorpha longicaudata, a parasitic wasp of the Caribbean fruit fly, Anastrepha suspensa and other tephritid fruit flies. The DlEPV virion is 250-300 nm in diameter, has a "bumpy" appearance and a unipartite double stranded DNA genome of 290-300 kb. DlEPV DNA restriction fragment profiles differed from those reported for Amsacta moorei EPV (AmEPV) and Melanoplus sanguinipes EPV (MsEPV), the only two EPVs whose genomes have been sequenced, and from those reported for vaccinia (Vac), a vertebrate poxvirus (chordopoxvirus, ChPV). Blast search and ClustalW alignment of the amino acids deduced from the 2316 nucleotides of a DlEPV DNA fragment cloned from an EcoR1 genomic library revealed 75-78% homology with the putative DNA-directed RNA polymerases of AmEPV, MsEPV, and two ChPV homologs of the Vac J6R gene. Of the deduced 772 amino acids in the DlEPV sequence, 28.4% are conserved/substituted among the four poxviruses aligned, 12.9% occur in at least one EPV, 6.5% in at least one ChPV, 3.1% in at least one EPV and one ChPV, and 49.1% occur only in DlEPV. Although the RI-36-1 fragment represents a portion of the gene, it contains nucleotides that encode the NADFDGDE consensus sequence of known DNA-directed RNA polymerases. Western blots using a mouse polyclonal anti-DlEPV serum recognized six major protein bands in combined fractions of sucrose-purified DlEPV, at least one band in homogenates of male and female wasps, and at least two bands in host hemolymph that contained DlEPV virions. A digoxigenin-labeled DlEPV genomic DNA probe recognized DNA in dot-blots of male and female wasps. These results confirm that DlEPV is a true EPV and probably a member of the Group C EPVs. Unlike other EPVs, DlEPV does not express the spheroidin protein. Since it also replicates in both the wasp and fly, members of two different insect Orders, DlEPV may represent a new EPV Group, or a subgroup of the Group C viruses.

Amino Acid Sequence↗
PubMed2006

iCR: a web tool to identify conserved targets of a regulatory protein across the multiple related prokaryotic species.

Gene regulatory circuits are often commonly shared between two closely related organisms. Our web tool iCR (identify Conserved target of a Regulon) makes use of this fact and identify conserved targets of a regulatory protein. iCR is a special refined extension of our previous tool PredictRegulon- that predicts genome wide, the potential binding sites and target operons of a regulatory protein in a single user selected genome. Like PredictRegulon, the iCR accepts known binding sites of a regulatory protein as ungapped multiple sequence alignment and provides the potential binding sites. However important differences are that the user can select more than one genome at a time and the output reports the genes that are common in two or more species. In order to achieve this, iCR makes use of Cluster of Orthologous Group (COG) indices for the genes. This tool analyses the upstream region of all user-selected prokaryote genome and gives the output based on conservation target orthologs. iCR also reports the Functional class codes based on COG classification for the encoded proteins of downstream genes which helps user understand the nature of the co-regulated genes at the result page itself. iCR is freely accessible at http://www.cdfd.org.in/icr/.

Bacterial Proteins↗
PubMed1989

Sequencing of megabase plus DNA by hybridization: theory of the method.

A mismatch-free hybridization of oligonucleotides containing from 11 to 20 monomers to unknown DNA represents, in essence, a sequencing of a complementary target. Realizing this, we have used probability calculations and, in part, computer simulations to estimate the types and numbers of oligonucleotides that would have to be synthesized in order to sequence a megabase plus segment of DNA. We estimate that 95,000 specific mixes of 11-mers, mainly of the 5'(A,T,C,G)(A,T,C,G)N8(A,T,C,G)3' type, hybridized consecutively to dot blots of cloned genomic DNA fragments would provide primary data for the sequence assembly. An optimal mixture of representative libraries in M13 vector, having inserts of (i) 7 kb, (ii) 0.5 kb genomic fragments randomly ligated in up to 10-kb inserts, and (iii) tandem "jumping" fragments 100 kb apart in the genome, will be needed. To sequence each million base pairs of DNA, one would need hybridization data from about 2100 separate hybridization sample dots. Inevitable gaps and uncertainties in alignment of sequenced fragments arising from nonrandom or repetitive sequence organization of complex genomes and difficulties in cloning "poisonous" sequences in Escherichia coli, inherent to large sequencing by any method, have been considered and minimized by choice of libraries and number of subclones used for hybridization. Because it is based on simpler biochemical procedures, our method is inherently easier to automate than existing sequencing methods. The sequence can be derived from simple primary data only by extensive computing. Phased experimental tests and computer simulations increasing in complexity are needed before accurate estimates can be made in terms of cost and speed of sequencing by the new approach. Nevertheless, sequencing by hybridization should show advantages over existing methods because of the inherent redundancy and parallelism in its data gathering.

Algorithms↗
PubMed2004

Molecular modeling of protein function regions.

Experimental protein structures often provide extensive insight into the mode and specificity of small molecule binding, and this information is useful for understanding protein function and for the design of drugs. We have performed an analysis of the reliability with which ligand-binding information can be deduced from computer model structures, as opposed to experimentally derived ones. Models produced as part of the CASP experiments are used. The accuracy of contacts between protein model atoms and experimentally determined ligand atom positions is the main criterion. Only comparative models are included (i.e., models based on a sequence relationship between the protein of interest and a known structure). We find that, as expected, contact errors increase with decreasing sequence identity used as a basis for modeling. Analysis of the causes of errors shows that sequence alignment errors between model and experimental template have the most deleterious effect. In general, good, but not perfect, insight into ligand binding can be obtained from models based on a sequence relationship, providing there are no alignment errors in the model. The results support a structural genomics strategy based on experimental sampling of structure space so that all protein domains can be modeled on the basis of 30% or higher sequence identity.

Alpha-Globulins↗
PubMed2004

Fast and reliable sexing of prosimian and human DNA.

Molecular sexing of mammals is normally done by PCR amplification of Y chromosomal fragments, or coamplification of homologous fragments from both sex chromosomes. Existing primers are often unreliable for distantly related species due to mutations in primer regions. Currently there are no published primers for the sexing of prosimian DNA. We show that an existing method (using the zinc finger protein) based on a size difference between the X and Y homologs does not work in prosimians. Multiple alignments of distantly related mammalian species from Genbank and genome databases enabled us to identify conserved regions in the amelogenin gene. Using these conserved regions, we can target species that have no sequence information. We designed a single, conserved primer pair that is useful for fast and reliable molecular sexing of prosimian primates. A single PCR yields two fragments in males and only one in females, which are easily separated with the use of agarose gels. Amplification of separable fragments was successful in seven species of lemurs, as well as humans.

Amelogenin↗
PubMed2006

Degenerate PCR method for identification of an antiapoptotic gene in BHV-1.

To investigate on the hypothetical presence of an antiapoptotic gene, we utilized the CODEHOP (COnsensus-DEgenerate Hybrid Oligonucleotide Primers) strategy amplifying unknown sequences from a background of genomic (bovine herpesvirus type-1) BHV-1 DNA. An alignment of carboxyl-terminal domains belonging to three proteins encoded by gamma34.5, MyD116 and GADD34 genes, was carried out to design degenerate PCR primers in highly conserved regions. This allowed the amplification of a 110 bp fragment. This fragment was subjected to automatic sequencing and DNA sequence analysis revealed that its position resided between the nt 14363 and the nt 14438 in bovine herpesvirus type-1 (BHV-1) Cooper strain sharing an identity of 86% (UL14). Transient transfections showed that UL14 protein is efficient in protecting MDBK and K562 cells from sorbitol induced apoptosis. The protein's anti-apoptotic function may derive from its heat shock protein-like properties.

Amino Acid Motifs↗
PubMed2001

Fine-mapping and construction of a bovine contig spanning the ovine callipyge locus.

The callipyge (CLPG) gene was fine-mapped by linkage analysis to a 4.6-cM chromosome interval on distal ovine OAR18q, flanked by microsatellite markers IDVGA30 and OY3. The OAR18q linkage map and human HSA14q transcript map were aligned by genotyping two bovine-hamster whole-genome radiation hybrid panels with the microsatellite markers, as well as with sequences corresponding to HSA 14q genes. Using Type I loci mapping to the IDVGA30-OY3 interval as anchor points, we have constructed a 1.4-Mb bovine BAC contig containing the IDVGA30-OY3 interval. We demonstrate that the IDVGA30-OY3 interval spans approximately 770 kb and contains at least four genes: YY1, WARS, DLK1, and GTL2.

Animals↗
PubMed1989

Complete nucleotide sequence of the rabbit beta-like globin gene cluster. Analysis of intergenic sequences and comparison with the human beta-like globin gene cluster.

The nucleotide sequence of the entire beta-like globin gene cluster of rabbits has been determined. This sequence of a continuous stretch of 44.5 x 10(3) base-pairs (bp) starts about 6 x 10(3) bp upstream from epsilon (the 5'-most gene) and ends about 12 x 10(3) bp downstream from beta (the 3'-most gene). Analysis of the sequence reveals that: (1) the sequence is relatively A + T rich (about 60%); (2) regions with high G + C content are associated with OcC repeats, a short interspersed repeated DNA in rabbits; (3) the distribution of polypurines, polypyrimidines and alternating purine/pyrimidine tracts is not random within the cluster; (4) most open reading frames are associated with known globin coding regions, OcC repeats or long interspersed repeats (L1 repeats); (5) the most prominent open reading frames are found in the L1 repeats; (6) different strand asymmetries in base composition are associated with embyronic and adult genes as well as the tandem L1 repeats at the 3' end of the cluster; and (7) essentially all the repeats appear to have been inserted by a transposon mechanism. A comparison of the sequence with itself by a dot-plot analysis has revealed nine new members of the OcC family of repeats in addition to the six previously reported. The OcC repeats tend to be clustered, particularly in the epsilon-gamma and gamma-psi delta intergenic regions. Dot-plot comparisons between the rabbit and the human clusters have revealed extensive sequence matches. Homology starts about 6 x 10(3) bp 5' to epsilon or as far upstream as the rabbit sequence is available. It continues throughout the entire cluster and stops about 0.7 x 10(3) bp 3' to beta, at which point several repeats have inserted in both rabbits and humans. Throughout the gene cluster, the homology is interrupted mainly by insertions or deletions in either the rabbit or the human genome. Almost all of the insertions are of known short or long repeated DNAs. The positions of the insertions are different in the two gene clusters, which indicates that both short and long repeats have been transposing throughout the genome for the time since the mammalian radiation. An alignment of rabbit and human sequences allows the calculation of the substitution rate around epsilon. Sequences far removed from the gene are evolving at a rate equivalent to the pseudogene rate, although some short regions show an apparently higher rate.(ABSTRACT TRUNCATED AT 400 WORDS)

Animals↗
PubMed2004

Molecular cloning and characterization of LCRG1 a novel gene localized to the tumor suppressor locus D17S800-D17S930.

Laryngeal carcinoma (LC) is the most commonly diagnosed malignancy and recently the incidence of this disease has increased. By means of the mRNA differential display method we identified a cDNA, Laryngeal carcinoma related gene 1 (LCRG1) which has significantly reduced expression in 40% (12/30) of primary LCs and in 6 of 10 various cancer cell lines. Northern Blot analysis showed that LCRG1 was expressed more abundantly in human heart, pancreas, skeletal muscle, and in murine testis, liver, brain and heart than in other tissues. The cDNA sequence of this gene is identical to part of the sequence of PAC HCIT75G16 clone (GenBank accession No. AC003042) from the chromosome band 17q12-21.1 which is one of the most common loss of heterozygosity (LOH) regions involved in LC, prostate cancer, etc. This gene is composed of six exons and spans about 60 kb of genomic DNA with a 3.4 kb mature transcript. The alignment of this gene with STS markers localized the gene to the previously identified tumor-suppressor locus D17S800-D17S930. The putative protein encoded by this gene is 288 amino acid with one potential site for phosphorylation by casein kinase II and no significant homology to any known proteins in the public databases. The primary tumor suppressive functions (proliferation rate,soft agar growth and tumor formation) were observed in a LC cell line (Hep-2) by lipofectin transfection. Together these data strongly suggest a potential role of LCRG1 contributing to the development of LC.

Adaptor Proteins, Signal Transducing↗