Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “indel”

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 55 records · Page 3Linked to original sources

Simultaneous genotyping of indels and SNPs by mass spectroscopy.

Nucleotide insertion/deletion polymorphisms (indels) in ApoE gene were precisely genotyped using artificial ribonucleases and MALDI-TOF MS. The RNA fragments for MS analysis were prepared by treating RNA specimens with our artificial ribonucleases, which consist of LuCl(3) (molecular scissors) and oligonucleotides bearing two acridine groups (RNA-activator for site-selective scission). RNA scission by Lu(III) ion always occurred at the phosphodiester linkages in front of the two acridines, even when the RNA specimens involved consecutive cytidine sequences of different lengths. Thus, even complicated mixtures of these indel specimens were completely genotyped by using only one acridine-bearing oligonucleotide and by subjecting the reaction mixture to single MS measurement. Moreover, single nucleotide polymorphism (SNP) in the consecutive sequences could be genotyped simultaneously with the indels.

Apolipoprotein E3↗

Mechanism of a genetic glissando: structural biology of indel mutations.

Insertions and deletions of bases in DNA (collectively termed "indels") are both common and biologically relevant, being associated with different human pathologies including cancer and diseases associated with expansions of repeats. Four decades of research have resulted in several hypotheses regarding how indels are generated during DNA synthesis and how they subsequently undergo or escape correction. Recent structural studies of DNA polymerases bound to mutagenic substrates have increased our understanding of how DNA polymerases cope with abnormal substrates. These structures provide insight into the molecular mechanisms underlying indel generation.

DNA↗

Using evolutionary Expectation Maximization to estimate indel rates.

MOTIVATION: The Expectation Maximization (EM) algorithm, in the form of the Baum-Welch algorithm (for hidden Markov models) or the Inside-Outside algorithm (for stochastic context-free grammars), is a powerful way to estimate the parameters of stochastic grammars for biological sequence analysis. To use this algorithm for multiple-sequence evolutionary modelling, it would be useful to apply the EM algorithm to estimate not only the probability parameters of the stochastic grammar, but also the instantaneous mutation rates of the underlying evolutionary model (to facilitate the development of stochastic grammars based on phylogenetic trees, also known as Statistical Alignment). Recently, we showed how to do this for the point substitution component of the evolutionary process; here, we extend these results to the indel process. RESULTS: We present an algorithm for maximum-likelihood estimation of insertion and deletion rates from multiple sequence alignments, using EM, under the single-residue indel model owing to Thorne, Kishino and Felsenstein (the 'TKF91' model). The algorithm converges extremely rapidly, gives accurate results on simulated data that are an improvement over parsimonious estimates (which are shown to underestimate the true indel rate), and gives plausible results on experimental data (coronavirus envelope domains). Owing to the algorithm's close similarity to the Baum-Welch algorithm for training hidden Markov models, it can be used in an 'unsupervised' fashion to estimate rates for unaligned sequences, or estimate several sets of rates for sequences with heterogenous rates. AVAILABILITY: Software implementing the algorithm and the benchmark is available under GPL from http://www.biowiki.org/

Algorithms↗

Positive selection for indel substitutions in the rodent sperm protein catsper1.

Catsper1 is a voltage-gated calcium channel located in the plasma membrane of the sperm tail and is necessary for sperm motility and fertility in mice. We here examine the evolutionary pattern of Catsper1 from nine species of the rodent subfamily Murinae of family Muridae. We show that the rate of insertion/deletion (indel) substitutions in exon 1 of the gene is 4-15 times that in introns or neutral genomic regions, suggesting the presence of strong positive selection that promotes fixations of indel mutations in exon 1. The number of indel polymorphisms within species appears higher than expected from interspecific comparisons, although there are too little data to provide a statistically significant conclusion. These results, together with an earlier report in primates, indicate that positive selection promoting length variation in Catsper1 may be widespread in mammals. A structural model of Catsper1 suggested the importance of the exon 1-encoded region in regulating channel inactivation, which may affect sperm mobility and sperm competition. Our findings provide a necessary foundation for future experimental investigations of Catsper1's function in sperm physiology and role in sperm competition using rodent models.

Amino Acid Sequence↗

SIMPROT: using an empirically determined indel distribution in simulations of protein evolution.

BACKGROUND: General protein evolution models help determine the baseline expectations for the evolution of sequences, and they have been extensively useful in sequence analysis and for the computer simulation of artificial sequence data sets. RESULTS: We have developed a new method of simulating protein sequence evolution, including insertion and deletion (indel) events in addition to amino-acid substitutions. The simulation generates both the simulated sequence family and a true sequence alignment that captures the evolutionary relationships between amino acids from different sequences. Our statistical model for indel evolution is based on the empirical indel distribution determined by Qian and Goldstein. We have parameterized this distribution so that it applies to sequences diverged by varying evolutionary times and generalized it to provide flexibility in simulation conditions. Our method uses a Monte-Carlo simulation strategy, and has been implemented in a C++ program named Simprot. CONCLUSION: Simprot will be useful for testing methods of analysis of protein sequence families particularly alignment methods, phylogenetic tree building, detection of recombination and horizontal gene transfer, and homology detection, where knowing the true course of sequence evolution is essential.

Amino Acid Substitution↗

MCALIGN2: faster, accurate global pairwise alignment of non-coding DNA sequences based on explicit models of indel evolution.

BACKGROUND: Non-coding DNA sequences comprise a very large proportion of the total genomic content of mammals, most other vertebrates, many invertebrates, and most plants. Unraveling the functional significance of non-coding DNA depends on how well we are able to align non-coding DNA sequences. However, the alignment of non-coding DNA sequences is more difficult than aligning protein-coding sequences. RESULTS: Here we present an improved pair-hidden-Markov-Model (pair HMM) based method for performing global pairwise alignment of non-coding DNA sequences. The method uses an explicit model of indel length frequency distribution which can be specified, and allows any time reversible model of nucleotide substitution. The method uses a deterministic global optimiser to find the alignment with the highest posterior probability. We test MCALIGN2 in simulations, and compare it to a previous Monte Carlo based method (MCALIGN), to the pair HMM method of Knudsen and Miyamoto, and to a heuristic method (AVID) that performed very well in a previous simulation study. We show that the pair HMM methods have excellent performance for all combinations of parameter values we have considered. MCALIGN2 is up to ten times faster than MCALIGN. MCALIGN2 is more accurate in resolving indels given an accurate explicit model than heuristic methods, but is computationally slower. CONCLUSION: MCALIGN2 produces better quality alignments by explicitly using biological knowledge about the indel length distribution and time reversible models of nucleotide substitution. As a result, it can outperform other available sequence alignment methods for the cases we have considered to align non-coding DNA sequences.

Algorithms↗

Incorporating indel channels into average-case analysis of seed-chain-extend.

MOTIVATION: Given a sequence s1 of n letters drawn independently and identically (i.i.d.) from an alphabet of size &#x3c3; and a mutated substring s2 of length m<n, we want to recover the mutation history that generated s2 from s1. Many modern sequence aligners for this task use seed-chain-extend with k-mer seeds. Previously, Shaw and Yu showed linear-gap cost chaining can produce a chain with 1-O(1m) recoverability, the proportion of the mutation history that is recovered, in O(mn2.43&#x3b8;&#x2009;log&#x2009;n) expected time for seed-chain-extend (assuming pre-seeded reference), where &#x3b8;<0.206 is the mutation rate under a substitution-only channel and s1 is uniformly random. A gap remains between theory and practice, as real genomes include insertions and deletions (indels). RESULTS: We introduce mathematical machinery to deal with the two new obstacles introduced by indel channels: the dependence of neighbouring anchors and the presence of anchors that are only partially correct. We prove that expected recoverability of an optimal chain is &#x2265;1-O(1m) and expected runtime is O(mn3.15&#xb7;&#x3b8;T&#x2009;log&#x2009;n), given the total mutation rate &#x3b8;T=&#x3b8;i+&#x3b8;d+&#x3b8;s (sum of substitution, insertion, and deletion rates) is &#x3b8;T&#x2264;0.159. We thus narrow (but not close) the gap between theory and practice. AVAILABILITY AND IMPLEMENTATION: https://github.com/Lazarus42/seed_chainer_indels.

INDEL Mutation↗

Indels and imperfect duplication have driven the evolution of human Complement Receptor 1 (CR1) and CR1-like from their precursor CR1 alpha: importance of functional sets.

This study examines the effects of duplication and insertions-deletions (indels) by comparing human complement receptor 1 (CR1) and human CR1-like (CR1L) with syntenic genes from four other vertebrates (chimpanzee, baboon, rat, and mouse). By phylogenetic analysis, the domains of these genes can be classified into 10 distinct subfamilies (a, b, c, d, e, f, g(-like), h, j, and k), which have been largely conserved throughout vertebrate and invertebrate evolution. In spite of many complex and diverse duplications and indels, the subfamily order of domains (a, j, e, f, b, k, d, g(-like)) has been maintained. The number of domain sets has increased progressively, thereby expanding the functional repertoire.

Amino Acid Sequence↗

Comparative radiation hybrid map of canine chromosome 1 incorporating SNP and indel polymorphisms.

We report a comparative map of canine chromosome 1 (CFA1) incorporating single nucleotide polymorphisms (SNPs) and insertion/deletion (indel) polymorphisms, developed by using cross-species primers, radiation hybrid analysis, and pool-and-sequence identification of genetic variations. Fifty-five genes were chosen with relatively even spacing (approximately 3 Mb between the human homologues) and were mapped to CFA1, with 49 of these being new assignments. Evolutionary chromosomal breakpoints between CFA1 and the corresponding human chromosomes (HSA6, HSA9, HSA18, and HSA19) were located within 1 to 5 Mb based upon the human genome sequence. The process of identifying the evolutionary chromosomal breakpoints between CFA1 and the relevant human chromosomes led to an improvement in the comparative maps of CFA7, CFA12, and CFA29 through the mapping of 21 additional genes. A manual pool-and-sequence method was used to identify 79 SNPs, 9 small indels, 7 simple tandem repeats, and 2 polymorphic SINE insertions within the genes mapped. The cross-species primers can also be used in the manner described here to improve the comparative maps for other mammalian species.

Animals↗

In polymorphic genomic regions indels cluster with nucleotide polymorphism: Quantum Genomics.

Previously, we have described polymorphic frozen blocks (PFBs) within the Major Histocompatibility Complex (MHC) as regions of several hundred kilobases characterised by high nucleotide diversity, little or no recombination, duplicated segments, disease susceptibility, and human endogenous retroviruses. The nucleotide diversity profile within these PFBs shows peaks and troughs outside of the Class I genes, reflecting other important genes (or sequences) in the region. Here we show that indel density is also clustered with similar peaks and troughs. In fact, SNPs and indels are co-located within PFBs.

Databases, Nucleic Acid↗

A "Long Indel" model for evolutionary sequence alignment.

We present a new probabilistic model of sequence evolution, allowing indels of arbitrary length, and give sequence alignment algorithms for our model. Previously implemented evolutionary models have allowed (at most) single-residue indels or have introduced artifacts such as the existence of indivisible "fragments." We compare our algorithm to these previous methods by applying it to the structural homology dataset HOMSTRAD, evaluating the accuracy of (1) alignments and (2) evolutionary time estimates. With our method, it is possible (for the first time) to integrate probabilistic sequence alignment, with reliability indicators and arbitrary gap penalties, in the same framework as phylogenetic reconstruction. Our alignment algorithm requires that we evaluate the likelihood of any specific path of mutation events in a continuous-time Markov model, with the event times integrated out. To this effect, we introduce a "trajectory likelihood" algorithm (Appendix A). We anticipate that this algorithm will be useful in more general contexts, such as Markov Chain Monte Carlo simulations.

Algorithms↗

Molecular signatures (unique proteins and conserved indels) that are specific for the epsilon proteobacteria (Campylobacterales).

BACKGROUND: The epsilon proteobacteria, which include many important human pathogens, are presently recognized solely on the basis of their branching in rRNA trees. No unique molecular or biochemical characteristics specific for this group are known. RESULTS: Comparative analyses of proteins in the genomes of Wolinella succinogenes DSM 1740 and Campylobacter jejuni RM1221 against all available sequences have identified a large number of proteins that are unique to various epsilon proteobacteria (Campylobacterales), but whose homologs are not detected in other organisms. Of these proteins, 49 are uniquely found in nearly all sequenced epsilon-proteobacteria (viz. Helicobacter pylori (26695 and J99), H. hepaticus, C. jejuni (NCTC 11168, RM1221, HB93-13, 84-25, CF93-6, 260.94, 11168 and 81-176), C. lari, C. coli, C. upsaliensis, C. fetus, W. succinogenes DSM 1740 and Thiomicrospira denitrificans ATCC 33889), 11 are unique for the Wolinella and Helicobacter species (i.e. Helicobacteraceae family) and many others are specific for either some or all of the species within the Campylobacter genus. The primary sequences of many of these proteins are highly conserved and provide novel resources for diagnostics and therapeutics. We also report four conserved indels (i.e. inserts or deletions) in widely distributed proteins (viz. B subunit of exinuclease ABC, phenylalanyl-tRNA synthetase, RNA polymerase beta '-subunit and FtsH protein) that are specific for either all epsilon proteobacteria or different subgroups. In addition, a rare genetic event that caused fusion of the genes for the largest subunits of RNA polymerase (rpoB and rpoC) in Wolinella and Helicobacter is also described. The inter-relationships amongst Campylobacterales as deduced from these molecular signatures are in accordance with the phylogenetic trees based on the 16S rRNA and concatenated sequences for nine conserved proteins. CONCLUSION: These molecular signatures provide novel tools for identifying and circumscribing species from the Campylobacterales order and its subgroups in molecular terms. Although sequence information for these signatures is presently limited to Campylobacterales species, it is likely that many of them will also be found in other epsilon proteobacteria. Functional studies on these proteins and conserved indels should reveal novel biochemical or physiological characteristics that are unique to these groups of epsilon proteobacteria.

Bacterial Proteins↗

A Plasmodium whole-genome synteny map: indels and synteny breakpoints as foci for species-specific genes.

Whole-genome comparisons are highly informative regarding genome evolution and can reveal the conservation of genome organization and gene content, gene regulatory elements, and presence of species-specific genes. Initial comparative genome analyses of the human malaria parasite Plasmodium falciparum and rodent malaria parasites (RMPs) revealed a core set of 4,500 Plasmodium orthologs located in the highly syntenic central regions of the chromosomes that sharply defined the boundaries of the variable subtelomeric regions. We used composite RMP contigs, based on partial DNA sequences of three RMPs, to generate a whole-genome synteny map of P. falciparum and the RMPs. The core regions of the 14 chromosomes of P. falciparum and the RMPs are organized in 36 synteny blocks, representing groups of genes that have been stably inherited since these malaria species diverged, but whose relative organization has altered as a result of a predicted minimum of 15 recombination events. P. falciparum-specific genes and gene families are found in the variable subtelomeric regions (575 genes), at synteny breakpoints (42 genes), and as intrasyntenic indels (126 genes). Of the 168 non-subtelomeric P. falciparum genes, including two newly discovered gene families, 68% are predicted to be exported to the surface of the blood stage parasite or infected erythrocyte. Chromosomal rearrangements are implicated in the generation and dispersal of P. falciparum-specific gene families, including one encoding receptor-associated protein kinases. The data show that both synteny breakpoints and intrasyntenic indels can be foci for species-specific genes with a predicted role in host-parasite interactions and suggest that, besides rearrangements in the subtelomeric regions, chromosomal rearrangements may also be involved in the generation of species-specific gene families. A majority of these genes are expressed in blood stages, suggesting that the vertebrate host exerts a greater selective pressure than the mosquito vector, resulting in the acquisition of diversity.

Animals↗

Indel mutation in transcription factor PabHLH2 regulates amygdalin accumulation and kernel bitterness in apricot.

Amygdalin, the phytochemical responsible for the characteristic bitterness of apricot (Prunus armeniaca L.) kernels, also exhibits significant bioactive properties and therapeutic potential. Genetic regulation of amygdalin content is therefore a key objective in apricot breeding programs aimed at quality improvement. In this study, we conducted quantitative trait loci (QTL) mapping to uncover the genetic basis of sweet-bitter differentiation in apricot kernels. We identified a 15-bp insertion/deletion (indel) polymorphism strongly related to kernel bitterness, with marker validation achieving 100% concordance across 601 apricot germplasm accessions. Notably, this polymorphic site is located within the helix-loop-helix (HLH) domain of the basic HLH (bHLH) transcription factor PabHLH2. Protein interaction analyses revealed that the 15-bp deletion variant impaired dimerization capacity, reducing transcriptional activation of downstream targets. Using yeast one-hybrid screening and dual-luciferase reporter assays, we identified PaCYP71AN24 and PaCYP79D16 as direct transcriptional targets of PabHLH2. Functional characterization further indicated that the PabHLH2a variant (harboring the 15-bp insertion) significantly enhanced the promoter activity of these cytochrome P450 genes compared with the deletion variant. Transient overexpression and silencing experiments in apricot kernels further confirmed that the 15-bp insertion positively regulates both PaCYP71AN24/PaCYP79D16 expression and prunasin accumulation, the immediate biosynthetic precursor of amygdalin. Overall, these findings provide mechanistic insights into the allelic variation underlying kernel bitterness and delineate the molecular cascade of amygdalin biosynthesis. The identified molecular markers and functional characterization establish a basis for marker-assisted breeding of low-amygdalin apricot cultivars, supporting the dual-purpose utilization of kernels in food and pharmaceutical industries.

Amygdalin↗

Algorithms to reconstruct past indels: The deletion-only parsimony problem.

Ancestral sequence reconstruction is an important task in bioinformatics, with applications ranging from protein engineering to the study of genome evolution. When sequences can only undergo substitutions, optimal reconstructions can be efficiently computed using well-known algorithms. However, accounting for indels in ancestral reconstructions is much harder. First, for biologically-relevant problem formulations, no polynomial-time exact algorithms are available. Second, multiple reconstructions are often equally parsimonious or likely, making it crucial to correctly display uncertainty in the results. Here, we consider a parsimony approach where only deletions are allowed, while addressing the aforementioned limitations. First, we describe an exact algorithm to obtain all the optimal solutions. The algorithm runs in polynomial time if only one solution is sought. Second, we show that all possible optimal reconstructions for a fixed node can be represented using a graph computable in polynomial time. While previous studies have proposed graph-based representations of ancestral reconstructions, this result is the first to offer a solid mathematical justification for this approach. Finally we provide arguments for the relevance of the deletion-only case for the general case.

Algorithms↗

Selective targeting of indel-inferred differences in spatial structures of highly homologous proteins.

Recent findings have shown that the protein elongation factor-1alpha (EF-1alpha) from the eukaryotic pathogen Leishmania donovani possesses virulence properties. This was unexpected, since it has greater than 80% sequence identity with its human homologue. Given that EF-1alpha is essential for cell survival, in principle, it can be considered an attractive drug target. However, the challenge is to be able to selectively target the protein so as not to affect function of the human homologue. While a limited number of discrete differences were scattered throughout the sequence, most of the difference between these 2 homologues could be attributed to a 12-amino acid insert present in human EF-1alpha and absent from the leishmania sequence. In the present study, we modeled the spatial differences in structures of human and L. donovani EF-1alpha's inferred by this insertion-deletion (or "indel"). The protein models were used to develop antibodies directed specifically toward the deletion region of the pathogen protein. The strategy described allowed successful selective targeting of this putative leishmania virulence factor while avoiding recognition of the highly similar human EF-1alpha homologue. These findings may establish a new strategy for the development of antagonists directed against certain pathogenic targets having close human homologues.

Amino Acid Sequence↗

Chlamydiae-specific proteins and indels: novel tools for studies.

Chlamydiae species are important human and animal pathogens. Their obligate intracellular mode of replication has precluded the use of genetic and molecular biological approaches for understanding their biology. Comparative genomics have identified many rare genetic changes consisting of whole proteins and conserved indels (i.e. inserts or deletions) in widely distributed proteins that are distinctive characteristics of either all, or various subgroups within, chlamydiae. Additionally, several interesting cases of the lateral transfer of genes from free-living bacteria to a common ancestor of chlamydiae, and from chlamydiae to Trypanosoma/Leishmania, have been identified. These novel signatures have possible applications for advancing our understanding of the chlamydiae.

Bacterial Proteins↗

An empirical test of the treatment of indels during optimization alignment based on the phylogeny of the genus Secale (Poaceae).

The ability of the program POY, implementing optimization alignment, to deal with major indels is explored and discussed in connection with a phylogenetic analysis of the genus Secale based on partial Adh1 sequences. The Adh1 sequences used span exon 2-4. Nearly all variation is found in intron 2 and intron 3, which form the basis for the phylogenetic analyses. Both in some ingroup and outgroup taxa intron 3 has a major duplication. Previous phylogenetic analyses have repeatedly confirmed monophyly of both Secale and Hordeum, the latter being part of the outgroup. However, optimization alignment only recovers both genera as monophyletic when knowledge of the duplication is incorporated in the analysis. The phylogenetic relationships within Secale are not clearly resolved. Subspecific taxa of Secale strictum have identical sequences and they are confined to a monophyletic group. However, the two subspecific taxa of Secale cereale do not form a monophyletic group, and the position of Secale sylvestre is uncertain.

Algorithms↗