Hundreds Of Nucleotides Research Articles

Unidirectional gene pairs in archaea and bacteria require overlaps or very short intergenic distances for translational coupling via termination-reinitiation and often encode subunits of heteromeric complexes.

Genomes of bacteria and archaea contain a much larger fraction of unidirectional (serial) gene pairs than convergent or divergent gene pairs. Many of the unidirectional gene pairs have short overlaps of -4 nt and -1 nt. As shown previously, translation of the genes in overlapping unidirectional gene pairs is tightly coupled. Two alternative models for the fate of the post-termination ribosome predict either that overlaps or very short intergenic distances are essential for translational coupling or that the undissociated post-termination ribosome can scan through long intergenic regions, up to hundreds of nucleotides. We aimed to experimentally resolve the contradiction between the two models by analyzing three native gene pairs from the model archaeon Haloferax volcanii and three native pairs from Escherichia coli. A two reporter gene system was used to quantify the reinitiation frequency, and several stop codons in the upstream gene were introduced to increase the intergenic distances. For all six gene pairs from two species, an extremely strong dependence of the reinitiation efficiency on the intergenic distance was unequivocally demonstrated, such that even short intergenic distances of about 20 nt almost completely abolished translational coupling. Bioinformatic analysis of the intergenic distances in all unidirectional gene pairs in the genomes of H. volcanii and E. coli and in 1,695 prokaryotic species representative of 49 phyla showed that intergenic distances of -4 nt or -1 nt (= short gene overlaps of 4 nt or 1 nt) were by far most common in all these groups of archaea and bacteria. A small set of genes in E. coli, but not in H. volcanii, had intergenic distances of around +10 nt. Our experimental and bioinformatic analyses clearly show that translational coupling requires short gene overlaps, whereas scanning of intergenic regions by the post-termination ribosome occurs rarely, if at all. Short overlaps are enriched among genes that encode subunits of heteromeric complexes, and co-translational complex formation requiring precise subunit stoichiometry likely confers an evolutionary advantage that drove the formation and conservation of overlapping gene pairs during evolution.

Multiple Non-Canonical Base-Stacking Interactions as One of the Major Determinants of RNA Tertiary Structure Organization.

Stacking interactions of heterocyclic bases of ribonucleotides are one of the most important factors in the organization of RNA secondary and tertiary structure. Most of these (canonical) interactions are formed between adjacent residues in RNA polynucleotide chains. However, with the accumulation of data on the atomic tertiary structures of various RNAs and their complexes with proteins, it has become clear that nucleotide residues that are not adjacent in the polynucleotide chains and are sometimes separated in the RNA primary structure by tens or hundreds of nucleotides can interact via (non-canonical) base stacking. This paper presents an exhaustive database of such nonadjacent base-stacking elements (NA-BSEs) and their environment in the macromolecules of natural and synthetic RNAs. Analysis of these data showed that NA-BSE-forming nucleotides, on average, account for about a quarter of all nucleotides in a particular RNA and, therefore, should be considered as bona fide motifs of the RNA tertiary structure. We also classified NA-BSEs by their location in RNA macromolecules. It was shown that the structure-forming role of NA-BSEs involves compact folding of single-stranded RNA loops, transformation of double-stranded bulges into imperfect helices, and binding of RNA regions distant in the primary and secondary RNA structure.

Cyclic Dinucleotide-Based Enantioselective Fluorination in Water.

The diverse structures of DNA serve as potent chiral scaffolds for DNA-based asymmetric catalysis, yet in most cases tens to hundreds of nucleotides in DNA hybrid catalysts hinder the deep insight into their structure-activity relationship. Owing to the structural simplicity and design flexibility of nucleotides, nucleotide-based catalysts have been emerging as a promising way to obtain fine structural information and understand the catalytic mechanisms. Herein, we found that a cyclic dinucleotide of cyclic di-AMP (c-di-AMP) and 1,10-phenanthroline copper(II) nitrate (Cu(phen)(NO3)2) are assembled to a c-di-AMP-based catalyst (c-di-AMP/Cu(phen)(NO3)2), which could fast achieve enantioselective fluorination in water with 90-99% yields and up to 90% enantiomeric excess (ee). The host-guest interaction between c-di-AMP and Cu(phen)(NO3)2 has been proposed mainly in a supramolecular interaction mode as evidenced by spectroscopic techniques of ultraviolet-visible, fluorescence, circular dichroism, and nuclear magnetic resonance. Cu(phen)(NO3)2 tightly binds to c-di-AMP with a binding constant of 1.7 ± 0.3 × 105 M-1, and the assembly of c-di-AMP/Cu(phen)(NO3)2 shows a modest rate enhancement to carbon-fluorine bond formations as supported by kinetic studies.

G-quadruplex DNA: A Longer Story.

G-quadruplexes (G4s) are distinctive four-stranded DNA or RNA structures found within cells that are thought to play functional roles in gene regulation and transcription, translation, recombination, and DNA damage/repair. While G4 structures can be uni-, bi-, or tetramolecular with respect to strands, folded unimolecular conformations are most significant in vivo. Unimolecular G4 can potentially form in sequences with runs of guanines interspersed with what will become loops in the folded structure: 5'GxLyGxLyGxLyGx, where x is typically 2-4 and y is highly variable. Such sequences are highly conserved and specifically located in genomes. In the folded structure, guanines from each run combine to form planar tetrads with four hydrogen-bonded guanine bases; these tetrads stack on one another to produce four strand segments aligned in specific parallel or antiparallel orientations, connected by the loop sequences. Three types of loops (lateral, diagonal, or "propeller") have been identified. The stacked tetrads form a central cavity that features strong coordination sites for monovalent cations that stabilize the G4 structure, with potassium or sodium preferred. A single monomeric G4 typically forms from a sequence containing roughly 20-30 nucleotides. Such short sequences have been the primary focus of X-ray crystallographic or NMR studies that have produced high-resolution structures of a variety of monomeric G4 conformations. These structures are often used as the basis for drug design efforts to modulate G4 function.We believe that the focus on monomeric G4 structures formed by such short sequences is perhaps myopic. Such short sequences for structural studies are often arbitrarily selected and removed from their native genomic sequence context, and then are often changed from their native sequences by base substitutions or deletions intended to optimize the formation of a homogeneous G4 conformation. We believe instead that G-quadruplexes prefer company and that in a longer natural sequence context multiple adjacent G4 units can form to combine into more complex multimeric G4 structures with richer topographies than simple monomeric forms. Bioinformatic searches of the human genome show that longer sequences with the potential for forming multiple G4 units are common. Telomeric DNA, for example, has a single-stranded overhang of hundreds of nucleotides with the requisite repetitive sequence with the potential for formation of multiple G4s. Numerous extended promoter sequences have similar potentials for multimeric G4 formation. X-ray crystallography and NMR methods are challenged by these longer sequences (>30 nt), so other tools are needed to explore the possible multimeric G4 landscape. We have implemented an integrated structural biology approach to address this challenge. This approach integrates experimental biophysical results with atomic-level molecular modeling and molecular dynamics simulations that provide quantitatively testable model structures. In every long sequence we have studied so far, we found that multimeric G4 structures readily form, with a surprising diversity of structures dependent on the exact native sequence used. In some cases, stable hairpin duplexes form along with G4 units to provide an even richer landscape. This Account provides an overview of our approach and recent progress and provides a new perspective on the G-quadruplex folding landscape.

The Processivity of Telomerase: Insights from Kinetic Simulations and Analyses.

Telomerases are moderately processive reverse transcriptases that use an integral RNA template to extend the 3′ end of linear chromosomes. Processivity values, defined as the probability of extension rather than dissociation, range from about 0.7 to 0.99 at each step. Consequently, an average of tens to hundreds of nucleotides are incorporated before the single-stranded sDNA product dissociates. The RNA template includes a six nucleotide repeat, which must be reset in the active site via a series of translocation steps. Nucleotide addition associated with a translocation event shows a lower processivity (repeat addition processivity, RAP) than that at other positions (nucleotide addition processivity, NAP), giving rise to a characteristic strong band every 6th position when the product DNA is analyzed by gel electrophoresis. Here, we simulate basic reaction mechanisms and analyze the product concentrations using several standard procedures to show how the latter can give rise to systematic errors in the processivity estimate. Complete kinetic analysis of the time course of DNA product concentrations following a chase with excess unlabeled DNA primer (i.e., a pulse-chase experiment) provides the most rigorous approach. This analysis reveals that the higher product concentrations associated with RAP arise from a stalling of nucleotide incorporation reaction during translocation rather than an increased rate constant for the dissociation of DNA from the telomerase.

New Insights into the Importance of Long Non-Coding RNAs in Lung Cancer: Future Clinical Approaches.

In mammals, a large part of the gene expression products come from the non-coding ribonucleotide sequences of the protein. These short and long sequences are within the range of tens to hundreds of nucleotides, encompassing more than 200 RNA molecules, and their function is known as the molecular structure of long non-coding RNA (lncRNA). LncRNA molecules are unique nucleotides that have a substantial role in epigenetic regulation, transcription, and post-transcriptional modifications in different ways. According to the results of recent studies, lncRNAs have been shown to assume various roles, including tumor suppression or oncogenic functions in common types of cancer such as lung and breast cancer. These non-coding RNAs (ncRNAs) play a pivotal role in activating transcription factors, managing the ribonucleoproteins, the framework for collecting co-proteins, intermittent processing regulations, chromatin status alterations, and maintaining the control within the cell. Cutting-edge technologies have been introduced to disclose several types of lncRNAs within the nucleus and the cytoplasm, which have accomplished important achievements that are applicable in medicine. Due to these efforts, various data centers have been created to facilitate and modify scientific information related to these molecules, including detection, classification, biological evolution, gene status, spatial structure, status, and location of these small molecules. In the present study, we attempt to present the impacts of these ncRNAs on lung cancer with an emphasis on their mechanisms and functions.

Finding recurrent RNA structural networks with fast maximal common subgraphs of edge-colored graphs.

RNA tertiary structure is crucial to its many non-coding molecular functions. RNA architecture is shaped by its secondary structure composed of stems, stacked canonical base pairs, enclosing loops. While stems are precisely captured by free-energy models, loops composed of non-canonical base pairs are not. Nor are distant interactions linking together those secondary structure elements (SSEs). Databases of conserved 3D geometries (a.k.a. modules) not captured by energetic models are leveraged for structure prediction and design, but the computational complexity has limited their study to local elements, loops. Representing the RNA structure as a graph has recently allowed to expend this work to pairs of SSEs, uncovering a hierarchical organization of these 3D modules, at great computational cost. Systematically capturing recurrent patterns on a large scale is a main challenge in the study of RNA structures. In this paper, we present an efficient algorithm to compute maximal isomorphisms in edge colored graphs. We extend this algorithm to a framework well suited to identify RNA modules, and fast enough to considerably generalize previous approaches. To exhibit the versatility of our framework, we first reproduce results identifying all common modules spanning more than 2 SSEs, in a few hours instead of weeks. The efficiency of our new algorithm is demonstrated by computing the maximal modules between any pair of entire RNA in the non-redundant corpus of known RNA 3D structures. We observe that the biggest modules our method uncovers compose large shared sub-structure spanning hundreds of nucleotides and base pairs between the ribosomes of Thermus thermophilus, Escherichia Coli, and Pseudomonas aeruginosa.

FLT3-ITD Allele Frequency Is an Independent Prognostic Factor for Poor Outcome in FLT3-ITD–Positive AML Patients

FMS-like tyrosine kinase 3 internal tandem duplication (FLT3-ITD) is a molecular genetic alteration significantly affecting the clinical outcome in patients with acute myeloid leukemia (AML). FLT3-ITD mutations are characterized by variable mutant-to-wild allelic ratios (ARs) and sizes of the duplicated sequences. The size of the inserted sequence may vary from a few to hundreds of nucleotides. The aim of this work was to determine the impact of FLT3-ITD ARs, FLT3-ITD allelic frequency (AF), and allele size in de novo AML. We studied 117 patients with FLT3-ITD gene mutation-positive AML, dividing them into those with low ARs and those with high ARs (>0.64) and examined their prognostic impact. High FLT3-ITD AR ≥ 0.64 and AF ≥ 0.5 were significantly associated with a lower overall survival compared with lower AR (median 0.625 vs. 1.020 months, respectively; P=.041) and AF (median 0.493 vs. 0.954 months, respectively; P=.009). NPM1 mutation had no favorable impact on the low-level FLT3-ITD group. Initial high total leukocyte count, FLT3-ITD AF, and splenomegaly are independent prognostic factors for poor outcome in FLT3-ITD-positive AML.

De novo adult acute myeloid leukemia with two new mutations in juxtatransmembrane domain of the FLT3 gene: a case report

BackgroundApproximately 30% of adult acute myeloid leukemia (AML) acquire within fms-like tyrosine kinase 3 gene (FLT3) internal tandem duplications (FLT3/ITDs) in their juxtamembrane domain (JMD). FLT3/ITDs range in size from three to hundreds of nucleotides, and confer an adverse prognosis. Studies on a possible relationship between of FLT3/ITDs length and clinical outcomes in those AML patients were inconclusive, yet.Case presentationHere we report a 54-year-old Arab male diagnosed with AML who had two FLT3-ITD mutations in addition to NPM1 mutation. Cytogenetic approaches (banding cytogenetics) and fluorescence in situ hybridization (FISH) using specific probes to detect translocations t(8;21), t(15;17), t(16;16), t(12;21), and deletion del(13q)) were applied to exclude chromosomal abnormalities. Molecular genetic approaches (polymerase chain reaction (PCR) and the Sanger sequencing) identified a yet unreported combination of two new mutations in FLT3-ITDs. The first mutation induced a frameshift in JMD, and the second led to a homozygous substitution of c.1836T>A (p.F612L) also in JMD. Additionally a NPM1 type A mutation was detected. The first chemotherapeutic treatment was successful, but 1 month after the initial diagnosis, the patient experienced a relapse and unfortunately died.ConclusionsTo the best of our knowledge, a combination of two FLT3-ITD mutations in JMD together with an NPM1 type A mutation were not previously reported in adult AML. Further studies are necessary to prove or rule out whether the size of these FLT3-ITDs mutations and potential other double mutations in FLT3-ITD are correlated with the observed adverse outcome.

A Spatially Explicit Model of Stabilizing Selection for Improving Phylogenetic Inference

Ultraconserved elements (UCEs) are stretches of hundreds of nucleotides with highly conserved cores flanked by variable regions. Although the selective forces responsible for the preservation of UCEs are unknown, they are nonetheless believed to contain phylogenetically meaningful information from deep to shallow divergence events. Phylogenetic applications of UCEs assume the same degree of rate heterogeneity applies across the entire locus, including variable flanking regions. We present a Wright–Fisher model of selection on nucleotides (SelON) which includes the effects of mutation, drift, and spatially varying, stabilizing selection for an optimal nucleotide sequence. The SelON model assumes the strength of stabilizing selection follows a position-dependent Gaussian function whose exact shape can vary between UCEs. We evaluate SelON by comparing its performance to a simpler and spatially invariant GTR+ model using an empirical data set of 400 vertebrate UCEs used to determine the phylogenetic position of turtles. We observe much improvement in model fit of SelON over the GTR+ model, and support for turtles as sister to lepidosaurs. Overall, the UCE-specific parameters SelON estimates provide a compact way of quantifying the strength and variation in selection within and across UCEs. SelON can also be extended to include more realistic mapping functions between sequence and stabilizing selection as well as allow for greater levels of rate heterogeneity. By more explicitly modeling the nature of selection on UCEs, SelON and similar approaches can be used to better understand the biological mechanisms responsible for their preservation across highly divergent taxa and long evolutionary time scales.

Interaction Networks of Ribosomal Expansion Segments in Kinetoplastids.

Expansion segments (ES) are insertions of a few to hundreds of nucleotides at discrete locations on eukaryotic ribosomal RNA (rRNA) chains. Some cluster around 'hot spots' involved in translation regulation and some may participate in biogenesis. Whether ES play the same roles in different organisms is currently unclear, especially since their size may vary dramatically from one species to another and very little is known about their functions. Most likely, ES variation is linked to adaptation to a particular environment. In this chapter, we compare the interaction networks of ES from four kinetoplastid parasites, which have evolved in diverse insect vectors and mammalian hosts: Trypanosoma cruzi, Trypanosoma brucei, Leishmania donovani and Leishmania major. Here, we comparatively analyze ribosome structures from these representative kinetoplastids and ascertain meaningful structural differences from mammalian ribosomes. We base our analysis on sequence alignments and three-dimensional structures of 80S ribosomes solved by cryo-electron microscopy (cryo-EM). Striking differences in size are observed between ribosomes of different parasites, indicating that not all ES are expanded equally. Larger ES are not always matched by large surrounding ES or proteins extensions in their vicinity, a particularity that may lead to clues about their biological function. ES display different species-specific patterns of conservation, which underscore the density of their interaction network at the surface of the ribosome. Making sense of the conservation patterns of ES is part of a global effort to lay the basis for functional studies aimed at discovering unique kinetoplastid-specific sites suitable for therapeutic applications against these human and often animal pathogens.

Genomic diversity and evolution, diagnosis, prevention, and therapeutics of the pandemic COVID-19 disease.

The coronavirus disease 19 (COVID-19) is a highly transmittable and pathogenic viral infection caused by a novel evolutionarily divergent RNA virus, the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The virus first emerged in Wuhan, China in December 2019, and subsequently spreaded around the world. Genomic analyses revealed that this zoonotic virus may be evolved naturally but not a purposefully manipulated laboratory construct. However, currently available data are not sufficient to precisely conclude the origin of this fearsome virus. Comprehensive annotations of the whole-genomes revealed hundreds of nucleotides, and amino acids mutations, substitutions and/or deletions at different positions of the ever changing SARS-CoV-2 genome. The spike (S) glycoprotein of SARS-CoV-2 possesses a functional polybasic (furin) cleavage site at the S1-S2 boundary through the insertion of 12 nucleotides. It leads to the predicted acquisition of 3-O-linked glycan around the cleavage site. Although real-time RT-PCR methods targeting specific gene(s) have widely been used to diagnose the COVID-19 patients, however, recently developed more convenient, cheap, rapid, and specific diagnostic tools targeting antigens or CRISPR-Cas-mediated method or a newly developed plug and play method should be available for the resource-poor developing countries. A large number of candidate drugs, vaccines and therapies have shown great promise in early trials, however, these candidates of preventive or therapeutic agents have to pass a long path of trials before being released for the practical application against COVID-19. This review updates current knowledge on origin, genomic evolution, development of the diagnostic tools, and the preventive or therapeutic remedies of the COVID-19. We also discussed the future scopes for research, effective management, and surveillance of the newly emerged COVID-19 disease.

Reassessment of Viroid RNA Cytosine Methylation Status at the Single Nucleotide Level.

Composed of a few hundreds of nucleotides, viroids are infectious, circular, non-protein coding RNAs able to usurp plant cellular enzymes and molecular machineries to replicate and move in their hosts. Several secondary and tertiary RNA structural motifs have been implicated in the viroid infectious cycle, but whether modified nucleotides, such as 5C-methylcytosine (m5C), also play a role has not been deeply investigated so far. Here, the possible existence of m5C in both RNA polarity strands of potato spindle tuber viroid and avocado sunblotch viroid -which are representative members of the nucleus- and chloroplast-replicating viroids, respectively- has been assessed at single nucleotide level. We show that a standard bisulfite protocol efficiently used for identifying m5C in cellular RNAs may generate false positive results in the case of the highly structured viroid RNAs. Applying a bisulfite conversion protocol specifically adapted to RNAs with high secondary structure, no m5C was identified in both polarity strands of both viroids, indicating that this specific nucleotide modification does not likely play a role in viroid biology.

Prognostic Impact of Insertion Site in Acute Myeloid Leukemia (AML) with FLT3 Internal Tandem Duplication: Results from the Ratify Study (Alliance 10603)

Prognostic Impact of Insertion Site in Acute Myeloid Leukemia (AML) with FLT3 Internal Tandem Duplication: Results from the Ratify Study (Alliance 10603)

Stress-Induced Low Complexity RNA Activates Physiological Amyloidogenesis

SUMMARYAmyloid bodies (A-bodies) are inducible membraneless nuclear compartments composed of heterogeneous proteins that adopt an amyloid-like state. A-bodies are seeded by noncoding RNA derived from stimuli-specific loci of the rDNA intergenic spacer (rIGSRNA). This raises the question of how rIGSRNA recruits a large population of diverse proteins to confer A-body identity. Here, we show that long low-complexity dinucleotide repeats operate as the architectural determinants of rIGSRNA. On stimulus, clusters of rIGSRNA with simple cytosine/ uracil (CU) or adenosine/guanine (AG) repeats spanning hundreds of nucleotides accumulate in the nucleolar area. The low-complexity sequences facilitate charge-based interactions with short cationic peptides to produce multiple nucleolar liquid-like foci. Local concentration of proteins with fibrillation propensity in these nucleolar foci induces the formation of an amyloidogenic liquid phase that seeds Abodies. These results demonstrate the physiological importance of low-complexity RNA and repetitive regions of genomethe genome often dismissed as “junk” DNA.

CasPER, a method for directed evolution in genomic contexts using mutagenesis and CRISPR/Cas9

Here we describe a method for robust directed evolution using mutagenesis of large sequence spaces in their genomic contexts. The method employs error-prone PCR and Cas9-mediated genome integration of mutant libraries of large-sized donor variants into single or multiple genomic sites with efficiencies reaching 98–99%. From sequencing of genome integrants, we determined that the mutation frequency along the donor fragments is maintained evenly and successfully integrated into the genomic target loci, indicating that there is no bias of mutational load towards the proximity of the double strand break. To validate the applicability of the method for directed evolution of metabolic gene products we engineered two essential enzymes in the mevalonate pathway of Saccharomyces cerevisiae with selected variants supporting up to 11-fold higher production of isoprenoids. Taken together, our method extends on existing CRISPR technologies by facilitating efficient mutagenesis of hundreds of nucleotides in cognate genomic contexts.

Non-enzymatic recombination of RNA: Ligation in loops

BackgroundWhile the RNA world hypothesis is widely accepted, it is still far from complete: the existence of self-replicating ribozyme, consisting of potentially hundreds of nucleotides, is a core assumption for the majority of RNA world models. The appearance of such long RNA molecules under prebiotic conditions is not self-evident. Recombination seems to be a plausible way of creating RNA diversity, resulting in the appearance of functional RNAs, capable of self-replicating. MethodsWe report here on the study of recombination process modelled with two 96 nts RNA fragments. Detection of recombination products was performed with RT-PCR followed by TA-cloning and Sanger sequencing. ResultsA wide range of recombinant products was detected. We found that (i) the most efficient ligation was observed for RNA species forming bulges or internal loops, with ligation partners located within the loop; (ii) a strong preference was observed for formation of a few types of major products with a large variety of minor products; (iii) ligation could occur with participation of either 2′,3′-cyclophosphate or 5′-ppp; (iv) the presence of key reaction components, i.e. 5′ppp-RNAs, enabled the formation of additional types of product; (v) molecular dynamics simulations of one of the most abundant products suggests that the ligation results in a preferable formation of 2′-5′- rather than 3′-5′-linkages. ConclusionsThe study demonstrates regularities of new RNA molecules formation with non-enzymatic recombination process. General significanceOur findings provide new data supporting the RNA World hypothesis and show the way of new RNA sequences emergence under prebiotic conditions.

Detection of homozygous deletions in tumor-suppressor genes ranging from dozen to hundreds nucleotides in cancer models.

Tumor-suppressor genes can be inactivated by several mechanisms and, in a majority of cases, both alleles need to be affected. One of the mechanisms of inactivation is due to deletions ranging from dozen to hundreds of nucleotides; such deletions are often missed by variant callers. HomDelDetect is a method to detect such homozygous deletions in cancer models, such as cancer cell lines and potentially patient tumor-derived xenografts. This method can be applied to partial exome, whole-exome sequencing, whole-genome sequencing, and RNA-seq data. We applied our method across a panel of CCLE cancer cell lines and observed good concordance with SNP array-based analysis and also detected deletions that have been missed by variant callers and by SNP arrays, demonstrating the ability of HomDelDetect to improve the annotations of tumor-suppressor genes in cancer models.

DNA repair: Clamping down on copy errors.

Repair enzymes must communicate across hundreds of nucleotides to undo errors made during DNA replication. Imaging reveals that the enzymes do this by forming a series of ring-like clamps that diffuse along the DNA.

SVM2Motif--Reconstructing Overlapping DNA Sequence Motifs by Mimicking an SVM Predictor.

Identifying discriminative motifs underlying the functionality and evolution of organisms is a major challenge in computational biology. Machine learning approaches such as support vector machines (SVMs) achieve state-of-the-art performances in genomic discrimination tasks, but—due to its black-box character—motifs underlying its decision function are largely unknown. As a remedy, positional oligomer importance matrices (POIMs) allow us to visualize the significance of position-specific subsequences. Although being a major step towards the explanation of trained SVM models, they suffer from the fact that their size grows exponentially in the length of the motif, which renders their manual inspection feasible only for comparably small motif sizes, typically k ≤ 5. In this work, we extend the work on positional oligomer importance matrices, by presenting a new machine-learning methodology, entitled motifPOIM, to extract the truly relevant motifs—regardless of their length and complexity—underlying the predictions of a trained SVM model. Our framework thereby considers the motifs as free parameters in a probabilistic model, a task which can be phrased as a non-convex optimization problem. The exponential dependence of the POIM size on the oligomer length poses a major numerical challenge, which we address by an efficient optimization framework that allows us to find possibly overlapping motifs consisting of up to hundreds of nucleotides. We demonstrate the efficacy of our approach on a synthetic data set as well as a real-world human splice site data set.