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,279 records · Page 71Linked to original sources
PubMed2003

CORG: a database for COmparative Regulatory Genomics.

Sequence conservation in non-coding, upstream regions of orthologous genes from man and mouse is likely to reflect common regulatory DNA sites. Motivated by this assumption we have delineated a catalogue of conserved non-coding sequence blocks and provide the CORG-'COmparative Regulatory Genomics'-database. The data were computed based on statistically significant local suboptimal alignments of 15 kb regions upstream of the translation start sites of, currently, 10 793 pairs of orthologous genes. The resulting conserved non-coding blocks were annotated with EST matches for easier detection of non-coding mRNA and with hits to known transcription factor binding sites. CORG data are accessible from the ENSEMBL web site via a DAS service as well as a specially developed web service (http://corg.molgen.mpg.de) for query and interactive visualization of the conserved blocks and their annotation.

Animals↗
PubMed1985

Significance of similarities in patterns: an application to beta interferon-related DNA on human chromosome 2.

The nucleotide sequence of a 14-kilobase (kb) region of the human beta interferon (IFN-beta)-related DNA locus on chromosome 2 (genomic DNA clone lambda B3) was determined and compared to that of the IFN-beta 1 gene by using the Sellers TT algorithm. This algorithm aligns segments of one sequence with similar segments in a second sequence. A strategy was developed for assessing the significance of similarities between DNA sequences based on a scheme that recognizes patterns or runs of identities within an alignment. The pattern score (II) thus obtained is an entropy-like measure. Numerically it is a reflection of the length of the second longest run of identity in an alignment plus a correction factor due to the other shorter identity runs in the alignment. When the IFN-beta 1 gene is compared to a random nucleotide sequence, the distribution of II scores in such comparisons fits a Gaussian function. This strategy has been used to identify seven segments along one strand of lambda B3 DNA that are related to segments in IFN-beta 1; these seven alignments have II scores greater than or equal to 3 standard deviations above the mean score obtained in comparisons between IFN-beta 1 and random nucleotide sequences. One of these alignments (section 7) has a II score 8.02 standard deviations above this mean score. The likelihood of finding an alignment statement as good as that in section 7 in a random sequence the length of the human genome is approximately 10(-7). Furthermore, the lambda B3 DNA sequence in section 7 selects the human IFN-beta 1 gene as the most significant alignment in computer searches of mammalian nucleotide sequence data bases.

Base Sequence↗
PubMed2026

Strategies for mosaic variant calling in brain disorders.

The human brain is a genomic mosaic, where postzygotic mutations arising from embryogenesis to senescence drive diverse neurodevelopmental and neurodegenerative diseases. Because of numerous sequencing artifacts at ultralow variant allele frequencies (VAFs), detecting these variants remains a significant analytical challenge. This review focuses on single-nucleotide variants and small indels, summarizing current strategies for aligning sampling methods, including bulk, laser capture microdissection, and single-cell genomics, with the expected clonal architecture of the brain. It emphasizes that mosaic detection sensitivity is fundamentally constrained by sequencing depth, since even the most advanced algorithms cannot identify variants not physically represented in the sequencing library. The review further recommends the selection of variant calling algorithms based on validated VAF detection performance, matching tools like MuTect2 and MosaicForecast to their optimal performance ranges. Furthermore, we discuss how multitissue sampling, as emphasized by the SMaHT project, addresses the matched-control dilemma and supports accurate variant classification via cross-tissue VAF gradients. Integrating these established pipelines with multiomics modalities, including transcriptomic and epigenetic data, could advance the field toward a functional understanding of how the somatic genome impacts human brain health and disease.

Humans↗
PubMed2005

RNA secondary structural alignment with conditional random fields.

MOTIVATION: The computational identification of non-coding RNA regions on the genome is currently receiving much attention. However, it is essentially harder than gene-finding problems for protein-coding regions because non-coding RNA sequences do not have strong statistical signals. Since comparative sequence analysis is effective for non-coding RNA detection, efficient computational methods are expected for structural alignment of RNA sequences. Several methods have been proposed to accomplish the structural alignment tasks for RNA sequences, and we found that one of the most important points is to estimate an accurate score matrix for calculating structural alignments. RESULTS: We propose a novel approach for RNA structural alignment based on conditional random fields (CRFs). Our approach has some specific features compared with previous methods in the sense that the parameters for structural alignment are estimated such that the model can most probably discriminate between correct alignments and incorrect alignments, and has the generalization ability so that a satisfiable score matrix can be obtained even with a small number of sample data without overfitting. Experimental results clearly show that the parameter estimation with CRFs can outperform all the other existing methods for structural alignments of RNA sequences. Furthermore, structural alignment search based on CRFs is more accurate for predicting non-coding RNA regions than the other scoring methods. These experimental results strongly support our discriminative method employing CRFs to estimate the score matrix parameters. AVAILABILITY: The program which is implemented in C++ is available at http://phmmts.dna.bio.keio.ac.jp/ under the GNU public license.

Base Sequence↗
PubMed2006

Identification of Salmonella enterica serovar Typhimurium using specific PCR primers obtained by comparative genomics in Salmonella serovars.

Salmonella enterica serovar Typhimurium is a major foodborne pathogen throughout the world. Until now, the specific target genes for the detection and identification of serovar Typhimurium have not been developed. To determine the specific probes for serovar Typhimurium, the genes of serovar Typhimurium LT2 that were expected to be unique were selected with the BLAST (Basic Local Alignment Search Tool) program within GenBank. The selected genes were compared with 11 genomic sequences of various Salmonella serovars by BLAST. Of these selected genes, 10 were expected to be specific to serovar Typhimurium and were not related to virulence factor genes of Salmonella pathogenicity island or to genes of the O and H antigens of Salmonella. Primers for the 10 selected genes were constructed, and PCRs were evaluated with various genomic DNAs of Salmonella and non-Salmonella strains for the specific identification of Salmonella serovar Typhimurium. Among all the primer sets for the 10 genes, STM4497 showed the highest degree of specificity to serovar Typhimurium. In this study, a specific primer set for Salmonella serovar Typhimurium was developed on the basis of the comparison of genomic sequences between Salmonella serovars and was validated with PCR. This method of comparative genomics to select target genes or sequences can be applied to the specific detection of microorganisms.

Base Sequence↗
PubMed2003

A new sequence distance measure for phylogenetic tree construction.

MOTIVATION: Most existing approaches for phylogenetic inference use multiple alignment of sequences and assume some sort of an evolutionary model. The multiple alignment strategy does not work for all types of data, e.g. whole genome phylogeny, and the evolutionary models may not always be correct. We propose a new sequence distance measure based on the relative information between the sequences using Lempel-Ziv complexity. The distance matrix thus obtained can be used to construct phylogenetic trees. RESULTS: The proposed approach does not require sequence alignment and is totally automatic. The algorithm has successfully constructed consistent phylogenies for real and simulated data sets. AVAILABILITY: Available on request from the authors.

Algorithms↗
PubMed2006

CoSMoS: Conserved Sequence Motif Search in the proteome.

BACKGROUND: With the ever-increasing number of gene sequences in the public databases, generating and analyzing multiple sequence alignments becomes increasingly time consuming. Nevertheless it is a task performed on a regular basis by researchers in many labs. RESULTS: We have now created a database called CoSMoS to find the occurrences and at the same time evaluate the significance of sequence motifs and amino acids encoded in the whole genome of the model organism Escherichia coli K12. We provide a precomputed set of multiple sequence alignments for each individual E. coli protein with all of its homologues in the RefSeq database. The alignments themselves, information about the occurrence of sequence motifs together with information on the conservation of each of the more than 1.3 million amino acids encoded in the E. coli genome can be accessed via the web interface of CoSMoS. CONCLUSION: CoSMoS is a valuable tool to identify highly conserved sequence motifs, to find regions suitable for mutational studies in functional analyses and to predict important structural features in E. coli proteins.

Amino Acid Motifs↗
PubMed2003

Transcriptional control of the human sodium-coupled neutral amino acid transporter system A gene by amino acid availability is mediated by an intronic element.

System A amino acid transporter (SNAT2) gene expression is up-regulated at the transcriptional level in response to amino acid deprivation. Functional analysis of genomic fragments 5' upstream of the transcription start site, for both human and mouse SNAT2 genes showed that these regions exhibit promoter activity, but were amino acid unresponsive. However, when the human and mouse constructs were extended to include intron 1, it was observed that the rate of transcription was increased following amino acid deprivation. Deletion analysis of the human gene identified an intron 1 sequence spanning 54 nucleotides that was sufficient for conferring amino acid-dependent regulation to a minimal SNAT2 promoter. Alignment of the corresponding region from the human, mouse, and rat genomes revealed three highly conserved sequences. From site-directed mutagenesis, it was concluded that one of these sites functions as an amino acid response element (AARE) to regulate transcription. The core sequence of this site is identical to the AARE in the human CHOP gene. The SNAT2 AARE, along with a nearby conserved CAAT box, has enhancer activity in that it functions in an orientation and position independent manner, and it confers regulated transcription to a heterologous promoter.

Amino Acid Transport System A↗
PubMed2026

Climate-Driven Niche Tracking and Genomic Resilience Shape Future Distribution of a Widespread Agricultural Weed.

Understanding how agriculturally important species respond to environmental change is critical for maintaining productivity, mitigating agroecosystem threats and sustaining resilience. While crops have traditionally been the focus in agroecosystems, agricultural weeds are integral components that often face even stronger selective pressures, making them powerful models for investigating ecological and evolutionary responses to climatic and human-mediated challenges. Insights from how weeds adapt rapidly under these pressures can inform strategies to improve agricultural outcomes, since both pests and crops evolve under the same multivariate selective pressures. Here, we integrate two centuries of distribution records with whole-genome sequencing from natural populations of the most damaging weed in Europe-Alopecurus myosuroides (blackgrass) - to examine its ecological and evolutionary responses in agroecosystems. Blackgrass largely maintained its historical climatic niche, expanding its range primarily by tracking environments analogous to those it historically occupied. Genome-wide analyses revealed a polygenic basis of environmental responses, with most loci linked to single environmental variables and a subset showing limited environmental pleiotropy, indicating modular adaptation to the complex selective pressures of managed agricultural landscapes. Coupling these genomic-environment relationships with projected climate change and genomic offset analyses indicated that most blackgrass populations will remain well aligned with future conditions. Our findings show that ecological niche tracking and polygenic adaptation allow agricultural weeds like blackgrass to persist under rapid environmental change, offering insights relevant not only for weed management but also for designing resilient cropping systems under future climates.

Plant Weeds↗
PubMed2006

Integration of microsatellite-based genetic maps for the turkey (Meleagris gallopavo).

Integration of turkey genetic maps and their associated markers is essential to increase marker density in support of map-based genetic studies. The objectives of this study were to integrate 2 microsatellite-based turkey genetic maps--the Roslin map and the University of Minnesota (UMN) map--by genotyping markers from the Roslin study on the mapping families of the UMN study. A total of 279 markers was tested, and 240 were subsequently screened for polymorphisms in the UMN/Nicholas Turkey Breeding Farms (NTBF) mapping families. Of the 240 markers, 89 were genetically informative and were used for genotyping the F2 offspring. Significant genetic linkages (log of odds > 3.0) were found for 84 markers from the Roslin study. BLASTn comparison of marker sequences with the draft assembly of the chicken genome found 263 significant matches. The combination of genetic and in silico mapping allowed for the alignment of all linkage groups of the Roslin map with those of the UMN map. With the addition of the markers from the Roslin map, 438 markers are now genetically linked in the UMN/NTBF families, and more than 1700 turkey sequences have now been assigned to likely positions in the chicken-genome sequence.

Animals↗
PubMed2005

Distribution and intensity of constraint in mammalian genomic sequence.

Comparisons of orthologous genomic DNA sequences can be used to characterize regions that have been subject to purifying selection and are enriched for functional elements. We here present the results of such an analysis on an alignment of sequences from 29 mammalian species. The alignment captures approximately 3.9 neutral substitutions per site and spans approximately 1.9 Mbp of the human genome. We identify constrained elements from 3 bp to over 1 kbp in length, covering approximately 5.5% of the human locus. Our estimate for the total amount of nonexonic constraint experienced by this locus is roughly twice that for exonic constraint. Constrained elements tend to cluster, and we identify large constrained regions that correspond well with known functional elements. While constraint density inversely correlates with mobile element density, we also show the presence of unambiguously constrained elements overlapping mammalian ancestral repeats. In addition, we describe a number of elements in this region that have undergone intense purifying selection throughout mammalian evolution, and we show that these important elements are more numerous than previously thought. These results were obtained with Genomic Evolutionary Rate Profiling (GERP), a statistically rigorous and biologically transparent framework for constrained element identification. GERP identifies regions at high resolution that exhibit nucleotide substitution deficits, and measures these deficits as "rejected substitutions". Rejected substitutions reflect the intensity of past purifying selection and are used to rank and characterize constrained elements. We anticipate that GERP and the types of analyses it facilitates will provide further insights and improved annotation for the human genome as mammalian genome sequence data become richer.

Animals↗
PubMed2004

Characterization of pancreatic ERj3p, a homolog of yeast DnaJ-like protein Scj1p.

We have previously identified in the human EST sequence data base four overlapping clones that could be aligned with both a predicted protein sequence, deduced from the C. elegans genomic sequence, and partial amino acid sequences, obtained for a protein from canine pancreatic microsomes. We suggested that these proteins are homologs of yeast microsomal and DnaJ-like protein Scj1p and termed them ERj3p. Here we verified the predicted protein sequence of human ERj3p by sequence analysis of the corresponding cDNA. Multiple alignment of related sequences identified these proteins as true homologs of yeast Scj1p. Biochemical analysis of the canine protein characterized ERj3p as a soluble glycoprotein of the pancreatic endoplasmic reticulum. This pancreatic DnaJ-like protein was shown to interact with lumenal DnaK-like proteins, such as BiP. Furthermore, we found that ERj3p interacts with SDF2L1 protein that may be involved in protein O-glycosylation. We propose that ERj3p represents a cochaperone of DnaK-like chaperones of the mammalian endoplasmic reticulum and is involved in folding and maturation of newly synthesized proteins.

Amino Acid Sequence↗
PubMed2001

Genome trees constructed using five different approaches suggest new major bacterial clades.

BACKGROUND: The availability of multiple complete genome sequences from diverse taxa prompts the development of new phylogenetic approaches, which attempt to incorporate information derived from comparative analysis of complete gene sets or large subsets thereof. Such attempts are particularly relevant because of the major role of horizontal gene transfer and lineage-specific gene loss, at least in the evolution of prokaryotes. RESULTS: Five largely independent approaches were employed to construct trees for completely sequenced bacterial and archaeal genomes: i) presence-absence of genomes in clusters of orthologous genes; ii) conservation of local gene order (gene pairs) among prokaryotic genomes; iii) parameters of identity distribution for probable orthologs; iv) analysis of concatenated alignments of ribosomal proteins; v) comparison of trees constructed for multiple protein families. All constructed trees support the separation of the two primary prokaryotic domains, bacteria and archaea, as well as some terminal bifurcations within the bacterial and archaeal domains. Beyond these obvious groupings, the trees made with different methods appeared to differ substantially in terms of the relative contributions of phylogenetic relationships and similarities in gene repertoires caused by similar life styles and horizontal gene transfer to the tree topology. The trees based on presence-absence of genomes in orthologous clusters and the trees based on conserved gene pairs appear to be strongly affected by gene loss and horizontal gene transfer. The trees based on identity distributions for orthologs and particularly the tree made of concatenated ribosomal protein sequences seemed to carry a stronger phylogenetic signal. The latter tree supported three potential high-level bacterial clades,: i) Chlamydia-Spirochetes, ii) Thermotogales-Aquificales (bacterial hyperthermophiles), and ii) Actinomycetes-Deinococcales-Cyanobacteria. The latter group also appeared to join the low-GC Gram-positive bacteria at a deeper tree node. These new groupings of bacteria were supported by the analysis of alternative topologies in the concatenated ribosomal protein tree using the Kishino-Hasegawa test and by a census of the topologies of 132 individual groups of orthologous proteins. Additionally, the results of this analysis put into question the sister-group relationship between the two major archaeal groups, Euryarchaeota and Crenarchaeota, and suggest instead that Euryarchaeota might be a paraphyletic group with respect to Crenarchaeota. CONCLUSIONS: We conclude that, the extensive horizontal gene flow and lineage-specific gene loss notwithstanding, extension of phylogenetic analysis to the genome scale has the potential of uncovering deep evolutionary relationships between prokaryotic lineages.

Bacteria↗
PubMed2004

Organization of the human aromatase p450 (CYP19) gene.

The human CYP19 (p450arom) gene is located in the 21.2 region on the long arm of chromosome 15 (15q21.2). This gene spans a region that consists of a 30 kb coding region and a 93 kb regulatory region ( approximately 123 kb total length). Its regulatory region contains at least 10 distinct promoters regulated in a tissue- or signaling pathway-specific manner. The Human Genome Project data published in 2000 enabled us to accurately align these promoters within the 93 kb regulatory region of the p450arom gene. Each promoter is regulated by a distinct set of regulatory sequences in DNA and transcription factors that bind to these specific sequences. In most vertebrates, p450arom expression is under the control of gonad- and brain-specific promoters. In humans, however, there are at least eight additional promoters that were apparently recruited throughout evolution, possibly via alterations in DNA. A critical mechanism that permits the use of such a large number of promoters seems to be the extremely promiscuous nature of the common splice acceptor site because, with activation of each promoter, an untranslated first exon is spliced onto this common junction immediately upstream of the translation start site in the coding region. These partially tissue-specific promoters are used in the gonads, bone, brain, vascular tissue, adipose tissue, skin, fetal liver, and placenta for estrogen biosynthesis necessary for human physiology. Ovary and testis use promoter II, which is located immediately upstream of the coding region. The adipose tissue in general, including adipose tissue of the disease-free breast, on the other hand, maintains low levels of aromatase expression primarily via promoter I.4, which lies 73 kb upstream of the common coding region. Promoters I.3 and II are used only minimally in normal breast adipose tissue. Promoter II and I.3 activities in breast cancer tissue, however, are strikingly increased. Additionally, the endothelial-type promoter I.7 is also upregulated in breast cancer. Therefore, breast tumor tissue takes advantage of four promoters (II, I.3, I.7, and I.4) for aromatase expression and estrogen production. The sum of p450arom mRNA species arising from these four promoters contributes significantly to elevated levels of total p450arom mRNA in breast cancer in contrast to the normal breast that uses promoter I.4. Because each mRNA species contains the identical coding region regardless of the variable untranslated first exon, the encoded protein functions as the aromatase enzyme in each case.

Aromatase↗
PubMed2005

Comparative plant genomics resources at PlantGDB.

PlantGDB (http://www.plantgdb.org/) is a database of plant molecular sequences. Expressed sequence tag (EST) sequences are assembled into contigs that represent tentative unique genes. EST contigs are functionally annotated with information derived from known protein sequences that are highly similar to the putative translation products. Tentative Gene Ontology terms are assigned to match those of the similar sequences identified. Genome survey sequences are assembled similarly. The resulting genome survey sequence contigs are matched to ESTs and conserved protein homologs to identify putative full-length open reading frame-containing genes, which are subsequently provisionally classified according to established gene family designations. For Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa), the exon-intron boundaries for gene structures are annotated by spliced alignment of ESTs and full-length cDNAs to their respective complete genome sequences. Unique genome browsers have been developed to present all available EST and cDNA evidence for current transcript models (for Arabidopsis, see the AtGDB site at http://www.plantgdb.org/AtGDB/; for rice, see the OsGDB site at http://www.plantgdb.org/OsGDB/). In addition, a number of bioinformatic tools have been integrated at PlantGDB that enable researchers to carry out sequence analyses on-site using both their own data and data residing within the database.

Computational Biology↗
PubMed2006

Detecting overlapping coding sequences in virus genomes.

BACKGROUND: Detecting new coding sequences (CDSs) in viral genomes can be difficult for several reasons. The typically compact genomes often contain a number of overlapping coding and non-coding functional elements, which can result in unusual patterns of codon usage; conservation between related sequences can be difficult to interpret--especially within overlapping genes; and viruses often employ non-canonical translational mechanisms--e.g. frameshifting, stop codon read-through, leaky-scanning and internal ribosome entry sites--which can conceal potentially coding open reading frames (ORFs). RESULTS: In a previous paper we introduced a new statistic--MLOGD (Maximum Likelihood Overlapping Gene Detector)--for detecting and analysing overlapping CDSs. Here we present (a) an improved MLOGD statistic, (b) a greatly extended suite of software using MLOGD, (c) a database of results for 640 virus sequence alignments, and (d) a web-interface to the software and database. Tests show that, from an alignment with just 20 mutations, MLOGD can discriminate non-overlapping CDSs from non-coding ORFs with a typical accuracy of up to 98%, and can detect CDSs overlapping known CDSs with a typical accuracy of 90%. In addition, the software produces a variety of statistics and graphics, useful for analysing an input multiple sequence alignment. CONCLUSION: MLOGD is an easy-to-use tool for virus genome annotation, detecting new CDSs--in particular overlapping or short CDSs--and for analysing overlapping CDSs following frameshift sites. The software, web-server, database and supplementary material are available at http://guinevere.otago.ac.nz/mlogd.html.

Algorithms↗
PubMed1999

The Tangled Web: Gene Genealogies and the Origin of Eukaryotes.

Accessing data from the genomes of organisms (individual genes) and analyzing these data using sophisticated alignment and phylogenetic methods led to the expectation that we would be able to paint a clear picture of the evolution of eukaryotes. Previous analyses based on morphology and ultrastructure failed to pinpoint both the sister taxon to eukaryotes and the branching order of eukaryotic lineages. However, the expectation that molecular data would provide resolution has not been met since a growing number of gene genealogies present conflicting hypotheses for the origin and diversification of eukaryotes. Instead of reconstructing a simple bifurcating tree of life, these gene genealogies have generated a complex picture of eukaryotic genomes whereby ancient lateral transfers (of individual genes or perhaps even entire genomes) has tangled the evolutionary history of eukaryotes. Resolution of these conflicting genealogies comes in recognizing that eukaryotes are chimeric, containing genetic information from multiple ancestral lineages.

Journal Article↗
PubMed1992

Identification and characterization of foot-and-mouth disease virus O1 Burgwedel/1987 as an intertypic recombinant.

The foot-and-mouth disease virus field isolate Burgwedel/1987 subtype O1 was found to differ genetically from the antigenically related strain O1 Kaufbeuren within the region encoding the non-structural proteins. This genetic difference was indicated by the RNase mismatch cleavage method and confirmed by nucleotide sequencing. An alignment of sequences encoding proteinase 3C of the Burgwedel isolate and several other virus strains identified this isolate as an intertypic recombinant; the parent strains were O1 Kaufbeuren and a subtype C1 strain. Recombination occurred between nucleotide positions 5493 and 5521, within the region encoding peptide 3B1. Thus, the 5' three-quarters of the O1 genome were fused to the 3'-terminal quarter of the C1 genome. Other contemporary isolates from the same district are not recombinants. Sequence alignment distinguished four patterns of proteinase 3C-coding sequences among the virus strains analysed: subtypes A12, C1 and O1 exhibit one pattern each, and another pattern is common to subtypes A5, A10 and O2.

Aphthovirus↗