SEARCH · PubMed Health

Results for “INDEL Mutation”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9Linked to original sources

PubMed2005

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↗

PubMed2025

Analysis of genetic differences underlying chilling stress tolerance using whole genome Re-Sequencing in walnut (Juglans regia L.).

Walnut (Juglans regia L.) is prized worldwide for both its nutritional value and economic importance, yet it remains vulnerable to cold stress, with significant differences in tolerance among varieties. This study combined physiological analyses with whole-genome resequencing (WGS) to evaluate the cold stress responses of two varieties, ‘Qingxiang’ and ‘Liaoning No.8’. Under chilling stress (0 °C), we measured electrolyte leakage and antioxidant enzyme activity, applying both exogenous methyl jasmonate (MeJA) and the jasmonate inhibitor DIECA. Genomic variations were analyzed using WGS. Results showed that ‘Liaoning No.8’ exhibited superior cold tolerance. Application of MeJA reduced electrolyte leakage by 37% and MDA accumulation by 52% on average, whereas DIECA exacerbated stress-related damage. WGS achieved 16.24–16.26× coverage and identified 2.73–2.78 million SNPs, 378–382k InDels, 25–26k SVs, and 7.2–7.9k CNVs. Twenty genes containing sequence variants showed transcriptional responses under cold stress that were significantly correlated with mutation density (r = 0.62, P < 0.01). One gene, XM_018985465.2, which lacked SNPs in ‘Liaoning No.8’, was expressed 4.2 times more in this variety, suggesting cis-regulatory influence. These findings highlight the role of jasmonic acid signaling in enhancing cold tolerance in walnut and offer genomic insights into its underlying adaptive mechanisms.

Juglans↗

PubMed1997

Molecular evolution of angiosperm mitochondrial introns and exons.

Numbers of substitutions per site for 15 protein-coding genes and six introns of the plant mitochondria were estimated to compare modes and tempos of evolution between exons and introns, and numbers of insertions-deletions per site also were investigated in introns. Intra-gene homogeneity of numbers of substitutions per site was assessed further among different taxa and between mitochondrial and nuclear paralogs translocated from the mitochondrial genome. Gene-to-gene differences in numbers of substitutions per site were found to be higher for nonsynonymous than synonymous sites, and this could be due to differential selection if mutation rate is assumed constant for the genome. Some mitochondrial genes have evolved as fast as chloroplast genes, thus faster than previously thought. For coxI, relative rate tests showed that woody taxa evolved slower than annuals at synonymous sites. Generation time, population size, and speciation rate are likely factors involved in this rate heterogeneity. Introns were less constrained than their adjacent exons for both overall numbers of substitutions per site and indels, but, on average, overall numbers of substitutions per site for introns were similar to numbers of synonymous substitutions per site for exons. Correlations were generally high between numbers of substitutions and numbers of indels per site for the same intron. Mitochondrial genes transferred to the nucleus had an accelerated rate of substitution per site, which was most significant at synonymous sites. These differences between paralogs in two different genomes are likely the result of different mutation rates.

DNA, Mitochondrial↗

PubMed2006

Evidence that protein length expansion and contraction is partly due to mutational events in premeiotic cells.

Studies on the rate of evolution of proteins typically concentrate on rates of change of orthologous amino acids rather than on changes in size (i.e., generation of nonorthologous domains). Recent work has focused attention on Ser/Thr-rich regions in yeast as these tend to undergo size changes rapidly, with size polymorphisms commonly being found, especially in proteins with cell-surface localization. The underlying mechanism generating the indels is presently unclear though, due to a lack of correlation with the location of meiotic double-strand breaks, it has, by exclusion, been conjectured to be replication slippage. Here we provide new evidence to support this possibility. Notably, we show that Ser/Thr-rich repeat regions are more generally associated with the location of Mre11p in premeiotic cells. This is to be expected if the repeats were produced by mutational events in mitotic cells possibly through replication slippage.

Cell Cycle↗

PubMed1997

Nucleotide polymorphism in the acidic chitinase locus (ChiA) region of the wild plant Arabidopsis thaliana.

To investigate DNA variation in natural plant populations, a 1.8-kb region of the acidic chitinase locus (ChiA)was analyzed for 17 ecotypes of Arabidopsis thaliana sampled worldwide and 3 Arabis species in Japan. As in the Adh region, dimorphism was detected throughout the investigated ChiA region, suggesting the possibility that dimorphic DNA variation exists in the entire nuclear genome of A. thaliana. The ChiA region was divided into two blocks by an intragenic recombination between two parental sequence types, which diverged 7.4 MYA under the assumption that nucleotide mutation rate per site per year is mu = 10(-9). Nucleotide diversity in the entire ChiA region was 0.0104. Tajima's test was significantly negative for both nucleotide and indel variations, which was manifested as an excess of unique polymorphisms. However, the level and pattern of polymorphism in the ChiA region were inconsistent with simple theoretical explanations. The HKA test detected no difference in the levels of intra- and interspecific variations between the ChiA and Adh regions. In the ChiA coding region, no difference in the patterns of synonymous and replacement variation was found in intra- and interspecific comparisons by the MK test. Although it was difficult to determine the exact genetic mechanism acting on the ChA locus, these results suggested that the ChA locus region was under the same genetic mechanism before and after the establishment of A. thaliana as a species.

Alcohol Dehydrogenase↗

PubMed2006

"Genome design" model: evidence from conserved intronic sequence in human-mouse comparison.

Introns are shorter in housekeeping genes than in tissue- or development-specific genes. Differing explanations have been offered for this phenomenon: selection for economy (in housekeeping genes), mutation bias or "genomic design." The large-scale implementation in this present paper of a rigorous local sequence alignment algorithm revealed an unprecedented fraction of evolutionarily conserved DNA in human-mouse introns ( approximately 60% of human and approximately 70% of mouse intron length remained after masking for lineage-specific repeats). The length distributions of both conserved and nonconserved regions are very broad but show peaks close to nucleosomal and di-nucleosomal DNA. Both the fraction of conserved sequence and its absolute length were higher in introns of tissue-specific genes than housekeeping genes. This difference remained after control for between-species identity of the conserved fraction, mutation rate, and GC content. In a more direct control, the product of the conserved sequence fraction and the between-species identity of this fraction (which can be considered to be the fraction of conserved nucleotides) was greater in introns of tissue-specific genes than housekeeping genes. Neither the fraction of intron length covered by repeats nor the balance of small insertions and deletions (indels) can explain the greater length of introns in tissue-specific genes. The length of the conserved intronic DNA in a gene is correlated with the number of functional domains in the protein encoded by that gene. These results suggest that the greater length of introns in tissue-specific genes is not due to selection for economy or mutation bias but instead is related to functional complexity (probably mediated by chromatin condensation), and that the evolution of the bulk of noncoding DNA is not completely neutral.

Algorithms↗

PubMed2005

Recognizing the pseudogenes in bacterial genomes.

Pseudogenes are now known to be a regular feature of bacterial genomes and are found in particularly high numbers within the genomes of recently emerged bacterial pathogens. As most pseudogenes are recognized by sequence alignments, we use newly available genomic sequences to identify the pseudogenes in 11 genomes from 4 bacterial genera, each of which contains at least 1 human pathogen. The numbers of pseudogenes range from 27 in Staphylococcus aureus MW2 to 337 in Yersinia pestis CO92 (e.g. 1-8% of the annotated genes in the genome). Most pseudogenes are formed by small frameshifting indels, but because stop codons are A + T-rich, the two low-G + C Gram-positive taxa (Streptococcus and Staphylococcus) have relatively high fractions of pseudogenes generated by nonsense mutations when compared with more G + C-rich genomes. Over half of the pseudogenes are produced from genes whose original functions were annotated as 'hypothetical' or 'unknown'; however, several broadly distributed genes involved in nucleotide processing, repair or replication have become pseudogenes in one of the sequenced Vibrio vulnificus genomes. Although many of our comparisons involved closely related strains with broadly overlapping gene inventories, each genome contains a largely unique set of pseudogenes, suggesting that pseudogenes are formed and eliminated relatively rapidly from most bacterial genomes.

Genome, Bacterial↗

PubMed2005

[Genetic analysis and mapping of the rice leafy-hull mutant Oslh].

Oslh (lh=leafy hull), in the japonica cultivar 9522 background, a mutant of Oryza sativa L. spp. japonica cv. 9522 identified from an M(2) population, was mutagenized by irradiation with (60)Co gamma-ray. The Oslh mutant plants flowered about 15 days later than the wild-type plants (Fig.1e). The paleas, lemmas and lodicules of the flowers of Oslh mutant were transformed into leaf-like structures (Fig.1b, d). Genetic analysis of the F(2) progeny from a cross between the Oslh mutant and wild-type japonica cv. 9522 revealed that the Oslh mutant arouse from a single recessive nuclear gene mutation of the cv. 9522. To map the Oslh locus, an F(2) population generated by crossing between Oslh (japonica) mutant and Guangluai4 (indica) was analyzed. The Oslh locus was mapped to the long arm of rice chromosome 3, between a SSR marker RM5475 and an InDel marker GY305, 2.9 and 1.5 cM away from these two markers respectively (Fig.4). These results are useful for further cloning and functional analysis of the OsLH gene.

Chromosome Mapping↗

PubMed2003

The role of context-dependent mutations in generating compositional and codon usage bias in grass chloroplast DNA.

The influence of local base composition on mutations in chloroplast DNA (cpDNA) is studied in detail and the resulting, empirically derived, mutation dynamics are used to analyze both base composition and codon usage bias. A 4 x 4 substitution matrix is generated for each of the 16 possible flanking base combinations (contexts) using 17,253 noncoding sites, 1309 of which are variable, from an alignment of three complete grass chloroplast genome sequences. It is shown that substitution bias at these sites is correlated with flanking base composition and that the A+T content of these flanking sites as well as the number of flanking pyrimidines on the same strand appears to have general influences on substitution properties. The context-dependent equilibrium base frequencies predicted from these matrices are then applied to two analyses. The first examines whether or not context dependency of mutations is sufficient to generate average compositional differences between noncoding cpDNA and silent sites of coding sequences. It is found that these two classes of sites exist, on average, in very different contexts and that the observed mutation dynamics are expected to generate significant differences in overall composition bias that are similar to the differences observed in cpDNA. Context dependency, however, cannot account for all of the observed differences: although silent sites in coding regions appear to be at the equilibrium predicted, noncoding cpDNA has a significantly lower A+T content than expected from its own substitution dynamics, possibly due to the influence of indels. The second study examines the codon usage of low-expression chloroplast genes. When context is accounted for, codon usage is very similar to what is predicted by the substitution dynamics of noncoding cpDNA. However, certain codon groups show significant deviation when followed by a purine in a manner suggesting some form of weak selection other than translation efficiency. Overall, the findings indicate that a full understanding of mutational dynamics is critical to understanding the role selection plays in generating composition bias and sequence structure.

Base Composition↗

PubMed2003

Evolutionary rate variation in anthocyanin pathway genes.

Over a broad taxonomic range that spans monocots and dicots, upstream enzymes of the anthocyanin pigment pathway have evolved less rapidly than downstream enzymes. In this article we show that this pattern is also evident within the genus Ipomoea. Specifically, the most upstream enzyme, chalcone synthase (CHS-D), evolves more slowly than the two most downstream enzymes, ancyocyanidin synthase (ANS) and UDP glucose flavonoid 3-oxy-glucosyltransferase (UFGT). This pattern appears not to be due to variation in mutation rates, because the CHS-D gene exhibits higher synonymous substitution rates than the genes for the other two enzymes. Codon-based tests for positive selection suggest that it has been negligible or absent in all three genes. In addition, the mean number of indel-creating events is four times as high in the downstream genes as in CHS-D. Unlike the downstream genes, CHS-D also exhibits evidence of codon bias. Together, the evidence suggests that the difference in nonsynonymous substitution rates between upstream and downstream genes is due to relaxed constraint on the downstream genes rather than a greater frequency of positively selected substitutions.

Acyltransferases↗

PubMed2004

[Nucleotide polymorphism in the drought induced transcription factor CBF4 region of Arabidopsis thaliana and its molecular evolution analyses].

Intraspecific nucleotide polymorphism in the drought induced transcription factor CBF4 region of Arabidopsis thaliana was analyzed with 17 core accessions growing in different ecoclimate. High density of single nucleotide polymorphism (SNP) and insertion/deletion (Indel) were found, on average 1 SNP per 35.8 bp and 1 Indel per 143 bp. Nucleotide polymorphism in non-coding region was three times higher than that in coding region. In coding region of CBF4, SNP frequency is one SNP per 96.4 bp, one nonsynonymous mutation was detected from 25 av, 203 av and 244 av accessions, which is the 205th site amino acid variation: gly <--> val caused by the 1034th site (corresponding to 19,696 site nucleotide of GenBank No. AB015478 as 1) nucleotide variation: G <--> T. Statistical result of nucleotide diversity showed that linkage disequilibrium (LD) existed in large-scale region of CBF4 and recombination event was also detected in 5' non-coding region. Identical to the results of other genes of Arabidopsis, different regions of the gene were seemingly under different selective pressures. Balancing selection resulted in high nucleotide diversity in 3' non-coding region, and the neutral mutation hypothesis can explain the DNA polymorphism in coding region, whereas, nature positive selection in the population affected nucleotide variation in 5' non-coding region of gene.

Arabidopsis↗

PubMed1996

Duplication and divergence of the genes of the alpha-esterase cluster of Drosophila melanogaster.

The alpha-esterase cluster of D. melanogaster contains 11 esterase genes dispersed over 60 kb. Embedded in the cluster are two unrelated open reading frames that have sequence similarity with genes encoding ubiquitin-conjugating enzyme and tropomyosin. The esterase amino acid sequences show 37-66% identity with one another and all but one have all the motifs characteristic of functional members of the carboxyl/cholinesterase multigene family. The exception has several frameshift mutations and appears to be a pseudogene. Patterns of amino acid differences among cluster members in relation to generic models of carboxyl/cholinesterase protein structure are broadly similar to those among other carboxyl/cholinesterases sequenced to date. However the alpha-esterases differ from most other members of the family in: their lack of a signal peptide; the lack of conservation in cysteines involved in disulfide bridges; and in four indels, two of which occur in or adjacent to regions that align with proposed substrate-binding sites of other carboxyl/cholinesterases. Phylogenetic analyses clearly identify three simple gene duplication events within the cluster. The most recent event involved the pseudogene which is located in an intron of another esterase gene. However, relative rate tests suggest that the pseudogene remained functional after the duplication event and has become inactive relatively recently. The distribution of indels also suggests a deeper node in the gene phylogeny that separates six genes at the two ends of the cluster from a block of five in the middle.

Amino Acid Sequence↗

PubMed2005

Mitochondrial genome sequences and molecular evolution of the Irish potato famine pathogen, Phytophthora infestans.

The mitochondrial genomes of haplotypes of the Irish potato famine pathogen, Phytophthora infestans, were sequenced. The genome sizes were 37,922, 39,870 and 39,840 bp for the type Ia, IIa and IIb mitochondrial DNA (mtDNA) haplotypes, respectively. The mitochondrial genome size for the type Ib haplotype, previously sequenced by others, was 37,957 bp. More than 90% of the genome contained coding regions. The GC content was 22.3%. A total of 18 genes involved in electron transport, 2 RNA-encoding genes, 16 ribosomal protein genes and 25 transfer RNA genes were coded on both strands with a conserved arrangement among the haplotypes. The type I haplotypes contained six unique open reading frames (ORFs) of unknown function while the type II haplotypes contained 13 ORFs of unknown function. Polymorphisms were observed in both coding and non-coding regions although the highest variation was in non-coding regions. The type I haplotypes (Ia and Ib) differed by only 14 polymorphic sites, whereas the type II haplotypes (IIa and IIb) differed by 50 polymorphic sites. The largest number (152) of polymorphic sites was found between the type IIb and Ia haplotypes. A large spacer flanked by the genes coding for tRNA-Tyr (trnY) and the small subunit RNA (rns) contained the largest number of polymorphic sites and corresponds to the region where a large indel that differentiates type II from type I haplotypes is located. The size of this region was 785, 2,666 and 2,670 bp in type Ia, IIa and IIb haplotypes, respectively. Among the four haplotypes, 81 mutations were identified. Phylogenetic and coalescent analysis revealed that although the type I and II haplotypes shared a common ancestor, they clearly formed two independent lineages that evolved independently. The type II haplotypes diverged earlier than the type I haplotypes. Thus our data do not support the previous hypothesis that the type II lineages evolved from the type I lineages. The type I haplotypes diverged more recently and the mutations associated with the evolution of the Ia and Ib types were identified.

Evolution, Molecular↗

PubMed2026

Dual-dimensional profiling of host genomic variations and HPV integration in PD-L1-stratified cervical cancer via Oxford Nanopore Technology.

BACKGROUND: The integration of human papillomavirus (HPV) DNA into the host genome is a key step in the development of HPV-associated cervical cancer (CC). However, the genomic characteristics of host genomic variations and HPV integration within the context of programmed death-ligand 1 (PD-L1) expression stratification have not been systematically investigated. METHODS: Whole-genome sequencing was performed using Oxford Nanopore Technology (ONT) on six samples (three from the high PD-L1 expression group and three from the low PD-L1 expression group). The characteristics of host genomic variations under different PD-L1 expression stratifications were explored, including structural variations (SV), copy number variations (CNV), single nucleotide polymorphisms (SNP), and insertion-deletions (Indel). Subsequently, the distribution features of HPV integration sites were analyzed, different integration types were identified, and pathway analysis was conducted. RESULTS: Whole-genome SV analysis revealed that the total number of SVs and the composition of mutation types were similar between the high and low PD-L1 expression groups, with insertions (INS) and deletions (DEL) predominating in both. These variations were primarily enriched in intergenic regions and introns. In the low PD-L1 expression group, integration events were observed at multiple chromosomal loci, with the most frequent integration occurring in the KLF5 gene region on chromosome 13. No frequently integrated loci were identified in the high PD-L1 expression group. Additionally, four distinct HPV integration breakpoint patterns were preliminarily identified and analyzed. CONCLUSION: PD-L1 expression stratification did not significantly alter the overall genomic instability of the host. However, differences were observed in the distribution patterns of HPV integration sites. These findings provide new insights into the genomic heterogeneity of CC under different PD-L1 expression backgrounds and may lay the groundwork for future research exploring stratified immunotherapy based on HPV integration features.

Humans↗

PubMed2005

Isolation and characterization of a rice mutant hypersensitive to Al.

Rice (Oryza sativa L.) is a highly Al-resistant species among small grain crops, but the mechanism responsible for the high Al resistance has not been elucidated. In this study, rice mutants sensitive to Al were isolated from M(3) lines derived from an Al-resistant cultivar, Koshihikari, irradiated with gamma-rays. Relative root elongation was used as a parameter for evaluating Al resistance. After initial screening plus two rounds of confirmatory testing, a mutant (als1) was isolated from a total of 560 lines. This mutant showed a phenotype similar to the wild-type plant in the absence of Al. However, in the presence of 10 microM Al, root elongation was inhibited 70% in the mutant, but only 8% in the wild-type plant. The mutant also showed poorer root growth in acid soil. The Al content of root apices (0-1 cm) was much lower in the wild-type plant. The sensitivity to other metals including Cd and La did not differ between the mutant and the wild-type plants. A small amount of citrate was secreted from the roots of the mutant in response to Al stress, but there was no difference from that secreted by the wild-type plant. Genetic analysis of F(2) populations between als1 and wild-type plants showed that the Al-resistant seedlings and Al-sensitive seedlings segregated at a 3 : 1 ratio, indicating that the high sensitivity to Al in als1 is controlled by a single recessive gene. The gene was mapped to the long arm of chromosome 6, flanked by InDel markers MaOs0619 and MaOs0615.

Aluminum↗

PubMed2001

Functions of the mismatch repair gene mutS from Acinetobacter sp. strain ADP1.

The genus Acinetobacter encompasses a heterogeneous group of bacteria that are ubiquitous in the natural environment due in part to their ability to adapt genetically to novel challenges. Acinetobacter sp. strain ADP1 (also known as strain BD413) is naturally transformable and takes up DNA from any source. Donor DNA can be integrated into the chromosome by recombination provided it possesses sufficient levels of nucleotide sequence identity to the recipient's DNA. In other bacteria, the requirement for sequence identity during recombination is partly due to the actions of the mismatch repair system, a key component of which, MutS, recognizes mismatched bases in heteroduplex DNA and, along with MutL, blocks strand exchange. We have cloned mutS from strain ADP1 and examined its roles in preventing recombination between divergent DNA and in the repair of spontaneous replication errors. Inactivation of mutS resulted in 3- to 17-fold increases in transformation efficiencies with donor sequences that were 8 to 20% divergent relative to the strain ADP1. Strains lacking MutS exhibited increased spontaneous mutation frequencies, and reversion assays demonstrated that MutS preferentially recognized transition mismatches while having little effect on the repair of transversion mismatches. Inactivation of mutS also abolished the marker-specific variations in transforming efficiency seen in mutS(+) recipients where transition and frameshift alleles transformed at eightfold lower frequencies than transversions or large deletions. Comparison of the MutS homologs from five individual Acinetobacter strains with those of other gram-negative bacteria revealed that a number of unique indels are conserved among the Acinetobacter amino acid sequences.

Acinetobacter↗

PubMed1999

Evolutionary implications of matK indels in Poaceae.

Insertion/deletion events (indels) and nucleotide substitutions at the extreme 3' end of the chloroplast gene matK have been identified that distinguish certain major lineages of grasses. A 1-bp (base pair) deletion creating a shift in the open reading frame (ORF) and a point mutation support the positions of Streptochaeta and Anomochloa as the two most basal lineages in Poaceae. Another 1-bp deletion resulting in early termination of the ORF is unique to Ehrharta, a member of the taxonomically disputable tribe Ehrharteae. A 6-bp insertion supports monophyly of subfamilies Panicoideae, Arundinoideae, Centothecoideae, and Chloridoideae (PACC). This marker appears useful in defining PACC clade members and may have potential in providing insight into the sister-group relationship between PACC and other lineages. Alignment of deduced amino acid sequences from bryophytes, gymnosperms, and angiosperms shows that this region is relatively conserved, but variation is notably higher in Poaceae. The evolutionary implications of these changes in grasses and other plant families are addressed.

Journal Article↗

PubMed2003

Rapid molecular evolution of CYCLOIDEA-like genes in Antirrhinum and its relatives.

The CYCLOIDEA (CYC) and DICHOTOMA (DICH) genes encode related TCP transcription factors that control floral asymmetry in Antirrhinum majus. Analysis of sequences from relatives of Antirrhinum suggested that CYC and DICH arose from a gene duplication in an ancestor of the tribe Antirrhineae and have subsequently evolved at similar rates. Coding regions outside the conserved functional TCP and R domains differed by numerous indels, suggesting rapid evolution and low constraint on amino acid sequence. An analysis of variability within the genus Antirrhinum revealed very similar CYC alleles in 17 representative species, consistent with most of the species having diverged within the last 1 myr. Whereas substitution mutations appear to have accumulated constantly, one Antirrhinum CYC allele provided evidence for sporadic and rapid accumulation of insertion mutations.

Antirrhinum↗