A HMM-type anvi’o artifact. This artifact is typically provided by the user for anvi’o to import into its databases, process, and/or use.
Back to the main page of anvi’o programs and artifacts.
An HMM source is a collection of one or more hidden Markov models. HMM sources can be used to identify and recover genes in a contigs-db that match to those described in the model.
Models in a given HMM source can be searched in a given contigs-db via the program anvi-run-hmms which would yield an hmm-hits artifact. An anvi’o installation will include multiple HMM sources by default. But HMMs for any set of genes can also be put together by the end user and run on any anvi’o contigs-db.
HMM hits in a contigs-db for a given hmm-source source will be accessible to anvi’o programs globally. Sequences that match to HMM hits can be recovered in an aligned or non-aligned fashion as fasta files for downstream analyses including phylogenomics, they can be displayed in anvi’o interfaces, reported in summary outputs, and so on.
The default HMM sources in anvi’o include:
Peptidase_A8occurred in Lee collection (as they were exceptionally redundant or rare among MAGs from various habitats), and includes
Ribosomal_L2to make it more compatible with Hug et al’s set of ribosomal proteins.
Similar to Ribosomal RNAs, anvi’o can also identify Transfer RNAs. Even though Transfer RNAs will also appear as an HMM source for all downstream analyses, their initial identification will require running anvi-scan-trnas program on a contigs-db.
The user can employ additional HMM sources to identify matching genes in a given contigs-db.
Any directory with expected files in it will serve as an HMM source:
anvi-run-hmms -c contigs-db \ --hmm-source /PATH/TO/USER-HMM-DIRECTORY/
Anvi’o will expect the HMM source directory to contain six files (see this for an example directory). These files are explicitly defined as follows:
genes.hmm, and finally (3) compress this file using
accession(gene accession number (can be anything unique)), and
hmmsource(source of HMM profiles listed in genes.hmm.gz). The list of gene names in this file must perfectly match to the list of gene names in genes.hmm.gz.
RNA. Context can be
CONTIG. The content of this file should be any combination of one alphabet and one context term. For instance, if the content of this file is
AA:GENE, anvi’o will search genes amino acid sequences, and so on. An exception is
AA:CONTIG, which is an improper target since anvi’o can’t translate contigs to amino acid sequences. See this for more details. Please note that HMMs that target
DNA:CONTIGwill result in new gene calls in the contigs database to describe their hits.
It is also possible to generate an anvi’o compatible HMMs directory for a given set of PFAM accession ids. For instance, the following command will result in a new directory that can be used immediately with the program anvi-run-hmms:
anvi-script-pfam-accessions-to-hmms-directory --pfam-accessions-list PF00705 PF00706 \ -O AD_HOC_HMMs
These IDs can be given through the command line as a list, or through an input file where every line is a unique accession id.
An example. Let’s assume we have a genome or a metagenome that looks like this:
And we wish to identify locations of genes that match to this model: http://pfam.xfam.org/family/PF06603
One can run this command:
anvi-script-pfam-accessions-to-hmms-directory --pfam-accessions-list PF06603 \ -O UpxZ
which would createa a directory called
UpxZ. Then, one would run this command to find matches to this model in a given contigs database:
anvi-run-hmms -c CONTIGS.db \ -H UpxZ/ \ --num-threads 4
Now it is possible to get the sequences matching to this model:
anvi-get-sequences-for-hmm-hits -c CONTIGS.db \ --hmm-source UpxZ \ -o UpxZ.fa
Contigs DB ...................................: Initialized: CONTIGS.db (v. 19) Hits .........................................: 8 hits for 1 source(s) Mode .........................................: DNA sequences Genes are concatenated .......................: False Output .......................................: UpxZ.fa
Or see where they are by visualizing the project using again:
anvi-interactive -p PROFILE.db \ -c CONTIGS.db
Edit this file to update this information.