Rational Design of Profile HMMs for Sensitive and Specific Sequence Detection with Case Studies Applied to Viruses, Bacteriophages, and Casposons

Profile hidden Markov models (HMMs) are a powerful way of modeling biological sequence diversity and constitute a very sensitive approach to detecting divergent sequences. Here, we report the development of protocols for the rational design of profile HMMs. These methods were implemented on TABAJARA...

Full description

Bibliographic Details
Main Authors: Liliane S. Oliveira, Alejandro Reyes, Bas E. Dutilh, Arthur Gruber
Format: Article
Language:English
Published: MDPI AG 2023-02-01
Series:Viruses
Subjects:
Online Access:https://www.mdpi.com/1999-4915/15/2/519
_version_ 1797617788198584320
author Liliane S. Oliveira
Alejandro Reyes
Bas E. Dutilh
Arthur Gruber
author_facet Liliane S. Oliveira
Alejandro Reyes
Bas E. Dutilh
Arthur Gruber
author_sort Liliane S. Oliveira
collection DOAJ
description Profile hidden Markov models (HMMs) are a powerful way of modeling biological sequence diversity and constitute a very sensitive approach to detecting divergent sequences. Here, we report the development of protocols for the rational design of profile HMMs. These methods were implemented on TABAJARA, a program that can be used to either detect all biological sequences of a group or discriminate specific groups of sequences. By calculating position-specific information scores along a multiple sequence alignment, TABAJARA automatically identifies the most informative sequence motifs and uses them to construct profile HMMs. As a proof-of-principle, we applied TABAJARA to generate profile HMMs for the detection and classification of two viral groups presenting different evolutionary rates: bacteriophages of the <i>Microviridae</i> family and viruses of the <i>Flavivirus</i> genus. We obtained conserved models for the generic detection of any <i>Microviridae</i> or <i>Flavivirus</i> sequence, and profile HMMs that can specifically discriminate <i>Microviridae</i> subfamilies or <i>Flavivirus</i> species. In another application, we constructed Cas1 endonuclease-derived profile HMMs that can discriminate CRISPRs and casposons, two evolutionarily related transposable elements. We believe that the protocols described here, and implemented on TABAJARA, constitute a generic toolbox for generating profile HMMs for the highly sensitive and specific detection of sequence classes.
first_indexed 2024-03-11T08:00:44Z
format Article
id doaj.art-6ea8278ddac74c389a5349c06d103afc
institution Directory Open Access Journal
issn 1999-4915
language English
last_indexed 2024-03-11T08:00:44Z
publishDate 2023-02-01
publisher MDPI AG
record_format Article
series Viruses
spelling doaj.art-6ea8278ddac74c389a5349c06d103afc2023-11-16T23:50:07ZengMDPI AGViruses1999-49152023-02-0115251910.3390/v15020519Rational Design of Profile HMMs for Sensitive and Specific Sequence Detection with Case Studies Applied to Viruses, Bacteriophages, and CasposonsLiliane S. Oliveira0Alejandro Reyes1Bas E. Dutilh2Arthur Gruber3Department of Parasitology, Instituto de Ciências Biomédicas, Universidade de São Paulo, São Paulo 05508-000, SP, BrazilMax Planck Tandem Group in Computational Biology, Department of Biological Sciences, Universidad de los Andes, Bogotá 111711, ColombiaInstitute of Biodiversity, Faculty of Biological Sciences, Cluster of Excellence Balance of the Microverse, Friedrich-Schiller-University Jena, 07743 Jena, GermanyDepartment of Parasitology, Instituto de Ciências Biomédicas, Universidade de São Paulo, São Paulo 05508-000, SP, BrazilProfile hidden Markov models (HMMs) are a powerful way of modeling biological sequence diversity and constitute a very sensitive approach to detecting divergent sequences. Here, we report the development of protocols for the rational design of profile HMMs. These methods were implemented on TABAJARA, a program that can be used to either detect all biological sequences of a group or discriminate specific groups of sequences. By calculating position-specific information scores along a multiple sequence alignment, TABAJARA automatically identifies the most informative sequence motifs and uses them to construct profile HMMs. As a proof-of-principle, we applied TABAJARA to generate profile HMMs for the detection and classification of two viral groups presenting different evolutionary rates: bacteriophages of the <i>Microviridae</i> family and viruses of the <i>Flavivirus</i> genus. We obtained conserved models for the generic detection of any <i>Microviridae</i> or <i>Flavivirus</i> sequence, and profile HMMs that can specifically discriminate <i>Microviridae</i> subfamilies or <i>Flavivirus</i> species. In another application, we constructed Cas1 endonuclease-derived profile HMMs that can discriminate CRISPRs and casposons, two evolutionarily related transposable elements. We believe that the protocols described here, and implemented on TABAJARA, constitute a generic toolbox for generating profile HMMs for the highly sensitive and specific detection of sequence classes.https://www.mdpi.com/1999-4915/15/2/519viral classificationviral detectionprofile HMMssequence classesmutual information theorysequence entropy
spellingShingle Liliane S. Oliveira
Alejandro Reyes
Bas E. Dutilh
Arthur Gruber
Rational Design of Profile HMMs for Sensitive and Specific Sequence Detection with Case Studies Applied to Viruses, Bacteriophages, and Casposons
Viruses
viral classification
viral detection
profile HMMs
sequence classes
mutual information theory
sequence entropy
title Rational Design of Profile HMMs for Sensitive and Specific Sequence Detection with Case Studies Applied to Viruses, Bacteriophages, and Casposons
title_full Rational Design of Profile HMMs for Sensitive and Specific Sequence Detection with Case Studies Applied to Viruses, Bacteriophages, and Casposons
title_fullStr Rational Design of Profile HMMs for Sensitive and Specific Sequence Detection with Case Studies Applied to Viruses, Bacteriophages, and Casposons
title_full_unstemmed Rational Design of Profile HMMs for Sensitive and Specific Sequence Detection with Case Studies Applied to Viruses, Bacteriophages, and Casposons
title_short Rational Design of Profile HMMs for Sensitive and Specific Sequence Detection with Case Studies Applied to Viruses, Bacteriophages, and Casposons
title_sort rational design of profile hmms for sensitive and specific sequence detection with case studies applied to viruses bacteriophages and casposons
topic viral classification
viral detection
profile HMMs
sequence classes
mutual information theory
sequence entropy
url https://www.mdpi.com/1999-4915/15/2/519
work_keys_str_mv AT lilianesoliveira rationaldesignofprofilehmmsforsensitiveandspecificsequencedetectionwithcasestudiesappliedtovirusesbacteriophagesandcasposons
AT alejandroreyes rationaldesignofprofilehmmsforsensitiveandspecificsequencedetectionwithcasestudiesappliedtovirusesbacteriophagesandcasposons
AT basedutilh rationaldesignofprofilehmmsforsensitiveandspecificsequencedetectionwithcasestudiesappliedtovirusesbacteriophagesandcasposons
AT arthurgruber rationaldesignofprofilehmmsforsensitiveandspecificsequencedetectionwithcasestudiesappliedtovirusesbacteriophagesandcasposons