AnnoView enables large-scale analysis, comparison, and visualization of microbial gene neighborhoods.

The analysis and comparison of gene neighborhoods is a powerful approach for exploring microbial genome structure, function, and evolution. Although numerous tools exist for genome visualization and comparison, genome exploration across large genomic databases or user-generated datasets remains a challenge. Here, we introduce AnnoView, a web server designed for interactive exploration of gene neighborhoods across the bacterial and archaeal tree of life. Our server offers users the ability to identify, compare, and visualize gene neighborhoods of interest from 30 238 bacterial genomes and 1672 archaeal genomes, through integration with the comprehensive Genome Taxonomy Database and AnnoTree databases. Identified gene neighborhoods can be visualized using pre-computed functional annotations from different sources such as KEGG, Pfam and TIGRFAM, or clustered based on similarity. Alternatively, users can upload and explore their own custom genomic datasets in GBK, GFF or CSV format, or use AnnoView as a genome browser for relatively small genomes (e.g. viruses and plasmids). Ultimately, we anticipate that AnnoView will catalyze biological discovery by enabling user-friendly search, comparison, and visualization of genomic data. AnnoView is available at http://annoview.uwaterloo.ca.

Nasopharyngeal metatranscriptomics reveals host-pathogen signatures of pediatric sinusitis.

Acute sinusitis (AS) is the fifth leading cause of antibiotic prescriptions in children. Distinguishing bacterial AS from common viral upper respiratory infections in children is crucial to prevent unnecessary antibiotic use but is challenging with current diagnostic methods. Despite its speed and cost, untargeted RNA sequencing of clinical samples from children with suspected AS has the potential to overcome several limitations of other methods. However, the utility of sequencing-based approaches in analysis of AS has not been fully explored. Here, we performed RNA-seq of nasopharyngeal samples from 221 children with clinically diagnosed AS to characterize their pathogen and host-response profiles. Results from RNA-seq were compared with those obtained using culture for three common bacterial pathogens and qRT-PCR for 12 respiratory viruses. Metatranscriptomic pathogen detection showed high concordance with culture or qRT-PCR, showing 87%/81% sensitivity (sens) / specificity (spec) for detecting bacteria, and 86%/92% (sens/spec) for viruses, respectively. We also detected an additional 22 pathogens not tested for in the clinical panel, and identified plausible pathogens in 11/19 (58%) of cases where no organism was detected by culture or qRT-PCR. We assembled genomes of 205 viruses across the samples including novel strains of coronaviruses, respiratory syncytial virus (RSV), and enterovirus D68. By analyzing host gene expression, we identified host-response signatures that distinguished bacterial and viral infections and correlated with pathogen abundance. Ultimately, our study demonstrates the potential of untargeted metatranscriptomics for in depth analysis of the etiology of AS, comprehensive host-response profiling, and using these together to work towards optimized patient care.

Genomic and transcriptomic characterization of delta SARS-CoV-2 infection in free-ranging white-tailed deer (Odocoileus virginianus).

White-tailed deer (WTD) are susceptible to SARS-CoV-2 and represent an important species for surveillance. Samples from WTD (n = 258) collected in November 2021 from Québec, Canada were analyzed for SARS-CoV-2 RNA. We employed viral genomics and host transcriptomics to further characterize infection and investigate host response. We detected Delta SARS-CoV-2 (B.1.617.2) in WTD from the Estrie region; sequences clustered with human sequences from October 2021 from Vermont, USA, which borders this region. Mutations in the S-gene and a deletion in ORF8 were detected. Host expression patterns in SARS-CoV-2 infected WTD were associated with the innate immune response, including signaling pathways related to anti-viral, pro- and anti-inflammatory signaling, and host damage. We found limited correlation between genes associated with innate immune response from human and WTD nasal samples, suggesting differences in responses to SARS-CoV-2 infection. Our findings provide preliminary insights into host response to SARS-CoV-2 infection in naturally infected WTD.

Genomic classification and antimicrobial resistance profiling of Streptococcus pneumoniae and Haemophilus influenza isolates associated with paediatric otitis media and upper respiratory infection.

Acute otitis media (AOM) is the most common childhood bacterial infectious disease requiring antimicrobial therapy. Most cases of AOM are caused by translocation of Streptococcus pneumoniae or Haemophilus influenzae from the nasopharynx to the middle ear during an upper respiratory tract infection (URI). Ongoing genomic surveillance of these pathogens is important for vaccine design and tracking of emerging variants, as well as for monitoring patterns of antibiotic resistance to inform treatment strategies and stewardship.In this work, we examined the ability of a genomics-based workflow to determine microbiological and clinically relevant information from cultured bacterial isolates obtained from patients with AOM or an URI. We performed whole genome sequencing (WGS) and analysis of 148 bacterial isolates cultured from the nasopharynx (N = 124, 94 AOM and 30 URI) and ear (N = 24, all AOM) of 101 children aged 6-35 months presenting with AOM or an URI. We then performed WGS-based sequence typing and antimicrobial resistance profiling of each strain and compared results to those obtained from traditional microbiological phenotyping.WGS of clinical isolates resulted in 71 S. pneumoniae genomes and 76 H. influenzae genomes. Multilocus sequencing typing (MSLT) identified 33 sequence types for S. pneumoniae and 19 predicted serotypes including the most frequent serotypes 35B and 3. Genome analysis predicted 30% of S. pneumoniae isolates to have complete or intermediate penicillin resistance. AMR predictions for S. pneumoniae isolates had strong agreement with clinical susceptibility testing results for beta-lactam and non beta-lactam antibiotics, with a mean sensitivity of 93% (86-100%) and a mean specificity of 98% (94-100%). MLST identified 29 H. influenzae sequence types. Genome analysis identified beta-lactamase genes in 30% of H. influenzae strains, which was 100% in agreement with clinical beta-lactamase testing. We also identified a divergent highly antibiotic-resistant strain of S. pneumoniae, and found its closest sequenced strains, also isolated from nasopharyngeal samples from over 15 years ago.Ultimately, our work provides the groundwork for clinical WGS-based workflows to aid in detection and analysis of H. influenzae and S. pneumoniae isolates.

Ancient Clostridium DNA and variants of tetanus neurotoxins associated with human archaeological remains.

The analysis of microbial genomes from human archaeological samples offers a historic snapshot of ancient pathogens and provides insights into the origins of modern infectious diseases. Here, we analyze metagenomic datasets from 38 human archaeological samples and identify bacterial genomic sequences related to modern-day Clostridium tetani, which produces the tetanus neurotoxin (TeNT) and causes the disease tetanus. These genomic assemblies had varying levels of completeness, and a subset of them displayed hallmarks of ancient DNA damage. Phylogenetic analyses revealed known C. tetani clades as well as potentially new Clostridium lineages closely related to C. tetani. The genomic assemblies encode 13 TeNT variants with unique substitution profiles, including a subgroup of TeNT variants found exclusively in ancient samples from South America. We experimentally tested a TeNT variant selected from an ancient Chilean mummy sample and found that it induced tetanus muscle paralysis in mice, with potency comparable to modern TeNT. Thus, our ancient DNA analysis identifies DNA from neurotoxigenic C. tetani in archaeological human samples, and a novel variant of TeNT that can cause disease in mammals.

Phylogenomic Analysis of 155 Helminth Species Reveals Widespread Absence of Oxygen Metabolic Capacity.

The terminal electron acceptor of most aerobic respiratory chains, cytochrome c oxidase (COX), has been highly conserved throughout evolution, from aerobic prokaryotes to complex eukaryotes. Oxygen metabolism in parasitic helminths differs significantly from that of most aerobic eukaryotes, as these organisms can switch between aerobic and anaerobic metabolisms throughout their life cycles. Early studies suggested a lack of COX activity in certain parasitic helminths, and the role of COX in helminth mitochondria remains unclear. To determine whether a functional COX is widely present in helminths, we analyzed the phylogenetic distribution of oxygen metabolism systems across 155 helminth genomes, investigating three distinct sets of protein-coding genes involved in different aspects of oxygen metabolism: COX and its assembly factors, peroxisomes, and the most abundant reactive oxygen species (ROS)-metabolizing proteins. While glycolytic and citric acid cycle enzymes are highly conserved in helminthic species, we observed an apparent widespread absence of essential COX genes across 52% of helminth species investigated. While the most common proteins involved in the defense against ROS are highly maintained across virtually all lineages, we also observed an apparent absence of essential peroxisomal protein-coding genes in 42% of species investigated. Our results suggest that a subset of parasitic helminths utilize oxygen differently from related, nonparasitic species such as Caenorhabditis elegans, with significant differences in their mitochondrial electron transport chains and peroxisomes. The identification of substantive differences between parasite and host metabolism offers a new avenue for the development of anthelmintic agents that could target these divergent pathways.

Functional Profiling and Evolutionary Analysis of a Marine Microalgal Virus Pangenome.

Phycodnaviridae are large double-stranded DNA viruses, which facilitate studies of host-virus interactions and co-evolution due to their prominence in algal infection and their role in the life cycle of algal blooms. However, the genomic interpretation of these viruses is hampered by a lack of functional information, stemming from the surprising number of hypothetical genes of unknown function. It is also unclear how many of these genes are widely shared within the clade. Using one of the most extensively characterized genera, Coccolithovirus, as a case study, we combined pangenome analysis, multiple functional annotation tools, AlphaFold structural modeling, and literature analysis to compare the core and accessory pangenome and assess support for novel functional predictions. We determined that the Coccolithovirus pangenome shares 30% of its genes with all 14 strains, making up the core. Notably, 34% of its genes were found in at most three strains. Core genes were enriched in early expression based on a transcriptomic dataset of Coccolithovirus EhV-201 algal infection, were more likely to be similar to host proteins than the non-core set, and were more likely to be involved in vital functions such as replication, recombination, and repair. In addition, we generated and collated annotations for the EhV representative EhV-86 from 12 different annotation sources, building up information for 142 previously hypothetical and putative membrane proteins. AlphaFold was further able to predict structures for 204 EhV-86 proteins with a modelling accuracy of good-high. These functional clues, combined with generated AlphaFold structures, provide a foundational framework for the future characterization of this model genus (and other giant viruses) and a further look into the evolution of the Coccolithovirus proteome.

PathFams: statistical detection of pathogen-associated protein domains.

BACKGROUND: A substantial fraction of genes identified within bacterial genomes encode proteins of unknown function. Identifying which of these proteins represent potential virulence factors, and mapping their key virulence determinants, is a challenging but important goal. RESULTS: To facilitate virulence factor discovery, we performed a comprehensive analysis of 17,929 protein domain families within the Pfam database, and scored them based on their overrepresentation in pathogenic versus non-pathogenic species, taxonomic distribution, relative abundance in metagenomic datasets, and other factors. CONCLUSIONS: We identify pathogen-associated domain families, candidate virulence factors in the human gut, and eukaryotic-like mimicry domains with likely roles in virulence. Furthermore, we provide an interactive database called PathFams to allow users to explore pathogen-associated domains as well as identify pathogen-associated domains and domain architectures in user-uploaded sequences of interest. PathFams is freely available at https://pathfams.uwaterloo.ca .

PhyloCorrelate: inferring bacterial gene-gene functional associations through large-scale phylogenetic profiling.

MOTIVATION: Statistical detection of co-occurring genes across genomes, known as 'phylogenetic profiling', is a powerful bioinformatic technique for inferring gene-gene functional associations. However, this can be a challenging task given the size and complexity of phylogenomic databases, difficulty in accounting for phylogenetic structure, inconsistencies in genome annotation and substantial computational requirements. RESULTS: We introduce PhyloCorrelate-a computational framework for gene co-occurrence analysis across large phylogenomic datasets. PhyloCorrelate implements a variety of co-occurrence metrics including standard correlation metrics and model-based metrics that account for phylogenetic history. By combining multiple metrics, we developed an optimized score that exhibits a superior ability to link genes with overlapping GO terms and KEGG pathways, enabling gene function prediction. Using genomic and functional annotation data from the Genome Taxonomy Database and AnnoTree, we performed all-by-all comparisons of gene occurrence profiles across the bacterial tree of life, totaling 154 217 052 comparisons for 28 315 genes across 27 372 bacterial genomes. All predictions are available in an online database, which instantaneously returns the top correlated genes for any PFAM, TIGRFAM or KEGG query. In total, PhyloCorrelate detected 29 762 high confidence associations between bacterial gene/protein pairs, and generated functional predictions for 834 DUFs and proteins of unknown function. AVAILABILITYAND IMPLEMENTATION: PhyloCorrelate is available as a web-server at phylocorrelate.uwaterloo.ca as well as an R package for analysis of custom datasets. We anticipate that PhyloCorrelate will be broadly useful as a tool for predicting function and interactions for gene families. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Slr4, a newly identified S-layer protein from marine Gammaproteobacteria, is a major biofilm matrix component.

S-layers are paracrystalline proteinaceous lattices that surround prokaryotic cells, forming a critical interface between the cells and their extracellular environment. Here, we report the discovery of a novel S-layer protein present in the Gram-negative marine organism, Pseudoalteromonas tunicata D2. An uncharacterized protein (EAR28894) was identified as the most abundant protein in planktonic cultures and biofilms. Bioinformatic methods predicted a beta-helical structure for EAR28894 similar to the Caulobacter S-layer protein, RsaA, despite sharing less than 20% sequence identity. Transmission electron microscopy revealed that purified EAR28894 protein assembled into paracrystalline sheets with a unique square lattice symmetry and a unit cell spacing of ~9.1 nm. An S-layer was found surrounding the outer membrane in wild-type cells and completely removed from cells in an EAR28894 deletion mutant. S-layer material also appeared to be "shed" from wild-type cells and was highly abundant in the extracellular matrix where it is associated with outer membrane vesicles and other matrix components. EAR28894 and its homologs form a new family of S-layer proteins that are widely distributed in Gammaproteobacteria including species of Pseudoalteromonas and Vibrio, and found exclusively in marine metagenomes. We propose the name Slr4 for this novel protein family.

Gene expression and in situ protein profiling of candidate SARS-CoV-2 receptors in human airway epithelial cells and lung tissue.

In December 2019, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) emerged, causing the coronavirus disease 2019 (COVID-19) pandemic. SARS-CoV, the agent responsible for the 2003 SARS outbreak, utilises angiotensin-converting enzyme 2 (ACE2) and transmembrane serine protease 2 (TMPRSS2) host molecules for viral entry. ACE2 and TMPRSS2 have recently been implicated in SARS-CoV-2 viral infection. Additional host molecules including ADAM17, cathepsin L, CD147 and GRP78 may also function as receptors for SARS-CoV-2.To determine the expression and in situ localisation of candidate SARS-CoV-2 receptors in the respiratory mucosa, we analysed gene expression datasets from airway epithelial cells of 515 healthy subjects, gene promoter activity analysis using the FANTOM5 dataset containing 120 distinct sample types, single cell RNA sequencing (scRNAseq) of 10 healthy subjects, proteomic datasets, immunoblots on multiple airway epithelial cell types, and immunohistochemistry on 98 human lung samples.We demonstrate absent to low ACE2 promoter activity in a variety of lung epithelial cell samples and low ACE2 gene expression in both microarray and scRNAseq datasets of epithelial cell populations. Consistent with gene expression, rare ACE2 protein expression was observed in the airway epithelium and alveoli of human lung, confirmed with proteomics. We present confirmatory evidence for the presence of TMPRSS2, CD147 and GRP78 protein in vitro in airway epithelial cells and confirm broad in situ protein expression of CD147 and GRP78 in the respiratory mucosa.Collectively, our data suggest the presence of a mechanism dynamically regulating ACE2 expression in human lung, perhaps in periods of SARS-CoV-2 infection, and also suggest that alternative receptors for SARS-CoV-2 exist to facilitate initial host cell infection.

Time Series Resolution of the Fish Necrobiome Reveals a Decomposer Succession Involving Toxigenic Bacterial Pathogens.

Despite progress understanding microbial communities involved in terrestrial vertebrate decomposition, little is known about the microbial decomposition of aquatic vertebrates from a functional and environmental context. Here, we analyzed temporal changes in the "necrobiome" of rainbow darters, which are common North American fish that are sensitive indicators of water quality. By combining 16S rRNA gene and shotgun metagenomic sequence data from four time points, we studied the progression of decomposers from both taxonomic and functional perspectives. The 16S rRNA gene profiles revealed strong community succession, with early decomposition stages associated with Aeromonas and Clostridium taxa and later stages dominated by members of the Rikenellaceae (i.e., Alistipes/Acetobacteroides genera). These results were reproducible and independent of environmental perturbation, given that exposure to wastewater treatment plant effluent did not substantially influence the necrobiome composition of fish or the associated water sample microbiota. Metagenomic analysis revealed significant changes throughout decomposition in degradation pathways for amino acids, carbohydrates/glycans, and other compounds, in addition to putrefaction pathways for production of putrescine, cadaverine, and indole. Binning of contigs confirmed a predominance of Aeromonas genome assemblies, including those from novel strains related to the pathogen Aeromonas veronii These bins of Aeromonas genes also encoded known hemolysin toxins (e.g., aerolysin) that were particularly abundant early in the process, potentially contributing to host cell lysis during decomposition. Overall, our results demonstrate that wild-caught fish have a reproducible decomposer succession and that the fish necrobiome serves as a potential source of putative pathogens and toxigenic bacteria.IMPORTANCE The microbial decomposition of animal tissues is an important ecological process that impacts nutrient cycling in natural environments. We studied the microbial decomposition of a common North American fish (rainbow darters) over four time points, combining 16S rRNA gene and shotgun metagenomic sequence data to obtain both taxonomic and functional perspectives. Our data revealed a strong community succession that was reproduced across different fish and environments. Decomposition time point was the main driver of community composition and functional potential; fish environmental origin (upstream or downstream of a wastewater treatment plant) had a secondary effect. We also identified strains related to the putative pathogen Aeromonas veronii as dominant members of the decomposition community. These bacteria peaked early in decomposition and coincided with the metagenomic abundance of hemolytic toxin genes. Our work reveals a strong decomposer succession in wild-caught fish, providing functional and taxonomic insights into the vertebrate necrobiome.

An assessment of genome annotation coverage across the bacterial tree of life.

Although gene-finding in bacterial genomes is relatively straightforward, the automated assignment of gene function is still challenging, resulting in a vast quantity of hypothetical sequences of unknown function. But how prevalent are hypothetical sequences across bacteria, what proportion of genes in different bacterial genomes remain unannotated, and what factors affect annotation completeness? To address these questions, we surveyed over 27 000 bacterial genomes from the Genome Taxonomy Database, and measured genome annotation completeness as a function of annotation method, taxonomy, genome size, 'research bias' and publication date. Our analysis revealed that 52 and 79 % of the average bacterial proteome could be functionally annotated based on protein and domain-based homology searches, respectively. Annotation coverage using protein homology search varied significantly from as low as 14 % in some species to as high as 98 % in others. We found that taxonomy is a major factor influencing annotation completeness, with distinct trends observed across the microbial tree (e.g. the lowest level of completeness was found in the Patescibacteria lineage). Most lineages showed a significant association between genome size and annotation incompleteness, likely reflecting a greater degree of uncharacterized sequences in 'accessory' proteomes than in 'core' proteomes. Finally, research bias, as measured by publication volume, was also an important factor influencing genome annotation completeness, with early model organisms showing high completeness levels relative to other genomes in their own taxonomic lineages. Our work highlights the disparity in annotation coverage across the bacterial tree of life and emphasizes a need for more experimental characterization of accessory proteomes as well as understudied lineages.

Exploring the Evolution of Virulence Factors through Bioinformatic Data Mining.

The molecular evolution of virulence factors is a central theme in our understanding of bacterial pathogenesis and host-microbe interactions. Using bioinformatics and genome data mining, recent studies have shed light on the evolution of important virulence factor families and the mechanisms by which they have adapted and diversified in function. This perspective highlights three complementary approaches useful for studying the molecular evolution of virulence factors: identification and analysis of virulence factor homologs, detection of adaptations or functional shifts, and computational prediction of novel virulence factor families. Each of these research directions is associated with distinct questions, approaches, and challenges for future work. Moving forward, bioinformatics will continue to play a critical role in exploring the evolution of virulence factors, including those that target humans. By reconstructing past processes and events, we will be able to better interpret newly sequenced microbial genomes and detect future pathoadaptations.

AnnoTree: visualization and exploration of a functionally annotated microbial tree of life.

Bacterial genomics has revolutionized our understanding of the microbial tree of life; however, mapping and visualizing the distribution of functional traits across bacteria remains a challenge. Here, we introduce AnnoTree-an interactive, functionally annotated bacterial tree of life that integrates taxonomic, phylogenetic and functional annotation data from over 27 000 bacterial and 1500 archaeal genomes. AnnoTree enables visualization of millions of precomputed genome annotations across the bacterial and archaeal phylogenies, thereby allowing users to explore gene distributions as well as patterns of gene gain and loss in prokaryotes. Using AnnoTree, we examined the phylogenomic distributions of 28 311 gene/protein families, and measured their phylogenetic conservation, patchiness, and lineage-specificity within bacteria. Our analyses revealed widespread phylogenetic patchiness among bacterial gene families, reflecting the dynamic evolution of prokaryotic genomes. Genes involved in phage infection/defense, mobile elements, and antibiotic resistance dominated the list of most patchy traits, as well as numerous intriguing metabolic enzymes that appear to have undergone frequent horizontal transfer. We anticipate that AnnoTree will be a valuable resource for exploring prokaryotic gene histories, and will act as a catalyst for biological and evolutionary hypothesis generation. AnnoTree is freely available at http://annotree.uwaterloo.ca.

The impact of cigarette smoke exposure, COPD, or asthma status on ABC transporter gene expression in human airway epithelial cells.

ABC transporters are conserved in prokaryotes and eukaryotes, with humans expressing 48 transporters divided into 7 classes (ABCA, ABCB, ABCC, ABCD, ABDE, ABCF, and ABCG). Throughout the human body, ABC transporters regulate cAMP levels, chloride secretion, lipid transport, and anti-oxidant responses. We used a bioinformatic approach complemented with in vitro experimental methods for validation of the 48 known human ABC transporters in airway epithelial cells using bronchial epithelial cell gene expression datasets available in NCBI GEO from well-characterized patient populations of healthy subjects and individuals that smoke cigarettes, or have been diagnosed with COPD or asthma, with validation performed in Calu-3 airway epithelial cells. Gene expression data demonstrate that ABC transporters are variably expressed in epithelial cells from different airway generations, regulated by cigarette smoke exposure (ABCA13, ABCB6, ABCC1, and ABCC3), and differentially expressed in individuals with COPD and asthma (ABCA13, ABCC1, ABCC2, ABCC9). An in vitro cell culture model of cigarette smoke exposure was able to recapitulate select observed in situ changes. Our work highlights select ABC transporter candidates of interest and a relevant in vitro model that will enable a deeper understanding of the contribution of ABC transporters in the respiratory mucosa in lung health and disease.

Metagenomic Sequencing of Wastewater from a South African Research Farm.

We sequenced wastewater effluent from the Agricultural Research Council-Animal Production in South Africa that conducts studies on livestock health and farm ecology. Thauera, Oscillibacter, and Pseudomonas were the most abundant genera within the community. Thirty-one different antibiotic resistance genes were identified, 10 of which are associated with tetracycline resistance.

Draft Genome Sequences of Two Novel Cellulolytic Streptomyces Strains Isolated from South African Rhizosphere Soil.

We report here the draft genome sequences of two novel strains of Streptomyces (NWU339 and NWU49) isolated from South African rhizosphere soils. Both strains were found to possess strong cellulolytic activity and contain numerous putative cellulase genes. Both genomes possess benzoate degradation pathways, while NWU49 contains the genomic potential for enediyne biosynthesis.

Novel function discovery through sequence and structural data mining.

Large-scale sequence and structural data is a goldmine of novel proteins, but how can this data be effectively mined for new functions? Here, we review protein function prediction methods and recent studies that apply these methods to discover new functionality. Core approaches include sequence-based homology detection, phylogenetic analysis, structural bioinformatics, and inference of functional associations using genomic context and related methods. With such a wide range of approaches, sequences may reveal new functionality regardless of their similarity to a characterized reference. Homologs of known function may be identified in unexpected species or associations. Detection of functional shifts in sequences may reveal new activities and specificities. New protein functions may also be predicted in uncharacterized sequences and structures. Finally, methods and data may be integrated and applied at increasingly large scales due to improved protein domain knowledge and structural coverage, which amplifies the ability to predict and discover novel protein functions.

MetAnnotate: function-specific taxonomic profiling and comparison of metagenomes.

BACKGROUND: Metagenomes provide access to the taxonomic composition and functional capabilities of microbial communities. Although metagenomic analysis methods exist for estimating overall community composition or metabolic potential, identifying specific taxa that encode specific functions or pathways of interest can be more challenging. Here we present MetAnnotate, which addresses the common question: "which organisms perform my function of interest within my metagenome(s) of interest?" MetAnnotate uses profile hidden Markov models to analyze shotgun metagenomes for genes and pathways of interest, classifies retrieved sequences either through a phylogenetic placement or best hit approach, and enables comparison of these profiles between metagenomes. RESULTS: Based on a simulated metagenome dataset, the tool achieves high taxonomic classification accuracy for a broad range of genes, including both markers of community abundance and specific biological pathways. Lastly, we demonstrate MetAnnotate by analyzing for cobalamin (vitamin B12) synthesis genes across hundreds of aquatic metagenomes in a fraction of the time required by the commonly used Basic Local Alignment Search Tool top hit approach. CONCLUSIONS: MetAnnotate is multi-threaded and installable as a local web application or command-line tool on Linux systems. Metannotate is a useful framework for general and/or function-specific taxonomic profiling and comparison of metagenomes.