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

FootPrinter3: phylogenetic footprinting in partially alignable sequences.

FootPrinter3 is a web server for predicting transcription factor binding sites by using phylogenetic footprinting. Until now, phylogenetic footprinting approaches have been based either on multiple alignment analysis (e.g. PhyloVista, PhastCons), or on motif-discovery algorithms (e.g. FootPrinter2). FootPrinter3 integrates these two approaches, making use of local multiple sequence alignment blocks when those are available and reliable, but also allowing finding motifs in unalignable regions. The result is a set of predictions that joins the advantages of alignment-based methods (good specificity) to those of motif-based methods (good sensitivity, even in the presence of highly diverged species). FootPrinter3 is thus a tool of choice to exploit the wealth of vertebrate genomes being sequenced, as it allows taking full advantage of the sequences of highly diverged species (e.g. chicken, zebrafish), as well as those of more closely related species (e.g. mammals). The FootPrinter3 web server is available at: http://www.mcb.mcgill.ca/~blanchem/FootPrinter3.

Animals↗

Cerebellar gene expression profiling and eQTL analysis in inbred mouse strains selected for ethanol sensitivity.

BACKGROUND: Inbred Long-Sleep (ILS) and Inbred Short-Sleep (ISS) mice exhibit striking differences in a number of alcohol and drug related behaviors. This study examined the expression levels of more than 39,000 transcripts in these strains in the cerebellum, a major target of ethanol's actions in the CNS, to find differentially expressed (DE) candidate genes for these phenotypes. METHODS: Genes that were differentially expressed between the strains were identified using oligonucleotide arrays as well as complimentary DNA arrays. Sequence alignment was used to locate DE genes in the mouse genome assembly. In silico expression QTL (eQTL) mapping was used to identify chromosomal regions likely to control the transcription level of DE genes, and the EASE program identified overrepresented functional themes. The genomic region immediately upstream of the cyclase associated protein homolog 1 (Cap1) gene was directly sequenced from PCR products. RESULTS: Nearly 300 genes were identified as differentially expressed between the cerebella of ILS and ISS. These genes and their corresponding eQTLs map to genomic regions linked to several phenotypes that differ between the ILS and ISS strains, including ethanol preference and cocaine-induced locomotor activation on Chromosomes 4 and 7 respectively. Eight genes were cross-platform validated, four of which are more highly expressed in ILS cerebellum. Three SNPs, one of which disrupts a predicted Sp1 binding site, were found in the upstream region of Cap1, a strong candidate for influencing ethanol phenotypes. CONCLUSIONS: Many of these DE genes are candidates to influence ethanol and drug regulated phenotypes because they either map to ethanol related QTLs in the genome or are linked to them through eQTL mapping. Genes involved in calcium ion binding and transcriptional regulation are overrepresented and therefore these gene classes may influence ethanol behaviors in mice and humans.

Adenylyl Cyclases↗

Identification and functional analysis of a putative non-hr origin of DNA replication from the Spodoptera littoralis type B multinucleocapsid nucleopolyhedrovirus.

A putative non-hr origin of DNA replication was identified in the Spodoptera littoralis multinucleocapsid nucleopolyhedrovirus (SpliNPV) genome by transient replication assays. The putative SpliNPV ori was mapped to the PstI-J fragment between 75.1-77.9 map units in the SpliNPV genome. While the DNA sequence of the putative SpliNPV ori aligned with regions within the non-hr oris of Autographa californica, Orgyia pseudotsugata and Spodoptera exigua multinucleocapsid nucleopolyhedroviruses, it has limited DNA sequence identity with these elements. The sequence of the putative SpliNPV non-hr ori fragment contains a unique distribution of imperfect palindromes, multiple direct repeats and putative transcription factor-binding sites. Transient expression assays indicated that the putative SpliNPV ori fragment repressed SpliNPV lef-3 promoter-mediated luciferase reporter gene expression. However, the putative SpliNPV ori fragment itself was capable of directing luciferase expression in the absence of a recognizable baculovirus promoter element in an orientation-independent fashion, suggesting that DNA sequence motifs within its sequence can activate transcription. Gel mobility shift analyses confirmed that proteins within nuclear extracts from both uninfected and virus-infected cells bound with specificity to the putative SpliNPV ori fragment.

Animals↗

Genetic diversity of small-ruminant lentiviruses: characterization of Norwegian isolates of Caprine arthritis encephalitis virus.

Small-ruminant lentiviruses (SRLVs), including Caprine arthritis encephalitis virus (CAEV) in goats and maedi-visna virus (MVV) in sheep, are lentiviruses that, despite overall similarities, show considerable genetic variation in regions of the SRLV genome. To gain further knowledge about the genetic diversity and phylogenetic relationships among field isolates of SRLVs occurring in geographically distinct areas, the full-length genomic sequence of a CAEV isolate (CAEV-1GA) and partial env sequences obtained from Norwegian CAEV-infected goats were determined. The genome of CAEV-1GA consisted of 8,919 bp. Alignment studies indicated significant diversity from published SRLV sequences. Deletions and hypervariability in the 5' part of the env gene have implications for the size of the proposed CAEV-1GA Rev protein and the encoded surface glycoprotein (SU). The variable regions in the C-terminal part of SU obtained from Norwegian CAEV isolates demonstrate higher sequence divergence than has been described previously for SRLVs. Phylogenetic analysis based on SU sequences gives further support for a unique group designation. The results described here reveal a distant genetic relationship between Norwegian CAEV and other SRLVs and demonstrate that there is more geographical heterogeneity among SRLVs than reported previously.

Amino Acid Sequence↗

Comparative DNA sequence analysis of wheat and rice genomes.

The use of DNA sequence-based comparative genomics for evolutionary studies and for transferring information from model species to crop species has revolutionized molecular genetics and crop improvement strategies. This study compared 4485 expressed sequence tags (ESTs) that were physically mapped in wheat chromosome bins, to the public rice genome sequence data from 2251 ordered BAC/PAC clones using BLAST. A rice genome view of homologous wheat genome locations based on comparative sequence analysis revealed numerous chromosomal rearrangements that will significantly complicate the use of rice as a model for cross-species transfer of information in nonconserved regions.

Chromosome Mapping↗

Identification of glycosylphosphatidylinositol-anchored proteins in Arabidopsis. A proteomic and genomic analysis.

In a recent bioinformatic analysis, we predicted the presence of multiple families of cell surface glycosylphosphatidylinositol (GPI)-anchored proteins (GAPs) in Arabidopsis (G.H.H. Borner, D.J. Sherrier, T.J. Stevens, I.T. Arkin, P. Dupree [2002] Plant Physiol 129: 486-499). A number of publications have since demonstrated the importance of predicted GAPs in diverse physiological processes including root development, cell wall integrity, and adhesion. However, direct experimental evidence for their GPI anchoring is mostly lacking. Here, we present the first, to our knowledge, large-scale proteomic identification of plant GAPs. Triton X-114 phase partitioning and sensitivity to phosphatidylinositol-specific phospholipase C were used to prepare GAP-rich fractions from Arabidopsis callus cells. Two-dimensional fluorescence difference gel electrophoresis and one-dimensional sodium dodecyl sulfate-polyacrylamide gel electrophoresis demonstrated the existence of a large number of phospholipase C-sensitive Arabidopsis proteins. Using liquid chromatography-tandem mass spectrometry, 30 GAPs were identified, including six beta-1,3 glucanases, five phytocyanins, four fasciclin-like arabinogalactan proteins, four receptor-like proteins, two Hedgehog-interacting-like proteins, two putative glycerophosphodiesterases, a lipid transfer-like protein, a COBRA-like protein, SKU5, and SKS1. These results validate our previous bioinformatic analysis of the Arabidopsis protein database. Using the confirmed GAPs from the proteomic analysis to train the search algorithm, as well as improved genomic annotation, an updated in silico screen yielded 64 new candidates, raising the total to 248 predicted GAPs in Arabidopsis.

Arabidopsis↗

Assignment of orthologous genes via genome rearrangement.

The assignment of orthologous genes between a pair of genomes is a fundamental and challenging problem in comparative genomics. Existing methods that assign orthologs based on the similarity between DNA or protein sequences may make erroneous assignments when sequence similarity does not clearly delineate the evolutionary relationship among genes of the same families. In this paper, we present a new approach to ortholog assignment that takes into account both sequence similarity and evolutionary events at a genome level, where orthologous genes are assumed to correspond to each other in the most parsimonious evolving scenario under genome rearrangement. First, the problem is formulated as that of computing the signed reversal distance with duplicates between the two genomes of interest. Then, the problem is decomposed into two new optimization problems, called minimum common partition and maximum cycle decomposition, for which efficient heuristic algorithms are given. Following this approach, we have implemented a high-throughput system for assigning orthologs on a genome scale, called SOAR, and tested it on both simulated data and real genome sequence data. Compared to a recent ortholog assignment method based entirely on homology search (called INPARANOID), SOAR shows a marginally better performance in terms of sensitivity on the real data set because it is able to identify several correct orthologous pairs that are missed by INPARANOID. The simulation results demonstrate that SOAR, in general, performs better than the iterated exemplar algorithm in terms of computing the reversal distance and assigning correct orthologs.

Algorithms↗

ESTMAP: a system for expressed sequence tags mapping on genomic sequences.

The completion of a number of large genome sequencing projects emphasizes the importance of protein-coding gene predictions. Most of the problems associated with gene prediction are caused by the complex exon-intron structures commonly found in eukaryotic genomes. However, information from homologous sequences can significantly improve the accuracy of the prediction. In particular, expressed sequence tags (ESTs) are very useful for this purpose, since currently existing EST collections are very large. We developed an ESTMAP system, which utilizes homology searches against a database of repetitive elements using the RepeatView program and the EST Division of GenBank using the BLASTN program. ESTMAP extracts "exact" matches with EST sequences (> 95% of homology) from BLASTN output file and predicts introns in DNA comparing ESTs and a query sequence. ESTMAP is implemented as a part of the WebGene system (http://www.cnr.it/webgene).

Base Sequence↗

Use of the Rhodopseudomonas palustris genome sequence to identify a single amino acid that contributes to the activity of a coenzyme A ligase with chlorinated substrates.

Rhodopseudomonas palustris strain RCB100 degrades 3-chlorobenzoate (3-CBA) anaerobically. We purified from this strain a coenzyme A ligase that is active with 3-CBA and determined its N-terminal amino acid sequence to be identical to that of a cyclohexanecarboxylate-CoA ligase encoded by aliA from the R. palustris strain (CGA009) that has been sequenced. Strain CGA009 differs from strain RCB100 in that it does not use 3-CBA as a sole carbon source. The aliA gene from the 3-CBA degrading strain differed by a single nucleotide from the aliA gene from strain CGA009, causing the substitution of a serine for a threonine at position 208. Both AliA enzymes, purified as His-tagged fusion proteins, had comparable activities with cyclohexanecarboxylate. However, AliA from the 3-CBA degrading strain was 10-fold more active with 3-CBA (kcat/Km of 4.3 x 10(4) M(-1) s(-1)) than the enzyme from the sequenced strain (kcat/Km 0.32 x 10(4) M(-1) s(-1)). The CGA009 enzyme was not sufficiently active with 3-CBA to complement an RCB100 aliA mutant for growth on this compound. Here, whole genome sequence information enabled us to identify a single nucleotide among 5.4 million nucleotides that contributes to the substrate preference of a coenzyme A ligase.

Acyl Coenzyme A↗

The alkyl tert-butyl ether intermediate 2-hydroxyisobutyrate is degraded via a novel cobalamin-dependent mutase pathway.

Fuel oxygenates such as methyl and ethyl tert-butyl ether (MTBE and ETBE, respectively) are degraded only by a limited number of bacterial strains. The aerobic pathway is generally thought to run via tert-butyl alcohol (TBA) and 2-hydroxyisobutyrate (2-HIBA), whereas further steps are unclear. We have now demonstrated for the newly isolated beta-proteobacterial strains L108 and L10, as well as for the closely related strain CIP I-2052, that 2-HIBA was degraded by a cobalamin-dependent enzymatic step. In these strains, growth on substrates containing the tert-butyl moiety, such as MTBE, TBA, and 2-HIBA, was strictly dependent on cobalt, which could be replaced by cobalamin. Tandem mass spectrometry identified a 2-HIBA-induced protein with high similarity to a peptide whose gene sequence was found in the finished genome of the MTBE-degrading strain Methylibium petroleiphilum PM1. Alignment analysis identified it as the small subunit of isobutyryl-coenzyme A (CoA) mutase (ICM; EC 5.4.99.13), which is a cobalamin-containing carbon skeleton-rearranging enzyme, originally described only in Streptomyces spp. Sequencing of the genes of both ICM subunits from strain L108 revealed nearly 100% identity with the corresponding peptide sequences from M. petroleiphilum PM1, suggesting a horizontal gene transfer event to have occurred between these strains. Enzyme activity was demonstrated in crude extracts of induced cells of strains L108 and L10, transforming 2-HIBA into 3-hydroxybutyrate in the presence of CoA and ATP. The physiological and evolutionary aspects of this novel pathway involved in MTBE and ETBE metabolism are discussed.

Amino Acid Sequence↗

Citrate synthase mutants of Agrobacterium are attenuated in virulence and display reduced vir gene induction.

A citrate synthase (CS) deletion mutant of Agrobacterium tumefaciens C58 is highly attenuated in virulence. The identity of the mutant was initially determined from its amino acid sequence, which is 68% identical to Escherichia coli and 77% identical to Brucella melitensis. The mutant lost all CS enzymatic activity, and a cloned CS gene complemented a CS mutation in Sinorhizobium. The CS mutation resulted in a 10-fold reduction in vir gene expression, which likely accounts for the attenuated virulence. When a plasmid containing a constitutive virG [virG(Con)] locus was introduced into this mutant, the level of vir gene induction was restored to nearly wild-type level. Further, the virG(Con)-complemented CS mutant strain induced tumors that were similar in size and number to those induced by the parental strain. The CS mutation resulted in only a minor reduction in growth rate in a glucose-salts medium. Both the CS mutant and the virG(Con)-complemented CS strain displayed similar growth deficiencies in a glucose-salts medium, indicating that the reduced growth rate of the CS mutant could not be responsible for the attenuated virulence. A search of the genome of A. tumefaciens C58 revealed four proteins, encoded on different replicons, with conserved CS motifs. However, only the locus that when mutated resulted in an attenuated phenotype has CS activity. Mutations in the other three loci did not result in attenuated virulence and any loss of CS activity, and none were able to complement the CS mutation in Sinorhizobium. The function of these loci remains unknown.

Agrobacterium tumefaciens↗

Electron microscopic analysis of partially replicated bacteriophage T7 DNA.

Partially replicated bacteriophage T7 DNA was isolated from Escherichia coli infected with UV-irradiated T7 bacteriophage and was analyzed by electron microscopy. The analysis determined the distribution of eye forms and forks in the partially replicated molecules. Eye forms and forks in unit length molecules were aligned with respect to the left end of the T7 genome, and segments were scored for replication in each molecule. The resulting histogram showed that only the left 25 to 30% of the molecules was replicated. Several different origins of DNA replication were used to initiate replication in the UV-irradiated experiments in which 32P-labeled progeny DNA from UV-irradiated phage was annealed with ordered restriction fragments of T7 DNA (K. B. Burck and R. C. Miller, Jr., Proc. Natl. Acad. Sci. U.S.A. 75:6144--6148, 1978). Both analyses support partial-replica hypotheses (N. A. Barricelli and A. H. Doermann, Virology 13:460--476, 1961; Doermann et al., J. Cell. comp. Physiol. 45[Suppl.]:51--74, 1955) as an explanation for the distribution of marker rescue frequencies during cross-reactivation; i.e., replication proceeds in a bidirectional manner from an origin to a site of UV damage, and those regions of the genome which replicate most efficiently are rescued most efficiently by a coinfecting phage. In addition, photoreactivation studies support the hypothesis that thymine dimers are the major UV damage blocking cross-reactivation in the right end of the T7 genome.

DNA Replication↗

Heterogeneity of the fibre sequence in subgenus C adenoviruses.

AIMS: To determine the nucleotide sequences of adenovirus (Ad) types 1 and 6 fibre genes; to clarify the molecular basis of the distinct haemagglutination properties of subgenus C Ads and their phylogenetic relations. METHODS: Human Ad1 and Ad6 fibre genes were sequenced from genomic DNA by direct sequencing. Primer selection was based on alignment of the fibre gene of human Ad serotypes Ad2 and Ad5. Fibre based subgenus C specific polymerase chain reaction (PCR) was performed to check for deletions in field isolates of Ad6, as revealed by sequence analysis of the Ad6 prototype. A phylogenetic tree was constructed from the predicted amino acid (AA) sequences of the fibre gene of important Ads. RESULTS: Ad1 and Ad6 comprise 1746 and 1584 nucleotides, encoding 582 and 528 AA, respectively. Ad6 showed deletions in motifs 15-17 (51 AA) of the shaft when compared with Ad1, Ad2, and Ad5. Subgenus C specific PCR with both prototype and field isolates also showed deletions in Ad6. In the shaft and knob, AA homology was 58.82-72.91% and 68.89-74.59%, respectively. The tail was 100% conserved. Phylogenetically, Ad1 and Ad6, including Ad2 and Ad5, formed a subgenus specific cluster, like other serotypes. CONCLUSIONS: The fibre gene (including the knob region) of subgenus C Ads is heterogeneous, providing the molecular basis for lack of crossreactivity in the haemagglutination inhibition test. This heterogeneity could be helpful in fibre based genotyping of subgenus C field isolates. Phylogeny might be useful for subgenus specific identification of important field strains.

Adenoviruses, Human↗

Species relationships between antifungal chitinase and nuclear rDNA (internal transcribed spacer) sequences in the genus Hordeum.

The sequences of the chitinase gene (Chi-26) and the internal transcribed spacer of 18S - 5.8S - 26S rDNA (ITS1) were determined to analyze the phylogenetic relationships among species representing the four basic genomes of the genus Hordeum. Grouping analysis based on data for Chi-26 gene sequences placed Hordeum secalinum (H genome) near the Hordeum murinum complex (Xu genome), and Hordeum bulbosum distant from the other species that carried the I genome. ITS sequence data showed the expected grouping based on the genome classification of the species studied. Different sequences of ITS were detected even in the genomes of the diploid species. The results are interpreted in terms of defective or unfinished concerted evolution processes in each taxon.

Antifungal Agents↗

Exploring alternative transcript structure in the human genome using blocks and InterPro.

Understanding how alternative splicing affects gene function is an important challenge facing modern-day molecular biology. Using homology-based, protein sequence analysis methods, it should be possible to investigate how transcript diversity impacts protein function. To test this, high-quality exon-intron structures were deduced for over 8000 human genes, including over 1300 (17 percent) that produce multiple transcript variants. A data mining technique (DiffMotif) was developed to identify genes in which transcript variation coincides with changes in conserved motifs between variants. Applying this method, we found that 30 percent of the multi-variant genes in our test set exhibited a differential profile of conserved InterPro and/or BLOCKS motifs across different mRNA variants. To investigate these, a visualization tool (ProtAnnot) that displays amino acid motifs in the context of genomic sequence was developed. Using this tool, genes revealed by the DiffMotif method were analyzed, and when possible, hypotheses regarding the potential role of alternative transcript structure in modulating gene function were developed. Examples of these, including: MEOX1, a homeobox-containing protein; AIRE, involved in auto-immune disease; PLAT, tissue type plasminogen activator; and CD79b, a component of the B-cell receptor complex, are presented. These results demonstrate that amino acid motif databases like BLOCKS and InterPro are useful tools for investigating how alternative transcript structure affects gene function.

Algorithms↗