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

Structural similarity to link sequence space: new potential superfamilies and implications for structural genomics.

The current pace of structural biology now means that protein three-dimensional structure can be known before protein function, making methods for assigning homology via structure comparison of growing importance. Previous research has suggested that sequence similarity after structure-based alignment is one of the best discriminators of homology and often functional similarity. Here, we exploit this observation, together with a merger of protein structure and sequence databases, to predict distant homologous relationships. We use the Structural Classification of Proteins (SCOP) database to link sequence alignments from the SMART and Pfam databases. We thus provide new alignments that could not be constructed easily in the absence of known three-dimensional structures. We then extend the method of Murzin (1993b) to assign statistical significance to sequence identities found after structural alignment and thus suggest the best link between diverse sequence families. We find that several distantly related protein sequence families can be linked with confidence, showing the approach to be a means for inferring homologous relationships and thus possible functions when proteins are of known structure but of unknown function. The analysis also finds several new potential superfamilies, where inspection of the associated alignments and superimpositions reveals conservation of unusual structural features or co-location of conserved amino acids and bound substrates. We discuss implications for Structural Genomics initiatives and for improvements to sequence comparison methods.

Amino Acid Sequence↗
PubMed1995

Adelaide River virus nucleoprotein gene: analysis of phylogenetic relationships of ephemeroviruses and other rhabdoviruses.

The nucleotide sequence of the Adelaide River virus (ARV) genome was determined from the 3' terminus to the end of the nucleoprotein (N) gene. The 3' leader sequence comprises 50 nucleotides and shares a common terminal trinucleotide (3' UGC-), a conserved U-rich domain and a variable AU-rich domain with other animal rhabdoviruses. The N gene comprises 1355 nucleotides from the transcription start sequence (AACAGG) to the poly(A) sequence [CATG(A)7] and encodes a polypeptide of 429 amino acids. The N protein has a calculated molecular mass of 49429 Da and a pI of 5.4 and, like the bovine ephemeral fever virus (BEFV) N protein, features a highly acidic C-terminal domain. Analysis of amino acid sequence relationships between all available rhabdovirus N proteins indicated that ARV and BEFV are closely related viruses (48.3% similarity) which share higher sequence similarity to vesiculoviruses than to lyssaviruses. Phylogenetic trees based on a multiple sequence alignment of all available rhabdovirus N protein sequences demonstrated clustering of viruses according to genome organization, host range and established taxonomic relationships.

Amino Acid Sequence↗
PubMed1998

Taxonomic characteristics of fijiviruses based on nucleotide sequences of the oat sterile dwarf virus genome.

Sequence determination of full-length cDNA clones of genomic segments 7-10 (S7-S10) of oat sterile dwarf fijivirus (OSDV) revealed that the 5' and 3' ends of the plus strands of these segments had the same conserved terminal sequences, 5' AACGAAAAA and UUUUUUUAGUC 3'. These sequences are similar, but not identical, to the conserved terminal nucleotide sequences of the genomic segments of rice black streaked dwarf fijivirus (RBSDV) and maize rough dwarf fijivirus (MRDV). The coding strands of S7 and S10 each contained two large nonoverlapping open reading frames (ORFs), as do RBSDV S7 and S9, MRDV S6 and S8 and Nilaparvata lugens reovirus (NLRV; a putative member of Fijivirus) S9. These results strongly suggest that the dicistronic nature of certain genomic segments is characteristic of fijiviruses. Computer analyses revealed sequence homology between RBSDV S7 ORF2, MRDV S6 ORF2 and OSDV S7 ORF2, suggesting that this protein is conserved among plant fijiviruses. No counterparts were found in the genome of NLRV, which is a nonphytopathogenic insect reovirus. Furthermore, phylogenetic trees derived from multiple sequence alignments of each of the homologous proteins from OSDV, RBSDV, MRDV and NLRV suggest that NLRV did not evolve from either Fijivirus group 2 (RBSDV and MRDV) or group 3 (OSDV).

Amino Acid Sequence↗
PubMed1995

Frequency and distribution of DNA uptake signal sequences in the Haemophilus influenzae Rd genome.

The naturally transformable, Gram-negative bacterium Haemophilus influenzae Rd preferentially takes up DNA of its own species by recognizing a 9-base pair sequence, 5'-AAGTGCGGT, carried in multiple copies in its chromosome. With the availability of the complete genome sequence, 1465 copies of the 9-base pair uptake site have been identified. Alignment of these sites unexpectedly reveals an extended consensus region of 29 base pairs containing the core 9-base pair region and two downstream 6-base pair A/T-rich regions, each spaced about one helix turn apart. Seventeen percent of the sites are in inverted repeat pairs, many of which are located downstream to gene termini and are capable of forming stem-loop structures in messenger RNA that might function as signals for transcription termination.

Base Composition↗
PubMed2003

Melanin biosynthesis in the fungus Curvularia lunata (teleomorph: Cochliobolus lunatus).

Curvularia lunata (teleomorph: Cochliobolus lunatus) is a known plant and human pathogen. Tricyclazole, a specific inhibitor of pentaketide melanin biosynthesis, blocked the biosynthesis of melanin in Curvularia lunata and caused the accumulation of the melanin metabolites flaviolin and 2-hydroxyjuglone. This showed that melanin in Curvularia lunata is produced by a pentaketide pathway from 1,8-dihydroxynaphthalene. The 1,3,8-trihydroxynaphthalene reductase (3HNR) gene, associated with the melanin pathway of Curvularia lunata, was identified and characterized. An alignment of 3HNR sequences enabled the design of primers covering conserved regions. A PCR-amplified fragment of Curvularia lunata genomic DNA was used for screening the cDNA library. Three independent cDNA clones revealed an 801-bp open reading frame encoding a 267 amino acid protein. The protein was expressed in Escherichia coli and purified to homogeneity. The predicted amino acid sequence of the 28.6-kDa protein demonstrated homology to other fungal 3HNR and other members of the short-chain dehydrogenase super family. Northern analyses revealed that 3HNR from Curvularia lunata is expressed synchronously with melanization after 3 days of Curvularia lunata growth in malt extract medium. No 3HNR reductase gene expression nor melanization was observed when Curvularia lunata was grown in yeast nitrogen base medium.

Ascomycota↗
PubMed2026

Molecular identification and diversity assessment of Tyrrhenian Romulea species (Iridaceae).

Taxonomic assignments based only on morphology are often insufficient for delimiting species, particularly in complexes shaped by hybridization and polyploidy, where species boundaries are unclear. This limitation hinders progress in ecological, biogeographic and conservation research. The genus Romulea, distributed across Africa and the Mediterranean Basin, exemplifies this challenge. Despite its remarkable diversity, Mediterranean Romulea has not received much attention from genetic and molecular studies. Here, we present the first multilocus genotype analysis of Mediterranean Romulea taxa, focusing on the Tyrrhenian biogeographic province. Using target-capture sequencing with the universal Angiosperms353 kit, we generated genomic data for 272 individuals representing 18 putative taxa. Our findings reveal genetic groups that align with current taxonomy, the existence of cryptic divergence, and highlight the role of hybridization. Furthermore, analysis of intra-individual genetic diversity suggests one or several allopolyploid origins for Mediterranean Romulea. Four taxa (R. assumptionis, R. revelieri, R. ligustica, R. rollii) are consistently well differentiated across nuclear and plastid datasets, supporting their recognition as distinct species. In contrast, the widespread species R. ramiflora and R. columnae contain well-differentiated groups that may represent cryptic speciation. Several other taxa, including R. x melitensis, R. corsica, and R. bulbocodium, exhibit genomic signatures consistent with hybrid origins. Plastid and nuclear variation patterns are consistent with a hypothesis of rapid radiation in the Tyrrhenian region. These results provide a primary genomic framework for the integrative taxonomy of Romulea.

Genetic Variation↗
PubMed2003

Whole-genome shotgun optical mapping of Rhodobacter sphaeroides strain 2.4.1 and its use for whole-genome shotgun sequence assembly.

Rhodobacter sphaeroides 2.4.1 is a facultative photoheterotrophic bacterium with tremendous metabolic diversity, which has significantly contributed to our understanding of the molecular genetics of photosynthesis, photoheterotrophy, nitrogen fixation, hydrogen metabolism, carbon dioxide fixation, taxis, and tetrapyrrole biosynthesis. To further understand this remarkable bacterium, and to accelerate an ongoing sequencing project, two whole-genome restriction maps (EcoRI and HindIII) of R. sphaeroides strain 2.4.1 were constructed using shotgun optical mapping. The approach directly mapped genomic DNA by the random mapping of single molecules. The two maps were used to facilitate sequence assembly by providing an optical scaffold for high-resolution alignment and verification of sequence contigs. Our results show that such maps facilitated the closure of sequence gaps by the early detection of nascent sequence contigs during the course of the whole-genome shotgun sequencing process.

Chromosomes, Bacterial↗
PubMed2004

Novel and deviant Walker A ATP-binding motifs in bacteriophage large terminase-DNA packaging proteins.

Bacteriophage terminases constitute a very interesting class of viral-coded multifunctional ATPase "motors" that apparently drive directional translocation of DNA into an empty viral capsid. A common Walker A motif and other conserved signatures of a critical ATPase catalytic center are identified in the N-terminal half of numerous large terminase proteins. However, several terminases, including the well-characterized lambda and SPP1 terminases, seem to lack the classic Walker A in the N-terminus. Using sequence alignment approaches, we discovered the presence of deviant Walker A motifs in these and many other phage terminases. One deviation, the presence of a lysine at the beginning of P-loop, may represent a 3D equivalent of the universally conserved lysine in the Walker A GKT/S signature. This and other novel putative Walker A motifs that first came to light through this study help define the ATPase centers of phage and viral terminases as well as elicit important insights into the molecular functioning of this fundamental motif in biological systems.

Adenosine Triphosphate↗
PubMed2004

MSARI: multiple sequence alignments for statistical detection of RNA secondary structure.

We present a highly accurate method for identifying genes with conserved RNA secondary structure by searching multiple sequence alignments of a large set of candidate orthologs for correlated arrangements of reverse-complementary regions. This approach is growing increasingly feasible as the genomes of ever more organisms are sequenced. A program called msari implements this method and is significantly more accurate than existing methods in the context of automatically generated alignments, making it particularly applicable to high-throughput scans. In our tests, it discerned clustalw-generated multiple sequence alignments of signal recognition particle or RNaseP orthologs from controls with 89.1% sensitivity at 97.5% specificity and with 74.4% sensitivity with no false positives in 494 controls. We used msari to conduct a comprehensive scan for secondary structure in mRNAs of coding genes, and we found many genes with known mRNA secondary structure and compelling evidence for secondary structure in other genes. msari uses a method for coping with sequence redundancy that is likely to have applications in a large set of other comparison-based search methods. The program is available for download from http://theory.csail.mit.edu/MSARi.

Algorithms↗
PubMedDate not supplied

StructSorter: a method for continuously updating a comprehensive protein structure alignment database.

Advances in protein crystallography and homology modeling techniques are producing vast amounts of high resolution protein structure data at ever increasing rates. As such, the ability to quickly and easily extract structural similarities is a key tool in discovering important functional relationships. We report on an approach for creating and maintaining a database of pairwise structure alignments for a comprehensive database comprising the PDB and homology models for the human and select pathogen genomes. Our approach consists of a novel, multistage method for determining pairwise structural similarity coupled with an efficient clustering protocol that approximates a full NxN assessment in a fraction of the time. Since biologists are commonly interested in recently released structures, and the homology models built from them, an automatically updating database of structural alignments has great value. Our approach yields a querying system that allows scientists to retrieve databank-wide protein structure similarities as easily as retrieving protein sequence similarities via BLAST or PSI-BLAST. Basic, noncommercial access to the database can be requested at https://tip.eidogen-sertanty.com/.

Databases, Protein↗
PubMed2006

A core linkage map of the bumblebee Bombus terrestris.

The bumblebee Bombus terrestris is an economically important pollinator and an emerging model species in quantitative and population genetics. We generated genetic linkage maps for 3 independent mapping populations of B. terrestris. The linkage map with the highest resolution had 21 linkage groups, which adequately represents the haploid chromosome number of B. terrestris (n = 18). This map can be considered saturated, with an average marker distance of 10.3 cM and an estimated genome coverage of 81%. Using flow cytometry, we have estimated the genome size of this species to be 625 Mb. With an estimated total recombination genome length of 2760 cM, this results in a ratio of 226 kb/cM between the physical and genetic genome sizes. A recurring set of microsatellites and amplified fragment length polymorphism (AFLP) markers allowed the alignment of 14 linkage groups between the 3 maps. We propose to adopt this core map as a reference tool for future genetic and molecular work in B. terrestris.

Animals↗
PubMed1999

Cloning and characterization of the murine genes for bHLH-ZIP transcription factors TFEC and TFEB reveal a common gene organization for all MiT subfamily members.

The microphthalmia-TFE (MiT) subfamily of basic helix-loop-helix leucine zipper (bHLH-ZIP) transcription factors, including TFE3, TFEB, TFEC, and Mitf, has been implicated in the regulation of tissue-specific gene expression in several cell lineages. In this report, we investigate the genomic organization and structural relatedness of MiT transcription factors. We characterized the gene for mTFEC, which covers a region of more than 50 kb and is composed of seven exons. Further, we cloned a cDNA for the murine TFEB homologue and characterized its genomic structure. The eight coding exons of mTFEB are distributed over a 6-kb region. A multiple alignment of amino acid sequences of known MiT subfamily members indicates undescribed, conserved N-terminal regions and common putative phosphorylation sites for TFE3, TFEB, and Mitf. Also, intron-exon borders for characterized MiT genes appear completely conserved. A new family member and closely related putative transcription factor in Caenorhabditis elegans was identified by database searches that show a similar genomic organization within the bHLH-ZIP region and the acidic domain. Evolutionary aspects and implications for structure-function relationships are discussed.

Amino Acid Sequence↗
PubMed1997

A tool for aligning very similar DNA sequences.

RESULTS: We have produced a computer program, named sim3, that solves the following computational problem. Two DNA sequences are given, where the shorter sequence is very similar to some contiguous region of the longer sequence. Sim3 determines such a similar region of the longer sequence, and then computes an optimal set of single-nucleotide changes (i.e. insertions, deletions or substitutions) that will convert the shorter sequence to that region. Thus, the alignment scoring scheme is designed to model sequencing errors, rather than evolutionary processes. The program can align a 100 kb sequence to a 1 megabase sequence in a few seconds on a workstation, provided that there are very few differences between the shorter sequence and some region in the longer sequence. The program has been used to assemble sequence data for the Genomes Division at the National Center for Biotechnology Information.

Aldehyde Dehydrogenase↗
PubMed2000

The 5S rRNA gene diversity in Kengyilia rigidula (Keng and S.L. Chen) J.L. Yang, Yen, and Baum (Poaceae: Triticeae): possible contribution of the H genome to the origin of Kengyilia.

Fifty-three units of 5S rDNA sequences from five accessions of Kengyilia rigidula, a member of the tribe Triticeae that also includes wheat, barley, rye, and their wild relatives, have been amplified by the polymerase chain reaction (PCR), cloned, and sequenced. The genome of K. rigidula consists of three haplomes, St, P, and Y. An evaluation of the aligned sequences of the diverse 53 different 5S DNA units yielded three 5S-unit classes. One unit class, Long S1, was assignable to the St haplome, one unit class, the Long P1, was assignable to the P haplome, and a third unit class, Long H1, was assignable to the H haplome. The last was expected to be assignable to the Y haplome, based on previous knowledge. Evolutionary scenarios are put forward to explain this finding. Among those possibilities is that the number of copies of units assignable to the Y haplome is very small and difficult to detect. Short units, reported earlier in K. alatavica, were not found in K. rigidula.

Base Sequence↗
PubMed1997

Analysis of the structure and expression of the chicken gene encoding a homolog of the human RREB-1 transcription factor.

Ras proteins are involved in a number of signal transduction pathways including the mitogen-activated kinase cascade. Activated MAPKs translocate to the nucleus and phosphorylate transcription factors such as c-myc, TCF and AP-1. Recently, a Ras-responsive element binding transcription factor, RREB-1, was cloned from a human medullary thyroid carcinoma cell line. RREB-1 is a zinc finger protein that binds to a Ras-responsive element in the promoter of the human calcitonin gene. We report the cloning of the chicken homologue to human RREB-1. Amino-acid alignment demonstrates that chicken and human RREB-1 are 53% identical and 69% similar. Genomic southern analysis indicates that chicken rreb-1 is a single-copy gene in the chicken genome. We demonstrate that chicken and human rreb-1 display the same tissue distribution, being expressed in all tissues examined except the brain. Interestingly, chicken RREB-1 has an extended N-terminus and contains 16 zinc fingers of the TFIIIA subclass, in comparison to human RREB-1 which was reported to contain only four zinc fingers. The size discrepancy between the two predicted gene products is further discussed. An unusual structural feature of RREB-1 is the widely spaced arrangement of the zinc fingers.

Amino Acid Sequence↗
PubMed2002

The human endogenous retrovirus family HERV-K(HML-3).

A substantial amount of the human genome is composed of human endogenous retroviruses (HERVs). Manifold HERV families have been identified, among them several so-called HERV-K(HML) families. Although the HERV-K(HML-2) family has been studied in detail, other HERV-K families are not as well characterized. We describe here the HERV-K HML-3 family in more detail. We estimate that there are about 140 proviral loci or remains of such per haploid genome. Most loci are severely mutated. Proviruses displaying larger deletions in gag and pol are common. A multiple alignment of 73 HERV-K(HML-3) sequences displays several potentially important differences compared with the HERVK9I sequence in Repbase. A consensus sequence with open reading frames for all retroviral genes was generated, for which intact dUTPase motifs and env gene variants with different coding capacities are observed. Phylogenetic analysis shows near-monophyly with distinction of two closely related subgroups. Proviruses formed about 36 million years ago. However, no continuous activity through primate evolution is indicated.

Base Sequence↗
PubMed1999

Benchmarking PSI-BLAST in genome annotation.

The recognition of remote protein homologies is a major aspect of the structural and functional annotation of newly determined genomes. Here we benchmark the coverage and error rate of genome annotation using the widely used homology-searching program PSI-BLAST (position-specific iterated basic local alignment search tool). This study evaluates the one-to-many success rate for recognition, as often there are several homologues in the database and only one needs to be identified for annotating the sequence. In contrast, previous benchmarks considered one-to-one recognition in which a single query was required to find a particular target. The benchmark constructs a model genome from the full sequences of the structural classification of protein (SCOP) database and searches against a target library of remote homologous domains (<20 % identity). The structural benchmark provides a reliable list of correct and false homology assignments. PSI-BLAST successfully annotated 40 % of the domains in the model genome that had at least one homologue in the target library. This coverage is more than three times that if one-to-one recognition is evaluated (11 % coverage of domains). Although a structural benchmark was used, the results equally apply to just sequence homology searches. Accordingly, structural and sequence assignments were made to the sequences of Mycoplasma genitalium and Mycobacterium tuberculosis (see http://www.bmm.icnet. uk). The extent of missed assignments and of new superfamilies can be estimated for these genomes for both structural and functional annotations.

Algorithms↗