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

PubMed2000

Complex mutations in a high proportion of microsatellite loci from the protozoan parasite Plasmodium falciparum.

Microsatellite loci are generally assumed to evolve via a stepwise mutational process and a battery of statistical techniques has been developed in recent years based on this or related mutation models. It is therefore important to investigate the appropriateness of these models in a wide variety of taxa. We used two approaches to examine mutation patterns in the malaria parasite Plasmodium falciparum: (i) we examined sequence variation at 12 tri-nucleotide repeat loci; and (ii) we analysed patterns of repeat structure and heterozygosity at 114 loci using data from 12 laboratory parasite lines. The sequencing study revealed complex patterns of mutation in five of the 12 loci studied. Alleles at two loci contain indels of 24 bp and 57 bp in flanking regions, while in the other three loci, blocks of imperfect microsatellites appear to be duplicated or inserted; these loci essentially consist of minisatellite repeats, with each repeat unit containing four to eight microsatellites. The survey of heterozygosity revealed a positive relationship between repeat number and microsatellite variability for both di- and trinucleotides, indicating a higher mutation rate in loci with longer repeat arrays. Comparisons of levels of variation in different repeat types indicate that the mutation rate of dinucleotide-bearing loci is 1.6-2.1 times faster than trinucleotides, consistent with the lower mean number of repeats in trinucleotide-bearing loci. However, despite the evidence that microsatellite arrays themselves are evolving in a manner consistent with stepwise mutation model in P. falciparum, the high frequency of complex mutations precludes the use of analytical tools based on this mutation model for many microsatellite-bearing loci in this protozoan. The results call into question the generality of models based on stepwise mutation for analysing microsatellite data, but also demonstrate the ease with which loci that violate model assumptions can be detected using minimal sequencing effort.

Animals↗

PubMed2004

Spiking of contemporary human template DNA with ancient DNA extracts induces mutations under PCR and generates nonauthentic mitochondrial sequences.

Proof of authenticity is the greatest challenge in palaeogenetic research, and many safeguards have become standard routine in laboratories specialized on ancient DNA research. Here we describe an as-yet unknown source of artifacts that will require special attention in the future. We show that ancient DNA extracts on their own can have an inhibitory and mutagenic effect under PCR. We have spiked PCR reactions including known human test DNA with 14 selected ancient DNA extracts from human and nonhuman sources. We find that the ancient DNA extracts inhibit the amplification of large fragments to different degrees, suggesting that the usual control against contaminations, i.e., the absence of long amplifiable fragments, is not sufficient. But even more important, we find that the extracts induce mutations in a nonrandom fashion. We have amplified a 148-bp stretch of the mitochondrial HVRI from contemporary human template DNA in spiked PCR reactions. Subsequent analysis of 547 sequences from cloned amplicons revealed that the vast majority (76.97%) differed from the correct sequence by single nucleotide substitutions and/or indels. In total, 34 positions of a 103-bp alignment are affected, and most mutations occur repeatedly in independent PCR amplifications. Several of the induced mutations occur at positions that have previously been detected in studies of ancient hominid sequences, including the Neandertal sequences. Our data imply that PCR-induced mutations are likely to be an intrinsic and general problem of PCR amplifications of ancient templates. Therefore, ancient DNA sequences should be considered with caution, at least as long as the molecular basis for the extract-induced mutations is not understood.

Base Sequence↗

PubMed2000

Extensive homoplasy, nonstepwise mutations, and shared ancestral polymorphism at a complex microsatellite locus in Lake Malawi cichlids.

Recent studies have suggested that size homoplasy is a prevalent feature of microsatellites and is expected to increase with time of divergence among populations and taxa. In this study, we performed sequence analysis of alleles from a complex microsatellite locus (Pzeb4, initially isolated from Pseudotropheus (Maylandia) zebra) from 1 midwater-feeding and 10 rock-dwelling cichlid fish species from Lake Malawi, East Africa, to investigate how widespread size homoplasy is among closely related taxa at this locus. All cichlid fishes endemic to this lake are believed to have originated within the last 700,000 years, and some species may be less than 200 years old. The number of eletromorphs found per species varied from 3 to 13. Sequence analysis of 95 cloned Pzeb4 PCR products (representing 18 electromorphs) revealed 13 new alleles. Ten of the 13 electromorphs (77%) were found to show size homoplasy due to either single nucleotide substitutions/indels or large indels. To investigate how well this locus fits the single-step mutation model (SMM), the minimum number of mutations required to explain the length differences between pairs of alleles was plotted against their size differences. Of the 300 comparisons, 166 (55.3%) corresponded to SMM expectations and 86 (28.7%) required a smaller number of mutations, and for 48 (16.0%) pairwise comparisons, a larger number of mutations were required to explain the length differences as compared with SMM expectations. Finally, a large deletion in the microsatellite sequence observed in the three rock-dwelling species Pseudotropheus lucerna, Pseudotropheus (Tropheops) 'band,' and Pseudotropheus (Tropheops) 'rust' and the midwater-feeding species Copadichromis sp. is believed to represent a shared ancestral polymorphism.

Africa↗

PubMed2000

Hairpins involving both inverted and direct repeats are associated with homoplasious indels in non-coding chloroplast DNA of Taraxacum (Lactuceae: Asteraceae).

Sequence variation in 2.2 kb of non-coding regions of the chloroplast genome of eight dandelions (Taraxacum: Lactuceae) from Asia and Europe is interpreted in the light of the phylogenetic signal of base substitutions vs. indels (insertions-deletions). The four non-coding regions displayed a total of approximately 30 structural mutations of which 9 are potentially phylogenetically informative. Insertions, deletions, and an inversion were found that involved consecutive stretches of up to 172 bases. When compared to phylogenetic relationships of the chloroplast genomes based on nucleotide substitutions only, many homoplasious indels (33%) were detected that differed considerably in length and did not comprise simple sequence repeats typically associated with replication slippage. Though many indels in the intergenic spacers were associated with direct repeats, frequently, the variable stretches participated in inverted repeat stabilized hairpins. In each intergenic spacer or intron examined, nucleotide stretches ranging from 30 to 60 bp were able to fold into stabilized secondary structures. When these indels were homoplasious, they always ranked among the most stabilized hairpins in the non-coding regions. The association of higher order structures that involve both classes of repeats and parallel structural mutations in hot spot regions of the chloroplast genome can be used to differentiate among mutations that differ in phylogenetic reliability.

Asteraceae↗

PubMed2026

Safe and Stable Germline Transmission of MSTN Mutations in Cattle.

With the global population expected to reach 10 billion by 2050, sustainable livestock production is critical. Gene editing of the myostatin (MSTN) gene represents a promising strategy to enhance muscle growth in cattle. In this study, MSTN-mutated founder (F0) cows were used to generate F1 offspring via ovum pick-up, in vitro fertilization, and embryo transfer. Four F1 calves were born, all confirmed to be heterozygous for the MSTN mutation. Long-term monitoring showed normal growth and no visible health abnormalities. Whole-genome sequencing identified SNPs, INDELs, and structural variants, most with minimal predicted functional effects. Proteomic profiling of Longissimus dorsi muscle quantified 2947 proteins, revealing only subtle expression differences between MSTN-mutated and wild-type cattle. These results demonstrate stable inheritance and confirm that MSTN editing does not disrupt genome integrity or protein expression. Overall, our findings support the safety and utility of MSTN gene editing to improve livestock productivity for future food security.

Animals↗

PubMed2000

The evolution of an alpha-esterase pseudogene inactivated in the Drosophila melanogaster lineage.

Previous analyses of the alpha-esterase cluster of Drosophila melanogaster revealed 10 active genes and the DmalphaE4a-Psi pseudogene. Here, we reconstruct the evolution of the pseudogene from the sequences of 12 alleles from widely scattered D. melanogaster populations and single alleles from Drosophila simulans and Drosophila yakuba. All of the DmalphaE4a-Psi alleles contain numerous inactivating mutations, suggesting that pseudogene alleles are fixed in natural populations. Several lines of evidence also suggest that DmalphaE4a is now evolving without selective constraint in the D. melanogaster lineage. There are three polymorphic indels which result in frameshifts; a key nucleotide of the intron splice acceptor is polymorphic; the neutral mutation parameter is the same for replacement and silent sites; one of the nonsilent polymorphisms results in a stop codon; only 1 of the 13 replacement polymorphisms is biochemically conservative; residues that are conserved among active esterases have different states in DmalphaE4a-Psi; and there are about half as many transitional polymorphisms as transversional ones. In contrast, the D. simulans and D. yakuba orthologs DsalphaE4a and DyalphaE4a do not have the inactivating mutations of DmalphaE4a-Psi and appear to be evolving under the purifying selection typical of protein- encoding genes. For instance, there have been more substitutions in the introns than in the exons, and more in silent sites than in replacement sites. Furthermore, most of the amino acid substitutions that have occurred between DyalphaE4a and DsalphaE4a are located in sites that typically vary among active alpha-esterases rather than those that are usually conserved. We argue that the original alphaE4a gene had a function which it has lost since the divergence of the D. melanogaster and D. simulans lineages.

Amino Acid Sequence↗

PubMed1993

Empirical and structural models for insertions and deletions in the divergent evolution of proteins.

The exhaustive matching of the protein sequence database makes possible a broadly based study of insertions and deletions (indels) during divergent evolution. In this study, the probability of a gap in an alignment of a pair of homologous protein sequences was found to increase with the evolutionary distance measured in PAM units (number of accepted point mutations per 100 amino acid residues). A relationship between the average number of amino acid residues between indels and evolutionary distance suggests that a unit 30 to 40 amino acid residues in length remains, on average, undisrupted by indels during divergent evolution. Further, the probability of a gap was found to be inversely proportional to gap length raised to the 1.7 power. This empirical law fits closely over the entire range of gap lengths examined. Gap length distribution is largely independent of evolutionary distance. These results rule out the widely used linear gap penalty as a satisfactory formula for scoring gaps when constructing alignments. Further, the observed gap length distribution can be explained by a simple model of selective pressures governing the acceptance of indels during divergent evolution. Finally, this model provides theoretical support for using indels as part of "parsing algorithms", important in the de novo prediction of the folded structure of proteins from the sequence data.

Amino Acid Sequence↗

PubMed2001

Genomic gigantism: DNA loss is slow in mountain grasshoppers.

Several studies have shown DNA loss to be inversely correlated with genome size in animals. These studies include a comparison between Drosophila and the cricket, Laupala, but there has been no assessment of DNA loss in insects with very large genomes. Podisma pedestris, the brown mountain grasshopper, has a genome over 100 times as large as that of Drosophila and 10 times as large as that of Laupala. We used 58 paralogous nuclear pseudogenes of mitochondrial origin to study the characteristics of insertion, deletion, and point substitution in P. pedestris and Italopodisma. In animals, these pseudogenes are "dead on arrival"; they are abundant in many different eukaryotes, and their mitochondrial origin simplifies the identification of point substitutions accumulated in nuclear pseudogene lineages. There appears to be a mononucleotide repeat within the 643-bp pseudogene sequence studied that acts as a strong hot spot for insertions or deletions (indels). Because the data for other insect species did not contain such an unusual region, hot spots were excluded from species comparisons. The rate of DNA loss relative to point substitution appears to be considerably and significantly lower in the grasshoppers studied than in Drosophila or Laupala. This suggests that the inverse correlation between genome size and the rate of DNA loss can be extended to comparisons between insects with large or gigantic genomes (i.e., Laupala and Podisma). The low rate of DNA loss implies that in grasshoppers, the accumulation of point mutations is a more potent force for obscuring ancient pseudogenes than their loss through indel accumulation, whereas the reverse is true for Drosophila. The main factor contributing to the difference in the rates of DNA loss estimated for grasshoppers, crickets, and Drosophila appears to be deletion size. Large deletions are relatively rare in Podisma and Italopodisma.

Animals↗

PubMed2004

Causes of size homoplasy among chloroplast microsatellites in closely related Clusia species.

Chloroplast DNA sequences and microsatellites are useful tools for phylogenetic as well as population genetic analyses of plants. Chloroplast microsatellites tend to be less variable than nuclear microsatellites and therefore they may not be as powerful as nuclear microsatellites for within-species population analysis. However, chloroplast microsatellites may be useful for phylogenetic analysis between closely related taxa when more conventional loci, such as ITS or chloroplast sequence data, are not variable enough to resolve phylogenetic relationships in all clades. To determine the limits of chloroplast microsatellites as tools in phylogenetic analyses, we need to understand their evolution. Thus, we examined and compared phylogenetic relationships of species within the genus Clusia, using both chloroplast sequence data and variation at seven chloroplast microsatellite loci. Neither ITS nor chloroplast sequences were variable enough to resolve relationships within some sections of the genus, yet chloroplast microsatellite loci were too variable to provide any useful phylogenetic information. Size homoplasy was apparent, caused by base substitutions within the microsatellite, base substitutions in the flanking regions, indels in the flanking regions, multiple microsatellites within a fragment, and forward/reverse mutations of repeat length resulting in microsatellites of identical base composition that were not identical by descent.

Base Sequence↗

PubMed2025

Mutagenic Impact and Evolutionary Influence of Chemoradiotherapy in Hematologic Malignancies.

UNLABELLED: Ionizing radiotherapy (RT) is a widely used treatment strategy for malignancies. In solid tumors, RT-induced double-strand breaks lead to the accumulation of insertion-deletions (indels; ID), and their repair by nonhomologous end joining has been linked to the ID8 mutational signature in surviving cells. However, the extent of RT-induced mutagenesis in hematologic malignancies and its impact on their mutational profiles and interplay with commonly used chemotherapies has not yet been explored. In this study, we interrogated 580 whole-genome sequence (WGS) samples from patients with large B-cell lymphoma, multiple myeloma, and myeloid neoplasms and identified ID8 only in relapsed disease. Yet ID8 was detected after exposure to both RT and mutagenic chemotherapy (i.e., platinum and melphalan). Using WGS of single-cell colonies derived from treated lymphoma cells, we revealed a dose-response relationship between RT and platinum and ID8. Finally, using ID8 as a genomic barcode, we demonstrate that a single RT-surviving cell may seed distant relapse. SIGNIFICANCE: RT and the ID8 indel signature are related, but their genomic impact on hematologic malignancies is unclear. Leveraging WGS, we linked ID8 to both RT and mutagenic chemotherapy and validated that platinum can induce ID8. We used ID8 as a genomic barcode to reveal that RT-resistant cells may seed systemic relapse.

Humans↗

PubMed1998

Evolution of a mammalian promoter through changes in patterns of transcription factor binding.

Previously, we identified evolutionarily derived changes in the pattern of nuclear factor binding to overlapping sites (termed A, B, and C, in proximal to distal order) within the promoter region of the murine D7Rp2e gene. In Mus domesticus, strong binding occurs at site B, with weak binding at sites A and C; in M. pahari, strong binding is observed at sites A and C, with weak binding at site B. The M. pahari-specific occupancy pattern is associated with repression of transcription. Presently, we have identified two interspecies changes, a single A --> G substitution within site B and a 7-bp indel within site C, that are primary determinants of the interspecies differences. These mutations alter both the pattern of binding site occupancy and the ability to repress transcription. The results are discussed in terms of a model involving a dynamic equilibrium between states of binding site occupancy, the nature of which can be modified during evolution.

Animals↗

PubMed2004

The development of ISSR-derived SCAR markers around the SEASONAL FLOWERING LOCUS (SFL) in Fragaria vesca.

Fragaria vesca is a short-lived perennial with a seasonal-flowering habit. Seasonality of flowering is widespread in the Rosaceae and is also found in the majority of temperate polycarpic perennials. Genetic analysis has shown that seasonal flowering is controlled by a single gene in F. vesca, the SEASONAL FLOWERING LOCUS ( SFL). Here, we report progress towards the marker-assisted selection and positional cloning of SFL, in which three ISSR markers linked to SFL were converted to locus-specific sequence-characterized amplified region (SCAR1-SCAR3) markers to allow large-scale screening of mapping progenies. We believe this is the first study describing the development of SCAR markers from ISSR profiles. The work also provides useful insight into the nature of polymorphisms generated by the ISSR marker system. Our results indicate that the ISSR polymorphisms originally detected were probably caused by point mutations in the positions targeted by primer anchors (causing differential PCR failure), by indels within the amplicon (leading to variation in amplicon size) and by internal sequence differences (leading to variation in DNA folding and so in band mobility). The cause of the original ISSR polymorphism was important in the selection of appropriate strategies for SCAR-marker development. The SCAR markers produced were mapped using a F. vesca f. vesca x F. vesca f. semperflorens testcross population. Marker SCAR2 was inseparable from the SFL, whereas SCAR1 mapped 3.0 cM to the north of the gene and SCAR3 1.7 cM to its south.

Base Sequence↗

PubMed2005

Sequence turnover and tandem repeats in cis-regulatory modules in drosophila.

The path by which regulatory sequence can change, yet preserve function, is an important open question for both evolution and bioinformatics. The recent sequencing of two additional species of Drosophila plus the wealth of data on gene regulation in the fruit fly provides new means for addressing this question. For regulatory sequences, indels account for more base pairs (bp) of change than substitutions (between Drosophila melanogaster and Drosophila yakuba), though they are fewer in number. Using Drosophila pseudoobscura as an out-group, we can distinguish insertions from deletions (with maximum parsimony criteria), and find a ratio between 1 and 5 (insertions to deletions) that is species dependent and much larger than the ratio of 1/8 for neutral sequences (Petrov and Hartl 1998). Because neutral sequence is rapidly cleared from the genome, most noncoding regions which preserve their length between D. melanogaster-D. pseudoobscura and have an excess of insertions over deletions should be functional. A fraction of 15%-18% (i.e., more than 20 standard deviations from random expectation) of the regulatory sequence is covered by low copy number tandem repeats whose repeating unit has an average length of 5-10 bp and which occur preferentially (25%-45% coverage) in indels. All indels may be due to tandem repeats if we extrapolate the detection efficiency of the repeat-finding algorithms using the observed point mutation rate between the species we compare. Sequence creation by local duplication accords with the tendency for multiple copies of transcription factor-binding sites to occur in regulatory modules. Thus, indel events and tandem repeats in particular need to be incorporated into models of regulatory evolution because they can alter the rate at which beneficial variants arise and should also influence bioinformatic algorithms that parse regulatory sequences into binding sites.

Animals↗

PubMed1997

Molecular evolution of the second ancient human mariner transposon, Hsmar2, illustrates patterns of neutral evolution in the human genome lineage.

A consensus sequence for the second ancient mariner identified in the human genome, Hsmar2, was constructed by majority rule from full-length and partial sequences of 44 of the +/-1000 copies in the genome. This 1300 base pair (bp) consensus has 31 bp imperfect terminal repeats (ITRs) and encodes a 351 amino acid (aa) mariner transposase. The sequence of this transposase has allowed classification of Hsmar2 as a basal lineage of the irritans subfamily of mariners, sharing at most 38% aa identity with other members of the subfamily. The individual copies in the human genome are all highly mutated from the consensus, having suffered numerous small and some large insertions and deletions (indels), including many insertions of S and J subfamily Alu elements. The copies differ, on average, from the consensus by 11.6%, have suffered 11.8 indels per kilobase (kb), and only 3.7% of the 30 hypermutable CpG dinucleotide pairs in the consensus remain intact. This level of divergence indicates that the ancestrally active Hsmar2 element represented by the consensus was present in the human genome lineage about 80 million years (Myr) ago. Each copy has apparently evolved since then largely independently of the others, and with little constraint on its transposase coding capacity. This pattern of molecular evolution fits the current model for mariner transposon evolution. These copies provide multiple independent datasets for evaluating the pattern of neutral evolution in the human genome, for example, they confirm that most indels are very short and that deletions are twice as common as insertions.

Amino Acid Sequence↗

PubMed2004

Molecular evolution and phylogenetic utility of the petD group II intron: a case study in basal angiosperms.

Sequences of spacers and group I introns in plant chloroplast genomes have recently been shown to be very effective in phylogenetic reconstruction at higher taxonomic levels and not only for inferring relationships among species. Group II introns, being more frequent in those genomes than group I introns, may be further promising markers. Because group II introns are structurally constrained, we assumed that sequences of a group II intron should be alignable across seed plants. We designed universal amplification primers for the petD intron and sequenced this intron in a representative selection of 47 angiosperms and three gymnosperms. Our sampling of taxa is the most representative of major seed plant lineages to date for group II introns. Through differential analysis of structural partitions, we studied patterns of molecular evolution and their contribution to phylogenetic signal. Nonpairing stretches (loops, bulges, and interhelical nucleotides) were considerably more variable in both substitutions and indels than in helical elements. Differences among the domains are basically a function of their structural composition. After the exclusion of four mutational hotspots accounting for less than 18% of sequence length, which are located in loops of domains I and IV, all sequences could be aligned unambiguously across seed plants. Microstructural changes predominantly occurred in loop regions and are mostly simple sequence repeats. An indel matrix comprising 241 characters revealed microstructural changes to be of lower homoplasy than are substitutions. In showing Amborella first branching and providing support for a magnoliid clade through a synapomorphic indel, the petD data set proved effective in testing between alternative hypotheses on the basal nodes of the angiosperm tree. Within angiosperms, group II introns offer phylogenetic signal that is intermediate in information content between that of spacers and group I introns on the one hand and coding sequences on the other.

Base Sequence↗

PubMed2026

Whole genome and exome sequencing of pancreatic neuroendocrine tumour to investigate PRRT response.

Patients with pancreatic neuroendocrine tumours (PNETs) often have similar baseline clinical characteristics, including grade and molecular imaging phenotype, yet have highly variable responses to peptide receptor radionuclide therapy (PRRT). To identify genomic alterations and mutational patterns associated with PRRT treatment response and acquired somatic changes following PRRT exposure, whole genome or exome sequencing was applied to 40 PNET samples from 32 patients, including eight paired pre- or post-PRRT samples. The genomic profile of tumours reflected the known mutational landscape of PNET with MEN1 (34%), ATRX/DAXX (47%) alterations and a recurrent pattern of aneuploidy (38%) detected. A recurrent PSIP1::TBL1X fusion of unknown function was also identified in four tumours. The disease control rate following PRRT using RECIST1.1 and molecular imaging criteria was 88% (28/32). No mutational features were found to be statistically associated with progression-free survival. There was no significant increase in tumour mutational burden in the post-PRRT tumours, nor recurrent emergent mutational changes in cancer driver genes to explain progression to higher-grade disease, when observed. However, a small indel signature (ID8) previously associated with DNA damage repair by non-homologous end joining (NHEJ) was higher in PRRT-exposed compared with PRRT-naive samples (23.8 vs 4.8%, respectively; P < 0.001). Thus, comprehensive DNA analysis of pancreatic NETs did not identify biomarkers predictive of PRRT response nor evidence for high-level PRRT-induced genomic instability or hypermutation, yet mutation signature analysis supports NHEJ as being important for DNA repair and survival of neuroendocrine cells following exposure to beta-particle radiation.

Humans↗

PubMed1993

Location of tolerated insertions/deletions in the structure of the maltose binding protein.

In a previous study [(1987) J. Mol. Biol. 194, 663-673], we isolated ten insertion/deletion mutants (indels) of the maltose binding protein for which the maltose binding constant was only a little or not at all affected. In this paper, we have localized these mutations in the recently solved three-dimensional structure. Contrary to the general expectation, most of the insertion/deletion modifications occurred within elements of secondary structure. An analysis of the inserted residues for three indels found within alpha helices allowed an interpretation regarding protein structure accommodation to such modifications.

ATP-Binding Cassette Transporters↗

PubMed2026

fRagmentomics: an R package for integrating cell-free DNA fragment features with mutational status to support liquid biopsy interpretation.

SUMMARY: Liquid biopsy offers a non-invasive approach to study tumor-derived genetic material circulating in plasma. Beyond genetic alterations, the fragmentomic features of cell-free DNA-such as fragment size, genomic position, and end-motifs-provide valuable insights into the biological and clinical context of DNA release. fRagmentomics is a user-friendly R package designed to characterize cfDNA fragments overlapping one or multiple small mutations of any type, starting from an aligned sequencing file (BAM). It supports multiple mutation input formats, accommodates one-based and zero-based genomic conventions, resolves mutation representation ambiguities, and accepts any reference file in FASTA format. For each fragment overlapping a mutation of interest, fRagmentomics outputs fragment-level features including its fragment size, end-motifs, and mutational status, along with additional fragment-level or read-level information. The package implements an indel-aware and optionally soft-clip-preserving fragment size computation that improves accuracy over conventional size estimates based solely on aligned positions. AVAILABILITY AND IMPLEMENTATION: fRagmentomics is licensed under GNU General Public License v3.0 and available at https://github.com/ElsaB-Lab/fRagmentomics, https://anaconda.org/elsab-lab/r-fragmentomics and https://bioconductor.org/packages/fRagmentomics, with documentation and a tutorial. CONTACT: yoann.pradat@gustaveroussy.fr, elsa.bernard@gustaveroussy.fr. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Software↗