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,459 records · Page 81Linked to original sources
PubMed2014

Inter-genomic sexual conflict drives antagonistic coevolution in harvester ants.

The reproductive interests of males and females are not always aligned, leading to sexual conflict over parental investment, rate of reproduction and mate choice. Traits that increase the genetic interests of one sex often occur at the expense of the other, selecting for counter-adaptations leading to antagonistic coevolution. Reproductive conflict is not limited to intraspecific interactions; interspecific hybridization can produce pronounced sexual conflict between males and females of different species, but it is unclear whether such conflict can drive sexually antagonistic coevolution between reproductively isolated genomes. We tested for hybridization-driven sexually antagonistic adaptations in queens and males of the socially hybridogenetic 'J' lineages of Pogonomyrmex harvester ants, whose mating system promotes hybridization in queens but selects against it in males. We conducted no-choice mating assays to compare patterns of mating behaviour and sperm transfer between inter- and intra-lineage pairings. There was no evidence for mate discrimination on the basis of pair type, and the total quantity of sperm transferred did not differ between intra- and inter-lineage pairs; however, further dissection of the sperm transfer process into distinct mechanistic components revealed significant, and opposing, cryptic manipulation of copulatory investment by both sexes. Males of both lineages increased their rate of sperm transfer to high-fitness intra-lineage mates, with a stronger response in the rarer lineage for whom mating mistakes are the most likely. By contrast, the total duration of copulation for intra-lineage mating pairs was significantly shorter than for inter-lineage crosses, suggesting that queens respond to prevent excessive sperm loading by prematurely terminating copulation. These findings demonstrate that sexual conflict can lead to antagonistic coevolution in both intra-genomic and inter-genomic contexts. Indeed, the resolution of sexual conflict may be a key determinant of the long-term evolutionary potential of host-dependent reproductive strategies, counteracting the inherent instabilities arising from such systems.

Animals↗
PubMed2000

Comparative protein structure modeling of genes and genomes.

Comparative modeling predicts the three-dimensional structure of a given protein sequence (target) based primarily on its alignment to one or more proteins of known structure (templates). The prediction process consists of fold assignment, target-template alignment, model building, and model evaluation. The number of protein sequences that can be modeled and the accuracy of the predictions are increasing steadily because of the growth in the number of known protein structures and because of the improvements in the modeling software. Further advances are necessary in recognizing weak sequence-structure similarities, aligning sequences with structures, modeling of rigid body shifts, distortions, loops and side chains, as well as detecting errors in a model. Despite these problems, it is currently possible to model with useful accuracy significant parts of approximately one third of all known protein sequences. The use of individual comparative models in biology is already rewarding and increasingly widespread. A major new challenge for comparative modeling is the integration of it with the torrents of data from genome sequencing projects as well as from functional and structural genomics. In particular, there is a need to develop an automated, rapid, robust, sensitive, and accurate comparative modeling pipeline applicable to whole genomes. Such large-scale modeling is likely to encourage new kinds of applications for the many resulting models, based on their large number and completeness at the level of the family, organism, or functional network.

Animals↗
PubMed2006

Multiple mapping method: a novel approach to the sequence-to-structure alignment problem in comparative protein structure modeling.

A major bottleneck in comparative protein structure modeling is the quality of input alignment between the target sequence and the template structure. A number of alignment methods are available, but none of these techniques produce consistently good solutions for all cases. Alignments produced by alternative methods may be superior in certain segments but inferior in others when compared to each other; therefore, an accurate solution often requires an optimal combination of them. To address this problem, we have developed a new approach, Multiple Mapping Method (MMM). The algorithm first identifies the alternatively aligned regions from a set of input alignments. These alternatively aligned segments are scored using a composite scoring function, which determines their fitness within the structural environment of the template. The best scoring regions from a set of alternative segments are combined with the core part of the alignments to produce the final MMM alignment. The algorithm was tested on a dataset of 1400 protein pairs using 11 combinations of two to four alignment methods. In all cases MMM showed statistically significant improvement by reducing alignment errors in the range of 3 to 17%. MMM also compared favorably over two alignment meta-servers. The algorithm is computationally efficient; therefore, it is a suitable tool for genome scale modeling studies.

Algorithms↗
PubMed2000

Evaluation of PSI-BLAST alignment accuracy in comparison to structural alignments.

The PSI-BLAST algorithm has been acknowledged as one of the most powerful tools for detecting remote evolutionary relationships by sequence considerations only. This has been demonstrated by its ability to recognize remote structural homologues and by the greatest coverage it enables in annotation of a complete genome. Although recognizing the correct fold of a sequence is of major importance, the accuracy of the alignment is crucial for the success of modeling one sequence by the structure of its remote homologue. Here we assess the accuracy of PSI-BLAST alignments on a stringent database of 123 structurally similar, sequence-dissimilar pairs of proteins, by comparing them to the alignments defined on a structural basis. Each protein sequence is compared to a nonredundant database of the protein sequences by PSI-BLAST. Whenever a pair member detects its pair-mate, the positions that are aligned both in the sequential and structural alignments are determined, and the alignment sensitivity is expressed as the percentage of these positions out of the structural alignment. Fifty-two sequences detected their pair-mates (for 16 pairs the success was bi-directional when either pair member was used as a query). The average percentage of correctly aligned residues per structural alignment was 43.5+/-2.2%. Other properties of the alignments were also examined, such as the sensitivity vs. specificity and the change in these parameters over consecutive iterations. Notably, there is an improvement in alignment sensitivity over consecutive iterations, reaching an average of 50.9+/-2.5% within the five iterations tested in the current study.

Algorithms↗
PubMed2006

Involvement of some large immunophilins and their ligands in the protection and regeneration of neurons: a hypothetical mode of action.

The powerful immunosuppressive drugs such as FK506 and its derivatives induce some regeneration and protection of neurons from ischaemic brain injury and some other neurological disorders. The drugs form complexes with diverse FKBPs but apparently the FKBP52/FK506 complex was shown to be involved in the protection and regeneration of neurons. We used several different sequence attributes in searching diverse genomic databases for similar motifs as those present in the FKBPs. A Fortran library of algorithms (Par_Seq) has been designed and used in searching for the similarity of sequence motifs extracted from the multiple sequence alignments of diverse groups of proteins (query motifs) and the target motifs which are encoded in various genomes. The following sequence attributes were used in the establishment of the degree of convergence between: (A) amino acid (AA) sequence similarity (ID) of the query/target motifs and (B) their: (1) AA composition (AAC); (2) hydrophobicity (HI); (3) Jensen-Shannon entropy; and (4) AA propensity to form a particular secondary structure. The sequence hallmark of two different groups of peptidylprolyl cis/trans isomerases (PPIases), namely tetratricopetide repeat (TPR) motifs, which are present in the heat-shock cyclophilins and in the large FK506-binding proteins (FKBPs) were used to search various genomic databases. The Par_Seq algorithm has revealed that the TPR motifs have similar sequence attributes as a number of hydrophobic sequence segments of functionally unrelated membrane proteins, including some of the TMs from diverse G protein-coupled receptors (GPCRs). It is proposed that binding of the FKBP52/FK506 complex to the membranes via the TPR motifs and its interaction with some membrane proteins could be in part responsible for some neuro-regeneration and neuro-protection of the brain during some ischaemia-induced stresses.

Algorithms↗
PubMed2002

Targeting novel folds for structural genomics.

The ultimate goal of structural genomics is to obtain the structure of each protein coded by each gene within a genome to determine gene function. Because of cost and time limitations, it remains impractical to solve the structure for every gene product experimentally. Up to a point, reasonably accurate three-dimensional structures can be deduced for proteins with homologous sequences by using comparative modeling. Beyond this, fold recognition or threading methods can be used for proteins showing little homology to any known fold, although this is relatively time-consuming and limited by the library of template folds currently available. Therefore, it is appropriate to develop methods that can increase our knowledge base, expanding our fold libraries by earmarking potentially "novel" folds for experimental structure determination. How can we sift through proteomic data rapidly and yet reliably identify novel folds as targets for structural genomics? We have analyzed a number of simple methods that discriminate between "novel" and "known" folds. We propose that simple alignments of secondary structure elements using predicted secondary structure could potentially be a more selective method than both a simple fold recognition method (GenTHREADER) and standard sequence alignment at finding novel folds when sequences show no detectable homology to proteins with known structures.

Animals↗
PubMed2002

Motif prediction in ribosomal RNAs Lessons and prospects for automated motif prediction in homologous RNA molecules.

The traditional way to infer RNA secondary structure involves an iterative process of alignment and evaluation of covariation statistics between all positions possibly involved in basepairing. Watson-Crick basepairs typically show covariations that score well when examples of two or more possible basepairs occur. This is not necessarily the case for non-Watson-Crick basepairing geometries. For example, for sheared (trans Hoogsteen/Sugar edge) pairs, one base is highly conserved (always A or mostly A with some C or U), while the other can vary (G or A and sometimes C and U as well). RNA motifs consist of ordered, stacked arrays of non-Watson-Crick basepairs that in the secondary structure representation form hairpin or internal loops, multi-stem junctions, and even pseudoknots. Although RNA motifs occur recurrently and contribute in a modular fashion to RNA architecture, it is usually not apparent which bases interact and whether it is by edge-to-edge H-bonding or solely by stacking interactions. Using a modular sequence-analysis approach, recurrent motifs related to the sarcin-ricin loop of 23S RNA and to loop E from 5S RNA were predicted in universally conserved regions of the large ribosomal RNAs (16S- and 23S-like) before the publication of high-resolution, atomic-level structures of representative examples of 16S and 23S rRNA molecules in their native contexts. This provides the opportunity to evaluate the predictive power of motif-level sequence analysis, with the goal of automating the process for predicting RNA motifs in genomic sequences. The process of inferring structure from sequence by constructing accurate alignments is a circular one. The crucial link that allows a productive iteration of motif modeling and realignment is the comparison of the sequence variations for each putative pair with the corresponding isostericity matrix to determine which basepairs are consistent both with the sequence and the geometrical data.

Base Pairing↗
PubMed2000

The chicken genome contains no HMG1 retropseudogenes but a functional HMG1 gene with long introns.

We have cloned the genomic sequence coding for the high mobility group 1 (HMG1) protein in chickens. Multiple sequence alignment shows that the chicken HMG1 gene is highly homologous to the human and the mouse HMG1 genes. The gene structure of chicken HMG1 is similar to that of the mouse and the human HMG1 genes, with the same exon-intron boundaries. However, in contrast to other avian genes that have shorter introns, the chicken HMG1 gene has introns that are twice as long as their mammalian homologues. In addition to the functional, intron-containing HMG1 gene, all mammalian genomes contain more than 50 copies of HMG1 retropseudogenes each, while in the chicken genome there are no HMG1 retropseudogenes. This finding suggests that the HMG1 retropseudogenes arose in mammals after their divergence away from the birds.

Amino Acid Sequence↗
PubMed2004

Finding functional sites in structural genomics proteins.

Assigning function to structures is an important aspect of structural genomics projects, since they frequently provide structures for uncharacterized proteins. Similarities uncovered by structure alignment can suggest a similar function, even in the absence of sequence similarity. For proteins adopting novel folds or those with many functions, this strategy can fail, but functional clues can still come from comparison of local functional sites involving a few key residues. Here we assess the general applicability of functional site comparison through the study of 157 proteins solved by structural genomics initiatives. For 17, the method bolsters confidence in predictions made based on overall fold similarity. For another 12 with new folds, it suggests functions, including a putative phosphotyrosine binding site in the Archaeal protein Mth1187 and an active site for a ribose isomerase. The approach is applied weekly to all new structures, providing a resource for those interested in using structure to infer function.

Algorithms↗
PubMed2004

Identifying gene regulatory elements by genome-wide recovery of DNase hypersensitive sites.

Analysis of the human genome sequence has identified approximately 25000-30000 protein-coding genes, but little is known about how most of these are regulated. Mapping DNase I hypersensitive (HS) sites has traditionally represented the gold-standard experimental method for identifying regulatory elements, but the labor-intensive nature of this technique has limited its application to only a small number of human genes. We have developed a protocol to generate a genome-wide library of gene regulatory sequences by cloning DNase HS sites. We generated a library of DNase HS sites from quiescent primary human CD4(+) T cells and analyzed approximately 5600 of the resulting clones. Compared to sequences from randomly generated in silico libraries, sequences from these clones were found to map more frequently to regions of the genome known to contain regulatory elements, such as regions upstream of genes, within CpG islands, and in sequences that align between mouse and human. These cloned sites also tend to map near genes that have detectable transcripts in CD4(+) T cells, demonstrating that transcriptionally active regions of the genome are being selected. Validation of putative regulatory elements was achieved by repeated recovery of the same sequence and real-time PCR. This cloning strategy, which can be scaled up and applied to any cell line or tissue, will be useful in identifying regulatory elements controlling global expression differences that delineate tissue types, stages of development, and disease susceptibility.

CD4-Positive T-Lymphocytes↗
PubMed2003

Increased coverage obtained by combination of methods for protein sequence database searching.

MOTIVATION: Sequence alignment methods that compare two sequences (pairwise methods) are important tools for the detection of biological sequence relationships. In genome annotation, multiple methods are often run and agreement between methods taken as confirmation. In this paper, we assess the advantages of combining search methods by comparing seven pairwise alignment methods, including three local dynamic programming algorithms (PRSS, SSEARCH and SCANPS), two global dynamic programming algorithms (GSRCH and AMPS) and two heuristic approximations (BLAST and FASTA), individually and by pairwise intersection and union of their result lists at equal p-value cut-offs. RESULTS: When applied singly, the dynamic programming methods SCANPS and SSEARCH gave significantly better coverage (p=0.01) compared to AMPS, GSRCH, PRSS, BLAST and FASTA. Results ranked by BLAST p-values gave significantly better coverage compared to ranking by BLAST e-values. Of 56 combinations of eight methods considered, 19 gave significant increases in coverage at low error compared to the parent methods at an equal p-value cutoff. The union of results by BLAST (p-value) and FASTA at an equal p-value cutoff gave significantly better coverage than either method individually. The best overall performance was obtained from the intersection of the results from SSEARCH and the GSRCH62 global alignment method. At an error level of five false positives, this combination found 444 true positives, a significant 12.4% increase over SSEARCH applied alone.

Algorithms↗
PubMed1999

Structure and genomic organization of a second cluster of immunoglobulin heavy chain gene segments in the channel catfish.

The structure, organization, and partial sequence of a 25-kb genomic region containing a second cluster of H chain gene segments in the channel catfish (Ictalurus punctatus) has been determined. Multiple VH gene segments, representing different VH families, are located upstream of a germline-joined VDJ. The VDJ segment has a split leader sequence and a single open reading consistent with that expressed in members of the VH1 family. Downstream of the germline-joined VDJ is a single JH segment and two pseudogene exons structurally similar to the Cmu1 and Cmu2 exons of the functional gene. Both pseudogene exons are multiply crippled with RNA splice sites destroyed, and open reading frames are interrupted by termination codons, insertions, and/or deletions. Sequence alignment of a 10.8-kb region within the second H chain cluster with the genomic sequence of the nine JH segments and the functional Cmu within the first H chain gene cluster indicates that the second H chain gene cluster probably arose by a massive duplication event. The JH region of the VDJ, the coding and flanking regions of the single JH segment, and the pseudogene Cmu exons were readily aligned with homologous segments in the first gene cluster. This duplication event may have extended to include the upstream VH segments. A member of the Tc1 mariner family of transposable elements is located downstream of the pseudogene Cmu2, which suggests that the transposition may have contributed to the evolution of the duplicated Cmu.

Amino Acid Sequence↗
PubMed1992

Parallelization of a local similarity algorithm.

The local similarity problem is to determine the similar regions within two given sequences. We recently developed a dynamic programming algorithm for the local similarity problem that requires only space proportional to the sum of the two sequence lengths, whereas earlier methods use space proportional to the product of the lengths. In this paper, we describe how to parallelize the new algorithm and present results of experimental studies on an Intel hypercube. The parallel method provides rapid, high-resolution alignments for users of our software toolkit for pairwise sequence comparison, as illustrated here by a comparison of the chloroplast genomes of tobacco and liverwort.

Algorithms↗
PubMed2005

Use of genomic history to improve phylogeny and understanding of births and deaths in a gene family.

Polyploidy events have played an important role in the evolution of angiosperm genomes. Here, we demonstrate how genomic histories can increase phylogenetic resolution in a gene family, specifically the expansin superfamily of cell wall proteins. There are 36 expansins in Arabidopsis and 58 in rice. Traditional sequence-based phylogenetic trees yield poor resolution below the family level. To improve upon these analyses, we searched for gene colinearity (microsynteny) between Arabidopsis and rice genomic segments containing expansin genes. Multiple rounds of genome duplication and extensive gene loss have obscured synteny. However, by simultaneously aligning groups of up to 10 potentially orthologous segments from the two species, we traced the history of 49 out of 63 expansin-containing segments back to the ancestor of monocots and eudicots. Our results indicate that this ancestor had 15-17 expansin genes, each ancestral to an extant clade. Some clades have strikingly different growth patterns in the rice and Arabidopsis lineages, with more than half of all rice expansins arising from two ancestral genes. Segmental duplications, most of them part of polyploidy events, account for 12 out of 21 new expansin genes in Arabidopsis and 16 out of 44 in rice. Tandem duplications explain most of the rest. We were also able to estimate a minimum of 28 gene deaths in the Arabidopsis lineage and nine in rice. This analysis greatly clarifies expansin evolution since the last common ancestor of monocots and eudicots and the method should be broadly applicable to many other gene families.

Amino Acid Sequence↗
PubMed2005

Molecular systematics of the gonorynchiform fishes (Teleostei) based on whole mitogenome sequences: implications for higher-level relationships within the Otocephala.

Although the order Gonorynchiformes includes only 31 species assigned to seven genera and four families, it exhibits a large variety of anatomical structures, making difficult the reconstruction of phylogenetic relationships among its representatives. Within the basal teleosts, the Gonorynchiformes belong to the Otocephala where they have been alternatively placed as the sister group of the Otophysi and of the Clupeiformes. In this context, we investigated the phylogeny of the Gonorynchiformes using whole mitogenome sequences from 40 species (six being newly determined for this study). Our taxonomic sampling included at least one species of each gonorynchiform genus and of each other major otocephalan lineage. Unambiguously aligned, concatenated mitogenomic sequences (excluding the ND6 gene and control region) were divided into five partitions (1st, 2nd, and 3rd codon positions, tRNA genes, and rRNA genes) and partitioned Bayesian analyses were conducted. The resultant phylogenetic trees were fully resolved, with most of the nodes well supported by the high posterior probabilities. As expected, the Otocephala were recovered as monophyletic. Within this group, the mitogenome data supported the monophyly of Alepocephaloidei, Gonorynchiformes, Otophysi, and Clupeiformes. The Gonorynchiformes and the Otophysi formed a sister group, rending the Ostariophysi monophyletic. This result conflicts with previous mitogenomic phylogenetic studies, in which a sister relationship was found between Clupeiformes and Gonorynchiformes. We discussed the possible causes of this incongruence. Within the Gonorynchiformes, the following original topology was found: (Gonorynchus (Chanos (Phractolaemus (Cromeria (Grasseichthys (Kneria, Parakneria)))))). We confirmed that the paedomorphic species Cromeria nilotica and Grasseichthys gabonensis belong to the family Kneriidae; however, the two species together did not form a monophyletic group. This result challenges the value of reductive or absent characters as synapomorphies in this group.

Animals↗
PubMed2006

Lineage structures in the genome sequences of three Epstein-Barr virus strains.

Whole genome sequences for three Epstein-Barr virus strains (B95-8, GD1 and AG876) were aligned and compared. In addition to known variable loci (including type-specific alleles for the EBNA2, EBNA3A, EBNA3B and EBNA3C genes, plus the EBNA1 and LMP1 genes), seven large-scale regions of lower-level diversity were identified with strains at each in two major groupings. All three possible patterns of strain associations were represented across the seven loci. Tree-building studies supported the existence of two distinct lineages in each case, and occurrence of recombination between lineages therefore has to be invoked to account for the observed genotypes of virus strains.

Base Sequence↗
PubMed2004

Duplication of DYS19 flanking regions in other parts of the Y chromosome.

During the testing of alternative primers for the Y chromosome short tandem repeat marker DYS19, a duplicated region of the Y chromosome was discovered. The duplicated sequence is contained within GenBank accession AC006335 and has a high degree of homology with the DYS19 flanking region (GenBank accession AC017019) but without the polymorphic TAGA repeat. Bioinformatic approaches have been taken to try and understand the implications of this homolog to enable improved primer design for DYS19. Sequence alignments and careful placement of primers in order to obtain specific amplification of the DYS19 locus are discussed in the context of all previously published primer sets. Since the DYS19 locus is part of the widely used minimal haplotype, its robust amplification is highly desirable particularly in multiplex reactions. The discovery of this duplicated region of the Y chromosome shows the value of newly available human genome sequence information for assay design and the importance of using sequence queries and alignments in the primer design process.

5' Flanking Region↗
PubMed1998

Sequence polymorphism in two novel Plasmodium vivax ookinete surface proteins, Pvs25 and Pvs28, that are malaria transmission-blocking vaccine candidates.

BACKGROUND: For many malarious regions outside of Africa, development of effective transmission-blocking vaccines will require coverage against both Plasmodium falciparum and P. vivax. Work on P. vivax transmission-blocking vaccines has been hampered by the inability to clone the vaccine candidate genes from this parasite. MATERIALS AND METHODS: To search for genes encoding the ookinete surface proteins from P. vivax, the DNA sequences of the eight known proteins in the P25 subfamily (Pfs25, Pgs25, Pys25, Pbs25) and in the P21/28 subfamily (Pfs28, Pgs28, Pys21, Pbs21) of zygote/ookinete surface proteins were aligned. Regions of highest identity were used to design degenerate PCR oligonucleotides. Genomic DNA from the Sal I strain of P. vivax and genomic and splinkerette DNA libraries were used as PCR templates. To characterize the polymorphisms of Pvs25 and Pvs28, these two genes were PCR amplified and the DNA sequences were determined from genomic DNA extracted from patients infected with P. vivax. RESULTS: Analysis of the deduced amino acid sequence of Pvs28 revealed a secretory signal sequence, four epidermal growth factor (EGF)-like domains, six copies of the heptad amino acid repeat (GSGGE/D), and a short hydrophobic region. Because the fourth EGF-like domain has four rather than six cysteines, the gene designated Pvs28 is the presumed homologue of P21/28 subfamily members. Analysis of the deduced amino acid sequence of Pvs25 revealed a similar structure to that of Pvs28. The presence of six rather than four cysteines in the fourth EGF-like domain suggested that Pvs25 is the homologue of P25 subfamily members. Several regions of genetic polymorphisms in Pvs25 and Pvs28 were identified in field isolates of P. vivax. CONCLUSIONS: The genes encoding two ookinete surface proteins, Pvs28 and Pvs25, from P. vivax have been isolated and sequenced. Comparison of the primary structures of Pvs25, Pvs28, Pfs25, and Pfs28 suggest that there are regions of genetic polymorphism in the P25 and P21/28 subfamilies.

Amino Acid Sequence↗