Skip to main content

Greedy de novo motif discovery to construct motif repositories for bacterial proteomes

Abstract

Background

Bacterial surfaces are complex systems, constructed from membranes, peptidoglycan and, importantly, proteins. The proteins play crucial roles as critical regulators of how the bacterium interacts with and survive in its environment. A full catalog of the motifs in protein families and their relative conservation grade is a prerequisite to target the protein-protein interaction that bacterial surface protein makes to host proteins.

Results

In this paper, we propose a greedy approach to identify conserved motifs in large sequence families iteratively. Each iteration discovers a motif de novo and masks all occurrences of that motif. Remaining unmasked sequences are subjected to the next round of motif detection until no more significant motifs can be found. We demonstrate the utility of the method through the construction of a proteome-wide motif repository for Group A Streptococcus (GAS), a significant human pathogen. GAS produce numerous surface proteins that interact with over 100 human plasma proteins, helping the bacteria to evade the host immune response. We used the repository to find that proteins part of the bacterial surface has motif architectures that differ from intracellular proteins.

Conclusions

We elucidate that the M protein, a coiled-coil homodimer that extends over 500 A from the cell wall, has a motif architecture that differs between various GAS strains. As the M protein is known to bind a variety of different plasma proteins, the results indicate that the different motif architectures are responsible for the quantitative differences of plasma proteins that various strains bind. The speed and applicability of the method enable its application to all major human pathogens.

Background

The rise of antibiotics resistant bacteria poses a major global health issue predicted to cause 10 million deaths per year in 2050, more than heart disease and cancer combined [1]. The increasing resistance to antibiotics necessitates the development of alternative treatment strategies. One promising alternative treatment strategy includes the disruption of protein binding interfaces between bacteria and human proteins to disarm bacterial defense systems [2]. Such strategies require high-confident identification of sequence motifs that correspond to a structural unit that are necessary for protein folding or binding of ligands and other proteins.

Motifs are short segments of a protein sequence which shows a level of conservation throughout a protein family and beyond. Conserved motifs can be extracted from multiple sequence alignment of proteins with similar functions in different species. While finding such motifs can provide insights for prediction of functional residues, identifying and understanding them is fundamental to discovering binding interfaces in protein complexes [3]. It is generally believed that the binding interfaces forming interactions to help bacteria evade the immune system or to obtain nutrients are comparatively more conserved compared to interactions that are benefiting the host, such as surface exposed epitope. Over time, this results in segments of exposed proteins that are significantly more conserved for functional reasons.

Disrupting the protein-protein interactions by targeting the conserved segments would potentially facilitate the host immune response [46]. However, the high variability of bacterial surface proteins makes it challenging to study them with traditional sequence analysis methods. InterPro for example [7] contains motifs for the anchor and the signal peptide whereas the rest of the protein sequence remains largely unannotated. Multiple-sequence alignment algorithms typically run into problems with the variable number of repeats and tends to produce highly gapped alignments. The rapid growth of known bacterial protein sequences presents an opportunity to identify protein-family specific motifs (in contrast to Interpro that attempts to find motifs common to multiple families).

Group A streptococcus (GAS) is one of the most important bacterial pathogens causing over 700 million mild infections such as tonsillitis, impetigo and erysipelas and, occasionally, severe invasive infections including sepsis, meningitis or necrotizing fasciitis with mortality rates up to 25% [8]. Surface proteins play important roles in the interaction with host proteins [9]. Several bacterial surface proteins interact with numerous of host proteins, forming complex protein-protein interaction networks.

One of the key surface proteins of S. pyogenes is the M protein, a coiled-coil homodimer that extends over 500 Å from the cell wall. The M protein is capable of binding several plasma proteins such as fibrinogen [6] and albumin [10, 11]. A crystal structure of M and fibrinogen was published in 2011 demonstrates that the M and fibrinogen form a cross-like complex structure. Further, the M protein is composed of several repeats that are present a variable number of times; some of these repeats overlap with protein-protein interactions binding interfaces [1215]. Accordingly, a comprehensive repository of the motifs in coiled-coil proteins and their relative conservation grade is a prerequisite to target the protein-protein interaction that bacterial surface protein makes to host proteins [16].

Here, we present a strategy to iteratively identify protein-family specific motifs from large genome resources, then mask all occurrences of these motifs until no more significant motifs can be found. We applied this strategy to a GAS strain as a model system. We constructed a compendium of almost 60000 motifs for GAS. Further, we demonstrate the power of the approach using the M protein and describe the motif resource in general terms.

Methods

Outline of the algorithm

The algorithm starts with a database of protein sequence families and sub-selects a user-defined number (here 100) of sequences for each family containing more than 100 sequences as shown in the pseudocode in Fig. 1. The main part of the algorithm finds the first motif in a sequence and mask all identified occurrences, then remove them from the sequence and produce one new sequence for each occurrence in the second iteration. After finishing this loop, all identified motifs are stored in the main repository and followed by architectural analysis by considering the occurrences of each motif in the entire genome and by computing the internal overlaps of that motif inside the family as well. In the last step, the results will be compared with InterPro to find common overlapping motifs to report and the new ones for further analysis. All the results as well as the main repository of all discovered motifs are stored and available in a SQLite table. SQLite is an embedded SQL database engine that implements a transactional SQL database engine. The code for SQLite is in the public domain and free for use.

Fig. 1
figure1

The pseudo-algorithm of de novo motif discovery approach

Construction of protein families

We selected a representative genome from an invasive M1 S. pyogenes isolated in Ontario, Canada. This sample is available with id 293653.4 from PatricBRC, the bacterial bioinformatics resource database [17]. This genome has 1931 coding sequences (CDS). We downloaded additional 70459 genomes from PatricBRC. This number here refers to the number of genomes available at PatricBRC that had both an.ffa protein fasta file and.cds files that contains a table which links the PatricBRC sequence accession number to the FigFam ID [18]. We used this resource to build one protein fasta file per FigFam ID filtering out duplicate entries. We constructed 1564 FIGfams families containing a total of 9,041,083 protein sequences of which 3,817,065 were unique at the amino acid level. This sequence resource was used as input to the workflow outlined in Fig. 2.

Fig. 2
figure2

The workflow of de novo motif discovery approach. All FIGfams families were downloaded for a user-specified organism and get through the core of the processing -MEME and FIMO- where all the motifs are discovered and masked in an iterative process. A sqlite repository stores all the motifs and motif architectures with required information like organism name, FIGfams family name, the start and stop points on the sequence and etc

MEME and FIMO

Figure 2 shows the general workflow of our approach, where we make use of MEME [19, 20] and FIMO [21] in the core part of the system to handle motif discovery and masking the multiple occurrences of each motif on the sequence. MEME is an open-source application which has been widely used for sequence motif discovery and analysis in both DNA and proteins. It is based on GLAM2 algorithm [22] and enables covering of motifs containing gaps. While MEME finds a single occurrence of a motif in the sequence, FIMO is able to consider the MEME’s output and define multiple occurrences for any individual gapped or un-gapped motifs. FIMO assign different scores for each matched sequence according to a dynamic programming approach [23] and then motif-specific q-values are computed based on a bootstrap procedure [24]. FIMO’s outputs are considered according to their p-values, and q-values make it possible to set a user-defined thresholds to cover only specific motif occurrences.

InterPro

InterProScan [7] is a reference resource that provides a functional analysis of protein sequences by classifying them into families and predicting the presence of domains and important sites. In order to achieve a general view of the coverage of our approach, we compared the generated de novo based motif repository of GAS with all GAS-related motifs in InterProScan.

Assigning proteins to cellular compartments

All proteins were assigned to one of seven compartments by using information from mass spectrometry experiments, annotations from several databases followed by manual curation. In short, we identified exposed, cell wall associated and secreted protein using data from Karlsson et al. [9]. Transmembrane proteins were identified using TMHMM [25]. DNA associated protein and transcription factors were identified using InterPro [7] and RegPrecise [26]. All other proteins were assigned to the intracellular compartment.

Software availability

MEME and FIMO 4.11.1 was used through out the project. The workflow is implemented in GC3pie [27] which makes it possible to parallelize over all available computational cores. All parts of the workflow are written in Python 2.7 and is wrapped by applicake, an open-source and free framework useful when designing workflows. The workflow is available through a singularity container [28] and the container together with the data and an ipython notebook contains instructions and examples to parse the data are provided online with this DOI: 10.5281/zenodo.1403142.

Results

Sub-selection

We analyzed a large sequence database of all GAS proteins containing 1564 FIGfams sequence families as outlined in the Methods section. The FIGfams contain a different number of sequences. This begs the question whether a subset of them would be sufficient to cover most of the motifs. We designed a general sub-selection test to reduce the number of sequences due to computational resource reasons. The sub-selection considers two different families and select a set of 2, 10, 20, 50, 100, 250, 500, and 1000 sequences randomly and repeat the whole analysis for 10 times. In each sub-selection test, we ran the workflow to find all the motifs, and we made an average of motif-coverage between all 10 repeats. Figure 3 demonstrates that sub-selection of 100 sequences is sufficient to cover the majority number of all motifs while reducing time and computational resources more than 20-fold (see Table 1).

Fig. 3
figure3

Sub-selection test on two sample families. Two sample families are selected and analyzed by sub-selection test. The bubble graph indicates that selection of 100 sequences has a good coverage while saving computational resources more than 20-fold

Table 1 Comparison of sub-selection test between different families with different number of selected sequences

MEME/FIMO

The workflow starts by entering the name of the desired organism and the q-value cut-off (optional) which are the only required inputs (Fig. 2). In the second phase, all FIGfams protein families related to the input organism are downloaded and stored in a database. Then, by considering the accessibility of computational resources, de novo motif discovery on protein families starts. Figure 4 shows two sample runs of the algorithm where MEME is applied to the sequence collection, restricting the number of identified motifs to one. Motif occurrences were discovered in the sequence collection using FIMO, and only occurrences with e-values of 1e-6 or lower were considered. The proteins were split using the number of occurrences and remaining parts longer than ten amino acids are carried forward to create a new merged sequence collection, mixed with full-length and partial proteins.

Fig. 4
figure4

Two rounds (first to the left, second to the right) of the algorithm is displayed. It starts with a collection of sequences and then discovers motifs in this collection using MEME. It uses FIMO to find additional occurrences of the motif within the sequence collection. In the second round, the motifs are masked (gray bars) before MEME is applied once more. The algorithm iterates through round 3 to N until no more motifs are found, or the sequence collection is fully annotated

The new sequence collection is used as the input for each iterative round of MEME, FIMO, split until no more significant motifs could be discovered, or all remaining sub-sequence were below ten amino acids. All the motif occurrences with corresponding features are stored in an SQLite database. To give further information to the user, known motifs are also integrated from [10] and the InterPro database and visualized using pViz.js [29].

Protein M1 Motif discovery

As an example of application on specific protein family, we collected a large sequence collection of M proteins from four sources: PatricBRC, genomes we have previously sequenced and assembled [30], the M database from CDC (Centers for Disease Control and prevention) and the UniProtKB/TrEMBL database. Any M protein sequence without motifs representing an anchor or a signal peptide was discarded, and the remaining sequences were reduced to a 98% sequence identity using CD-HIT [31]. In total, the algorithm ended after 18 rounds, resulting in 20 motifs from the M protein sequence collection. The SF370 M1 protein reference [32] contained motifs m01-m03, m05-m08, m11-m14, m16-m17 but not m04, m09-10, m15, and m18-20. Additional file 1 contains the logo of all discovered motifs.

Figure 5 shows the general motif architecture as the output of the algorithm. Note that an architecture (motif pattern) shows the distribution of motifs over the entire protein family. By considering such representation, it is possible to show the general motif pattern that most of the proteins in the family follow. So, the architectural motif view helps to find potential protein-protein interaction binding sites as the majority member of the family desire to follow such pattern. Accordingly, we found a total of 123 motif architectures, and of these, 85% (104) are associated with a single serotype.

Fig. 5
figure5

The auto-generated result of our approach on M1 protein. a: Binding interfaces of fibrinogen according to the reference crystal structure (PDB id 2XNX). [6]. b: M1 domains proposed in [10]. c: The output results of our approach

Here in M protein, architecture [m01:3, m02:1, m03:1, m05:1, m06:4, m07:3, m08:1, m11:1, m12:1, m14:1, m16:1, m17:1] is the architecture that exists in several serotypes (emm52, emm23, emm16, emm83, emm10). For the M1 proteins, we identified seven motif architectures (Table 2). For emm1 in Fig. 5, we see that m08 and m02 are the first and second part of the YSIRK signal peptide. m03 largely overlap with the anchor region. m01 and m06 correspond to the C repeats, m07 overlaps with the B repeats although we fail to identify the second and third B-domain. m13 finds the S region and m14 overlap partly with the A domain. The D domain is largely split into several motifs - m12, m18 followed by m05.

Table 2 Architectures identified for M1 protein

Analysis of conserved motif in GAS genome

We evaluated identified motifs separately based on protein families in different cellular compartments (Table 3). The main idea is to provide a general comparison between protein families in different cellular compartments in terms of motif-based conservation grade which helps to discover the general evolutionary pressure on cellular compartments and further distinguishing potential drug targets inside and outside the cell. To do that, one should consider the fact that the number of motifs in each compartment is a function of sequence length and the average is dependent on sequence variability. Such dependency affects the comparison between different protein families led to results that are biased against sequence length. To address these issues, we represent motif architecture per sequence and most importantly per family. Accordingly, each protein or its related family can have one or several architectures based on motif variability on that family. Consequently, protein families with few architectures indicates higher sequence conservation inside the family and generally shows that the family has more conserved motifs to do special cellular functions. In this way, by comparing the number of architecture in two different protein families, it is possible to state which family is more conserved.

Table 3 Protein characterization in different cellular compartments

As shown in Table 4, most variable proteins in GAS are transmembrane and secreted proteins which are less conserved and have a more diversified interaction with host proteins. Most conserved proteins are DNA-related and transcription factors together with intracellular proteins that have special machinery roles inside the cell. Transmembrane proteins which play crucial role as the transportation system on the bacterial surface are also more evolved according to the evolutionary pressure. In general as Fig. 6 indicates, we can conclude that the evolutionary pressure is lower on intracellular proteins compared to the surface and secreted proteins.

Fig. 6
figure6

Motif-based architectural comparison between different cellular compartments. The bar-plot to the left shows the general comparison while different cellular compartments are separated in a schematic cell view to the right. Transmembrane proteins (TM) and secreted proteins have the least number of architecture per family. It shows that motif diversity in these compartments are high and changing by the time. In contrary DNA-related proteins and Transcription Factors (TF) show more conserve motifs on their sequence with having the highest number of architecture per family

Table 4 Comprehensive comparison based on number of motifs per sequence and architecture in different cellular compartments

Comparison with InterPro

To compare our results to InterPro, we analyzed and filtered motifs based on their signature from InterProScan which revealed that 11996 distinct motifs related to GAS are not recognized/discovered by InterProScan (71.15% of all discovered motifs) while there are many important motifs also in common (28.85%). Table 5 contains the list of most commonly overlapping motifs with special InterPro description which discovered by our approach.

Table 5 The InterProScan results that are most commonly overlapping with a motif

Discussion

Conserved protein sequence domains, also referred to as motifs, play an important role in protein function, protein structure and protein-protein interactions. Motifs are the results of several evolutionary processes where, for example, a part of a protein is evolving at a different rate compared to other parts of the same protein. Identifying motifs are fundamental to understanding protein function and to discover putative binding interfaces.

Motifs can both be used to shed light on the evolutionary process underpinning the development of a protein family with respect to the protein’s function over time; it can also be used to produce a simplified view on the protein as a series of conserved motifs that together specify a proteins motif architecture. Although several approaches have been developed to address motif discovery on protein sequences, most are either focused on a given motif or finding motifs, such as signal peptides, that can be found in a general population of protein sequences.

Here, we developed a de novo motif discovery approach and applied to protein families that share a common ancestral protein; this resulted in a repository of motifs over an entire organism. This approach is focused on understanding the evolutionary processes that have acted on that protein family in a comparatively short evolutionary time. We developed and designed this approach as a software package which is written in Python and distributed via singularity containers [28] making it easy to install and use. We demonstrated the approach on GAS, an important human pathogen with a mortality rate of 25% at invasive infections. We also characterized the proteome-wide motif repository by comparing it to InterPro; furthermore, we analyzed the motif architecture for these proteins and discovered that the number of sequence per architecture is different for different cellular compartments.

Given the speed and flexibility of our approach, we believe it will be useful in breaking analyzing surface protein of pathogens as these proteins are under high selective pressure and therefore cannot be analyzed using more traditional approaches such as multiple-sequence alignments (MSAs). Our attempts to use various MSA algorithms failed due to high sequence variability in regions between motifs and the varying number of motifs. Also, motif searching approaches failed and only identified a small subset of the motifs that our approach discovered.

Conclusion

In this paper, we demonstrate a proof-of-principle approach to parsing large sequence families into motifs using a denovo-based greedy approach. This simple approach can easily handle situations where parts of proteins are repeated or re-arranged, and this can be time-consuming using other approaches. While this general approach can be applied to any bacteria, we used GAS as a model system to make a comprehensive motif repository of its proteins. We further analyzed M1 protein, one of the most important virulence factor of S. pyogenes to show the motif-based architectural analysis. We observe that we over-parse some domains, but also observe that many of these large domains are only partly conserved over the sequence collection. The result indicates that many of the newly discovered motifs are not always present together with adjacent motifs, indicating that they might have different and independent functions. Interestingly, many of our newly discovered motifs are not found in any of the emm1 strains, and some of these might be responsible for binding other ligands.

Abbreviations

CDC:

Centers for disease control and prevention

CDS:

Coding sequences

GAS:

Group A Streptococcus

MSA:

Multiple-sequence alignments

TM:

Transmembrane

TF:

Transcription factor

References

  1. 1

    O’Neill J. Antimicrobial resistance: tackling a crisis for the health and wealth of nations. Rev Antimicrob Resist. 2014.

  2. 2

    Forthal DN. Functions of antibodies. Microbiol Spectr. 2014; 2(4).

  3. 3

    Bork P, Koonin EV. Protein sequence motifs. Curr Opin Struct Biol. 1996; 6(3):366–76. https://doi.org/10.1016/S0959-440X(96)80057-1.

    CAS  Article  PubMed  Google Scholar 

  4. 4

    Ghosh P. Variation, indispensability, and masking in the m protein. Trends Microbiol. 2018; 26(2):132–44. https://doi.org/10.1016/j.tim.2017.08.002.

    CAS  Article  PubMed  Google Scholar 

  5. 5

    Charlotta S, Fredric C, Gunnar L. Binding of human plasma proteins to streptococcus pyogenes m protein determines the location of opsonic and non-opsonic epitopes. Mol Microbiol; 59(1):20–30. https://doi.org/10.1111/j.1365-2958.2005.04913.x. https://onlinelibrary.wiley.com/doi/pdf/10.1111/j.1365-2958.2005.04913.x.

    CAS  Article  PubMed  Google Scholar 

  6. 6

    Macheboeuf P, Buffalo C, Fu C-y, Zinkernagel AS, Cole JN, Johnson JE, Nizet V, Ghosh P. Streptococcal m1 protein constructs a pathological host fibrinogen network. Nature. 2011; 472(7341):64–68.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  7. 7

    Jones P, Binns D, Chang H-Y, Fraser M, Li W, McAnulla C, McWilliam H, Maslen J, Mitchell A, Nuka G, Pesseat S, Quinn AF, Sangrador-Vegas A, Scheremetjew M, Yong S-Y, Lopez R, Hunter S. Interproscan 5: genome-scale protein function classification. Bioinformatics. 2014; 30(9):1236–40. https://doi.org/10.1093/bioinformatics/btu031.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  8. 8

    Mitchell TJ. The pathogenesis of streptococcal infections: from tooth decay to meningitis. Nat Rev Microbiol. 2003; 1:219.

    CAS  Article  PubMed  Google Scholar 

  9. 9

    Karlsson C, Malmström L, Aebersold R, Malmström J. Proteome-wide selected reaction monitoring assays for the human pathogen streptococcus pyogenes. Nat Commun. 2012; 3:1301.

    Article  PubMed  PubMed Central  Google Scholar 

  10. 10

    Akesson P, Schmidt KH, Cooney J, Björck Larsrck L. M1 protein and protein h: Iggfc- and albumin-binding streptococcal surface proteins encoded by adjacent genes. Biochem J. 1994; 300(3):877–86.

    Article  PubMed  PubMed Central  Google Scholar 

  11. 11

    Hauri S, Khakzad H, Happonen L, Teleman J, Malmström J, Malmström L. Rapid determination of quaternary protein structures in complex biological samples. Nat Commun. 2019; 10(1):192. https://doi.org/10.1038/s41467-018-07986-1.

  12. 12

    McMillan DJ, Drèze P-A, Vu T, Bessen DE, Guglielmini J, Steer AC, Carapetis JR, Van Melderen L, Sriprakash KS, Smeesters PR, the M Protein Study Group. Updated model of group a streptococcus m proteins based on a comprehensive worldwide study. Clin Microbiol Infect Off Publ Eur Soc Clin Microbiol Infect Dis. 2013; 19(5):222–9. https://doi.org/10.1111/1469-0691.12134.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  13. 13

    Buffalo CZ, Bahn-Suh AJ, Hirakis SP, Biswas T, Amaro RE, Nizet V, Ghosh P. Conserved patterns hidden within group a streptococcus m protein hypervariability are responsible for recognition of human c4b-binding protein. Nat Microbiol. 2016; 1:16155. https://doi.org/10.1038/nmicrobiol.2016.155.

  14. 14

    Dale JB, Smeesters PR, Courtney HS, Penfound TA, Hohn CM, Smith JC, Baudry JY. Structure-based design of broadly protective group a streptococcal m protein-based vaccines. Vaccine. 2017; 35(1):19–26. https://doi.org/10.1016/j.vaccine.2016.11.065.

    CAS  Article  PubMed  Google Scholar 

  15. 15

    Sanderson-Smith M, De Oliveira DMP, Guglielmini J, McMillan DJ, Vu T, Holien JK, Henningham A, Steer AC, Bessen DE, Dale JB, Curtis N, Beall BW, Walker MJ, Parker MW, Carapetis JR, Van Melderen L, Sriprakash KS, Smeesters PR. A systematic and functional classification of streptococcus pyogenes that serves as a new tool for molecular typing and vaccine development. J Infect Dis. 2014; 210(8):1325–38. https://doi.org/10.1093/infdis/jiu260.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  16. 16

    Sjöholm K, Kilsgård O, Teleman J, Happonen L, Malmström L, Malmström J. Targeted proteomics and absolute protein quantification for the construction of a stoichiometric host-pathogen surface density model. Mol Cell Proteome. 2017; 16(4 suppl 1):29–41.

    Article  Google Scholar 

  17. 17

    Wattam AR, Davis JJ, Assaf R, Boisvert S, Brettin T, Bun C, Conrad N, Dietrich EM, Disz T, Gabbard JL, Gerdes S, Henry CS, Kenyon RW, Machi D, Mao C, Nordberg EK, Olsen GJ, Murphy-Olson DE, Olson R, Overbeek R, Parrello B, Pusch GD, Shukla M, Vonstein V, Warren A, Xia F, Yoo H, Stevens RL. Improvements to patric, the all-bacterial bioinformatics database and analysis resource center. Nucleic Acids Res. 2017; 45(D1):535–42.

    Article  Google Scholar 

  18. 18

    Meyer F, Overbeek R, Rodriguez A. Figfams: yet another set of protein families. Nucleic Acids Res. 2009; 37(20):6643–54. https://doi.org/10.1093/nar/gkp698.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  19. 19

    Bailey TL, Boden M, Buske FA, Frith M, Grant CE, Clementi L, Ren J, Li WW, Noble WS. Meme suite: tools for motif discovery and searching. Nucleic Acids Res. 2009; 37(suppl 2):202–8.

    Article  Google Scholar 

  20. 20

    Bailey TL, Elkan C. Fitting a Mixture Model by Expectation Maximization to Discover Motifs in Biopolymers, vol. 2. Stanford, California: AAAI Press; 1994, pp. 28–36. Proceedings of the Second International Conference on Intelligent Systems for Molecular Biology.

    Google Scholar 

  21. 21

    Grant CE, Bailey TL, Noble WS. Fimo: scanning for occurrences of a given motif. Bioinformatics. 2011; 27(7):1017–8.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  22. 22

    Frith MC, Saunders NFW, Kobe B, Bailey TL. Discovering sequence motifs with arbitrary insertions and deletions. PLOS Comput Biol. 2008; 4(5):1–12. https://doi.org/10.1371/journal.pcbi.1000071.

    Article  PubMed  PubMed Central  Google Scholar 

  23. 23

    Staden R. Staden: Searching for Motifs in Nucleic Acid Sequences. Totowa, NJ: Springer; 1994, pp. 93–102. https://doi.org/10.1385/0-89603-276-0:93.

    Chapter  Google Scholar 

  24. 24

    Storey JD. A direct approach to false discovery rates. J R Stat Soc Ser B Stat Methodol. 2002; 64(3):479–98. https://doi.org/10.1111/1467-9868.00346.

    Article  Google Scholar 

  25. 25

    Krogh A, Larsson B, von Heijne G, Sonnhammer ELL. Predicting transmembrane protein topology with a hidden markov model: application to complete genomes11edited by f. cohen. J Mol Biol. 2001; 305(3):567–80. https://doi.org/10.1006/jmbi.2000.4315.

    CAS  Article  PubMed  Google Scholar 

  26. 26

    Novichkov PS, Laikova ON, Novichkova ES, Gelfand MS, Arkin AP, Dubchak I, Rodionov DA. Regprecise: a database of curated genomic inferences of transcriptional regulatory interactions in prokaryotes. Nucleic Acids Res. 2010; 38(suppl_1):111–8.

    Article  Google Scholar 

  27. 27

    Maffioletti S, Murri R. Gc3pie: A python framework for high-throughput computing. Proc EGI Community Forum 2012/EMI Second Tech Conf (EGICF12-EMITC2). 2012:143.

  28. 28

    Kurtzer GM, Sochat V, Bauer MW. Singularity: Scientific containers for mobility of compute. PLoS ONE. 2017; 12(5):1–20. https://doi.org/10.1371/journal.pone.0177459.

    Article  PubMed  PubMed Central  Google Scholar 

  29. 29

    Mukhyala K, Masselot A. Visualization of protein sequence features using javascript and svg with pviz.js. Bioinformatics. 2014; 30(23):3408–9.

    CAS  Article  PubMed  Google Scholar 

  30. 30

    Malmström L, Bakochi A, Svensson G, Kilsgård O, Lantz H, Petersson AC, Hauri S, Karlsson C, Malmström J. Quantitative proteogenomics of human pathogens using dia-ms. J Proteomics. 2015; 129:98–107. https://doi.org/10.1016/j.jprot.2015.09.012. Special Issue : Computational Proteomics.

    Article  PubMed  Google Scholar 

  31. 31

    Fu L, Niu B, Zhu Z, Wu S, Li W. Cd-hit: accelerated for clustering the next-generation sequencing data. Bioinformatics. 2012; 28(23):3150–2.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  32. 32

    Ferretti JJ, McShan WM, Ajdic D, Savic DJ, Savic G, Lyon K, Primeaux C, Sezate S, Suvorov AN, Kenton S, Lai HS, Lin SP, Qian Y, Jia HG, Najar FZ, Ren Q, Zhu H, Song L, White J, Yuan X, Clifton SW, Roe BA, McLaughlin R. Complete genome sequence of an m1 strain of streptococcus pyogenes. Proc Natl Acad Sci. 2001; 98(8):4658–63. https://doi.org/10.1073/pnas.071559398. http://www.pnas.org/content/98/8/4658.full.pdf.

    CAS  Article  Google Scholar 

Download references

Acknowledgements

We thank S3IT, University of Zurich for computational infrastructure and technical support.

Funding

This work was supported by the Swiss National Science Foundation (grant no. SNF 200021 160188) and the Knut And Alice Wallenberg Foundation (grant no. KAW 2016.0023). Publication of this article was sponsored by the Knut And Alice Wallenberg Foundation (grant no. KAW 2016.0023).

Availability of data and materials

Not applicable.

About this supplement

This article has been published as part of BMC Bioinformatics Volume 20 Supplement 4, 2019: Methods, tools and platforms for Personalized Medicine in the Big Data Era (NETTAB 2017). The full contents of the supplement are available online at https://bmcbioinformatics.biomedcentral.com/articles/supplements/volume-20-supplement-4.

Author information

Affiliations

Authors

Contributions

HK, JM and LM designed the study, analyzed the data and wrote the manuscript. HK and LM implemented the software. All of the authors have read and approved the final manuscript.

Corresponding author

Correspondence to Lars Malmström.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional file

Additional file 1

Supplementary Table 1. The logo of all 20 motifs obtained from our approach for M protein family is listed in a multi-page table in the supplementary material. (PDF 1075 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver(http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Khakzad, H., Malmström, J. & Malmström, L. Greedy de novo motif discovery to construct motif repositories for bacterial proteomes. BMC Bioinformatics 20, 141 (2019). https://doi.org/10.1186/s12859-019-2686-8

Download citation

Keywords

  • De novo motif discovery
  • Infectious diseases
  • Group A streptococcus