Skip to main content

The taxonomic distribution of histamine-secreting bacteria in the human gut microbiome

Abstract

Background

Biogenic histamine plays an important role in immune response, neurotransmission, and allergic response. Although endogenous histamine production has been extensively studied, the contributions of histamine produced by the human gut microbiota have not been explored due to the absence of a systematic annotation of histamine-secreting bacteria.

Results

To identify the histamine-secreting bacteria from in the human gut microbiome, we conducted a systematic search for putative histamine-secreting bacteria in 36,554 genomes from the Genome Taxonomy Database and Unified Human Gastrointestinal Genome catalog. Using bioinformatic approaches, we identified 117 putative histamine-secreting bacteria species. A new three-component decarboxylation system including two colocalized decarboxylases and one transporter was observed in histamine-secreting bacteria among three different phyla. We found significant enrichment of histamine-secreting bacteria in patients with inflammatory bowel disease but not in patients with colorectal cancer suggesting a possible association between histamine-secreting bacteria and inflammatory bowel disease.

Conclusions

The findings of this study expand our knowledge of the taxonomic distribution of putative histamine-secreting bacteria in the human gut.

Peer Review reports

Background

Histamine is a health-relevant biogenic amine that plays important physiological roles in vascular permeability, mucus secretion, and neurotransmission via immunomodulation [1,2,3]. The enzyme histidine decarboxylase produces histamine via the decarboxylation of the amino acid histidine [4]. Two major families of histidine decarboxylase have been identified: pyridoxal-5′-phosphate (PLP) dependent histidine decarboxylases, which require PLP as a cofactor; and pyruvoyl-dependent histidine decarboxylases, which require a covalently bonded pyruvoyl moiety instead of PLP [5]. The decarboxylation of histidine occurs in the bacterial cytoplasm, therefore the histidine/histamine antiporter, which transports histidine into the cell and exports histamine out of the cell, is necessary for histidine decarboxylation to occur [6, 7].

While endogenously-produced histamine has been extensively studied, studies on exogenously-produced histamine have focused mostly on food-borne poisoning via histamine present in fish and dairy products [8, 9]. For example, food contaminated with a high concentration of histamine can cause neurological, gastrointestinal, and respiratory disorders [10, 11]. Histamine accumulated in food products is primarily derived from histamine-secreting bacteria, which generate histamine and other amines to maintain neutral cytoplasmic pH, allowing them to survive in acidic conditions [12]. For example, Lactobacillus vaginalis was found to produce histamine to maintain appropriate cytosolic pH in acidic conditions [13]. Pathogenic enteric bacteria also utilize amino acid decarboxylases to survive passage through the highly acidic gastric environment before reaching the gut. For example, the decarboxylation of L-arginine to agmatine offers Escherichia coli a robust acid-resistance mechanism [14, 15]. Alterations in microbiome taxonomic composition have been linked to many inflammatory diseases such as Inflammatory Bowel Disease (IBD) and asthma [16, 17]. For example, a recent study showed that histamine-secreting bacteria were increased in the gut of asthma patients [18]. Studies also demonstrated that the histamine secreted by the microbes influences immune responses within the gut via the host histamine receptor 2 and could exhibit anti-inflammatory effects [19,20,21].

However, our understanding of histamine-secreting bacteria and exogenous histamine production in the human gut is incomplete and limited to only a few cultured species [8, 22]. Despite the sequencing of numerous bacterial genomes, we lack a clear understanding of which gut bacteria are capable of producing and secreting histamine [23]. Therefore, a systematic functional annotation of histamine-secreting bacteria is needed to provide a comprehensive understanding of the abundance and prevalence of histamine-secreting bacteria in the human gut and to assess the role of bacterially-derived histamine in human health and disease. In this study, we conducted a systematic in-silico search of putative histamine-secreting bacteria in the 31,910 genomes from the Genome Taxonomy Database (GTDB) [24, 25] and the 4644 genomes in the Unified Human Gastrointestinal Genome database (UHGG) [26]. We further analyzed the relative abundance of putative histamine-secreting bacteria in metagenomic sequencing data from colorectal cancer (CRC) and inflammatory bowel disease (IBD) cohorts. Here, we aim to increase our understanding of the prevalence and abundance of histamine-secreting bacteria in the human gut microbiome and probe whether the abundance of histamine-secreting bacteria is altered in inflammatory disease.

Results

Histamine-secreting bacteria are sporadically distributed across six bacterial phyla

We manually collected and curated profile hidden markov models (HMM) for gene families involved in the histidine decarboxylation pathway and performed a systematic search on the proteomes from 31,910 genomes in GTDB. We identified 97 putative histamine-secreting bacteria with at least one histidine decarboxylase (hdcA) and one amino acid transporter (hdcP) localized in a gene cluster. Of these, 57 species contained the pyruvoyl-dependent histidine decarboxylase, 39 species contained the PLP-dependent histidine decarboxylase, and one species, Plesiomonas shigelloides, contained both pathways (Table S1). Plesiomonas shigelloides is a known histamine producer that causes acute diarrhea in humans after seafood consumption [27].

The putative histamine-secreting bacteria we identified were distributed across 6 phyla. Among the identified putative histamine-secreting bacteria, the species containing the pyruvoyl-dependent histidine decarboxylase were found in 6 phyla: Proteobacteria, Actinobacteriota, Bacteroidota, Firmicutes, Fusobacteriota, Verrucomicrobiota (Fig. 1A). Species containing the PLP-dependent histidine decarboxylase were found only in the phylum Proteobacteria (Fig. 1B). Proteobacteria expansion is considered to be a potential indicator of gut dysbiosis, a characteristic of proinflammatory diseases such as IBD [28]. Thirty-seven of the putative histamine-secreting bacteria we computationally identified in GTDB were experimentally confirmed to be histamine producers in prior studies (Table S2), validating the search strategy used in this study. Compared to previously identified species, we have identified additional putative histamine-secreting species in the Phyla Fusobacteriota and Verrucomicrobiota including Cetobacterium somerae, Fusobacterium ulcerans, Fusobacterium ulcerans A, Fusobacterium varium, Fusobacterium varium A, and 21–14–0-10-35-9 sp002773835 (in Verrucomicrobiota). To the best of our knowledge, no species in these two phyla has been experimentally verified to produce histamine.

Fig. 1
figure 1

Distribution of putative histamine-secreting bacteria across the Genome Taxonomy Database (GTDB) collection. Phylogenetic trees of (A) pyruvoyl-dependent and (B) PLP-dependent putative histamine-secreting bacteria across the 31,910 genomes of the representative GTDB collection were shown. The tip label colors represent different phyla. Stacked bar charts aligned to tree tips represent the percentages of genomes with and without HDC clusters in the species. The following red and blue numbers indicate the number of genomes with HDC present and absent, respectively

Among the putative histamine-secreting bacteria in GTDB, five genera including Bacteroides, Clostridium, Bifidobacterium, Fusobacterium, and Lactobacillus are common in the human gut microbiota [29,30,31,32]. Therefore, we searched for histamine-secreting bacteria in the UHGG database following the same procedure as in the GTDB database. Among the 4644 species in UHGG, we identified 44 putative histamine-secreting bacteria. Of these, 35 species encode a pyruvoyl-dependent decarboxylase, 8 species encode a PLP-dependent decarboxylase, and Plesiomonas shigelloides encodes both pyruvoyl-dependent and PLP-dependent decarboxylases (Fig. S1, Table S3). All phyla, 21 out of 24 genera, and 33 out of 44 species of putative histamine-secreting bacteria found in UHGG were also identified in GTDB. Among the 11 species exclusively found in UHGG, five are unassigned new species and unique in UHGG, while the other 6 species are likely due to the strain-level difference between the species collected in the GTDB and UHGG database.

We further examined the validity of the decarboxylases with structural modeling. We performed 3D protein homology modeling of pyruvoyl-dependent hdcA based on a high-resolution crystal structure (Lactobacillus saerimneri PDB: 1PYA) [33]. The backbone RMSD scores of pyruvoyl-dependent hdcA models were within 1 Å when comparing with the crystalline template (Fig. S2), indicating significant structural similarity between the putative histamine-secreting bacteria and the known histamine-secreting bacterium Lactobacillus saerimneri. Figure S2A shows the similarity between the model of hdcA of Clostridium perfringens and the template Lactobacillus saerimneri 30a (PDB: 1PYA) with a backbone RMSD of only 0.61 Å. The experimentally confirmed key residues in the active site of hdcA in Lactobacillus saerimneri 30a were illustrated in Fig. S2B including the Ser-81 and Ser-82 as the autocleavage pair [5], Ile-59 as the “lid” on the substrate-binding pocket [6], the Asp-198 and Asp-53 pair as the pH-regulating bridge [7], and the Tyr-62 and Asp-63 as the ligand-binding residues [33, 34]. The high structural similarity and identical key residues, as shown in Fig. S2C, strengthens the evidence of the genomic potential to produce histamine in these bacteria. We did not perform modeling on hdcP and PLP-dependent hdcA due to the lack of crystal structures and established enzymatic mechanisms for these proteins.

The ability to produce and secrete histamine is not a core function in all identified species. Out of the 97 species with a putative histamine-secretion operon in GTDB, 21 species belong to a clade where more than 50% of the strains lack the operon. These strain-specific histamine-secreting bacteria were predominantly observed in the phyla Firmicutes and Proteobacteria. Within the Proteobacteria, the genomic potential for histamine-secretion was found to be strain-specific in the families Pasteurellaceae, Vibrionaceae, and Enterobacteriaceae. In the Firmicutes, the genomic potential of histamine secretion was also found highly strain-specific in the Bacillaceae, Enterococcaceae, Lactobacillaceae, Staphylococcaceae, and Streptococcaceae families (Table S1). Species such as Streptococcus thermophilus, Staphylococcus warneri, Lactobacillus parabuchneri, and Lactobacillus reuteri have been previously reported as strain-specific in terms of histamine production [35,36,37,38], which agrees with the results of our computational search.

The strain-specificity of histamine-secretion may be due to the frequent loss of the histamine-secretion genes during evolution. Strain-specific gene loss is a common reason for the functional variation between different strains of the same species [39,40,41,42]. The strain-specificity could also be attributed to the horizontal gene transfer of genes encoding histamine-secretion pathways by mobile genetic elements. For example, while the histidine decarboxylase gene cluster is located on the chromosome of the Lactobacillus reuteri JCM 1112 strain, it is located on the pLRI01 plasmid of the Lactobacillus reuteri I5007 strain. The histidine decarboxylase gene cluster on the plasmid can be mobilized via conjugation and might have been horizontally transferred to Clostridium thermobutyricum (GCF_000371465.1). In Clostridium thermobutyricum, the histidine decarboxylase cluster is both phylogenetically closer (Fig. S3) and syntenically more similar (Fig. 2) to the histidine decarboxylase cluster in the plasmid of Lactobacillaceae than that in other Clostridiaceae. Such mobility of histidine decarboxylase gene clusters was also reported in a study on the species of Lactobacillus parabuchneri [44], supporting the hypothesis that mobile genetic elements could play an important role in the strain-specificity of histamine-secreting traits.

Fig. 2
figure 2

Representative histidine decarboxylase gene clusters for putative histamine-secreting bacteria. Two types of bacterial histidine decarboxylases were labeled as hdcA pyruvoyl and hdcA PLP respectively. Two arginine decarboxylases were identified in some histidine decarboxylase species including aaxB and adiA. Three different types of amino acid antiporters in histidine decarboxylase clusters (hdcP) were labelled as gadC, aaxC and adiC. The phylogenetic tree was a subtree extracted from a pre-built GTDB phylogenomic tree based on 120 bacterial marker genes. Gene clusters along with the tree were visualized by ggtree (version 2.5.1) [43]. The break symbol (double slash mark) was placed between two distant histidine decarboxylase gene clusters for species Plesiomonas shigelloides

New gene clusters observed in the putative histamine-secreting bacteria

Three different types of antiporters (hdcP) were found in the histidine decarboxylase gene clusters. These antiporters were annotated as the glutamate/gamma-aminobutyrate antiporter “gadC”, and the arginine/agmatine antiporters “aaxC” and “adiC”. The gadC and aaxC genes were observed with the pyruvoyl-dependent hdcA whereas adiC was found with the PLP-dependent hdcA. The only exception was seen in Plesiomonas shigelloides, which has both the pyruvoyl-dependent and PLP-dependent pathways, where the adiC antiporter was found in the gene cluster of both the PLP-dependent and pyruvoyl-dependent hdcA. At the phylum level, in Bacteroidota, Firmicutes, and Fusobacteriota, histidine decarboxylase was adjacent to aaxC; in Verrucomicrobiota, histidine decarboxylase was adjacent to gadC; and in Actinobacteroiodota, histidine decarboxylase was adjacent to either aaxC or gadC (Fig. S3).

Duplications of the antiporters were observed in some Proteobacteria where antiporter genes were found both upstream and downstream of the PLP-dependent hdcA. In addition to hdcA and hdcP, hdcB, which catalyzes the maturation of pyruvoyl-dependent hdcA, was occasionally found in close proximity to hdcA as seen in Streptococcus thermophilus CHCC1524 [45]. The genes, hisRS, which encode histidyl-tRNA synthetase were also found in the histidine decarboxylase cluster [46,47,48]. Throughout our search, the hdcA/antiporter/hdcB cluster was observed in Firmicutes such as Clostridium thermobutyricum and Streptococcus thermophilus. The hdcA/hdcB/antiporter/hisRS cluster was also observed in Firmicutes such as Tetragenococcus halophilus and Lactobacillus reuteri. The hdcA/antiporter/hisRS cluster was observed in Proteobacteria such as Klebsiella aerogenes and Shewanella marina.

In several phyla of the pyruvoyl-dependent histamine-secreting bacteria, we found a two-decarboxylase, one-transporter three-component decarboxylation system (Fig. 2) composed of arginine decarboxylase, histidine decarboxylase, and an antiporter. This arginine/histidine system was similar to the lysine/ornithine three-component decarboxylation system in Lactobacillus saerimneri 30a [49], the first reported bacterial three-component decarboxylation system. A difference between the arginine/histidine system found in histamine-secreting bacteria and the lysine/ornithine system found in Lactobacillus saerimneri 30a is that gene components of the histidine/arginine system were adjacent to each other, while the genes for the lysine/ornithine system were 24 kb apart.

The two decarboxylases in the three-component histidine/arginine system can utilize different enzymatic mechanisms and cofactors. Similar to histidine decarboxylase, two types of arginine decarboxylases have been previously identified, a PLP-dependent arginine decarboxylase (adiAYC) and a pyruvoyl-dependent arginine decarboxylase (aaxABC) [50]. In the histidine/arginine systems, the histidine decarboxylases were pyruvoyl-dependent while the arginine decarboxylases were either pyruvoyl-dependent (aaxB) or PLP-dependent (adiA). For example, in Actinobacteriota and Fusobacteriota, the adiA (PLP-dependent)/antiporter/hdcA (pyruvoyl-dependent) cluster was observed, while in Bacteroidota, both the adiA/antiporter/hdcA (Ornithobacterium rhinotracheale) and aaxB/antiporter/hdcA (Bacteroides stercorirosoris) clusters were observed.

Putative histamine-secreting bacteria were significantly enriched in inflammatory bowel disease patients but not in colorectal cancer patients compared to healthy controls

Analyzing 2451 stool metagenomes from IBD patients and 506 stool metagenomes from CRC patients, we found the putative histamine-secreting bacteria identified in UHGG were significantly enriched in IBD (Ulcerative colitis and Crohn’s disease combined) patients but not in CRC patients. We identified UHGG species with differential abundance between the patients and healthy controls in eight studies, including four IBD studies and four CRC studies. To avoid bias caused by detection methods, we used three methods, namely DESeq2, MaAsLin2, and LEfSe to perform the differential abundance analysis. We then performed one-tailed two-proportion Z-tests with continuity correction to test the hypothesis that histamine secreting species have a higher probability to be positively associated with disease. In the IBD studies, we found that putative histamine-secreting bacteria were significantly enriched when compared to all gut species through Z-tests with p < 0.05 in all four studies and by all three statistical methods (Fig. 3, Table S4-S5). The results from DESeq2, MaAsLin2, and LEfSe generally agreed well with each other, with DESeq2 being the most relaxed method in determining enrichment. Unlike IBD studies, no significant differential abundance was observed in putative histamine-secreting bacteria when compared to the overall gut species in any CRC studies by any method.

Fig. 3
figure 3

Differential abundance analysis of putative histamine-secreting bacteria. The analysis was performed on stool metagenomic samples between patients with colorectal cancer (CRC), inflammatory bowel disease (IBD) and healthy controls by DESeq2 (A), MaAslin2 (B) and LefSe (C). Relative abundances (counts) were calculated using Kraken2 (2.0.8-beta) (see methods section for more details). The criteria for enriched species for DEseq2 are log2FC > 1 (the Integrative Human Microbiome Project, i.e. HMP2 > 0.5) and q < 0.05, for MaAsLin2 are coef > 0.2 (HMP2 > 0) and q < 0.05, and for LefSe are LDA > 2 and q < 0.001. The statistical differences between all bacteria and histamine-secreting bacteria (HSB) were determined by one-tailed two-proportion Z-test with continuity correction and the p values were converted to asterisks (n.s. for p > 0.05; * for p ≤ 0.05; ** for p ≤ 0.01; *** for p ≤ 0.001 and **** for p ≤ 0.0001). Number of samples by disease state for each study are: Feng 2015 [51]: control = 61, CRC = 46; Vogtmann 2016 [52]: control = 52, CRC = 52; Yu 2017 [53]: control = 75, CRC = 53; Zeller 2014 [54]: control = 66, CRC = 91; Hall 2017 [55]: control = 74, IBD = 188; Franzosa 2019 (discovery cohorts) [56]: control = 34, IBD = 121; HMP2 [57, 58]: control = 429, IBD = 1209; Nielsen 0214 [59]: control = 248, IBD = 148. See also Table S4 for the enrichment of individual species and Table S5 for the detailed results of Z-test

The enrichment of histamine-secreting species was not attributed to a single taxon and we found the enriched species in IBD studies to be highly cohort-specific. For example, in the cohort in Hall et al. 2017, histamine-secreting bacteria were mainly enriched in Actinobacteriota and slightly enriched in Firmicutes. For Franzosa et al. 2019, histamine-secreting bacteria were mainly enriched in Firmicutes and Proteobacteria (Fig. 3, Table S4). For HMP2, histamine-secreting bacteria were mainly enriched in Bacteroidiota and Proteobacteria and for Nielsen et al. 2014, histamine-secreting bacteria were enriched in Actinobacteriota (Fig. 3, Table S4). In addition, we analyzed the abundance of histidine decarboxylase operon in IBD and CRC studies. Consistent with the results of differential abundance analysis of putative histamine-secreting bacteria, a similar trend was found that HDC gene clusters were enriched in IBD studies. Our result indicated the operon abundance was significantly enriched in all of the IBD studies including Franzosa et al. 2019, Hall et al. 2017, HMP2 and Nielson et al. 2014, and one of the CRC studies, namely Yu et al. 2017 (Fig. S4).

Discussion

In this study, we performed a systematic search for histamine secreting operons in species from the human gut microbiota. In total, we identified 117 species with the genomic potential to secrete histamine. We identified a few discrepancies in the histamine-secreting bacteria between our systematic search and previous experimental studies (Table S2). Several reasons could contribute to the discrepancies. First, the putative histamine-secreting operon was highly strain-specific in a number of histamine-secreting species [36]. The difference could be due to different strains used in this study and previous work. Second, we focused on the well-known pyruvoyl-dependent and PLP-dependent histidine decarboxylases and this approach might miss newly discovered histidine decarboxylases [60]. What’s more, the techniques to detect histamine produced by histamine-secreting bacteria varied between studies [8, 61]. For example, in Citrobacter freundii, 21 out of 21 strains were classified as positive for histamine-secretion by Niven’s agar method but only 1 out of 21 were classified as positive for histamine-secretion by PCR [61]. The Niven’s agar is a differential growth medium with a pH indicator. Bacteria that grow and increase the pH of the Niven’s agar are inferred to be as histamine-secreting bacteria. This method has a high false-positive rate largely due to the fact that bacteria can secrete alkaline other than histamine. The PCR based method is subject to primer bias which can lead to false negative results. The possible inaccuracy in the histamine detection methods has likely led to non-histamine-secreting bacteria mislabeled as histamine-secreting and vice versa. Further studies may include experimental evaluation of selected representatives of the putative histamine-secreting bacteria identified in this study, especially at some health-related taxonomic levels where no histamine-secreting bacteria have been described and explored before.

In this study, we identified that putative histamine-secreting bacteria were significantly enriched in IBD patients. Combined with the previous knowledge that histamine-secreting bacteria were found to be increased in the gut communities of asthma patients [18], it suggests a possible association between histamine-secreting bacteria and inflammatory diseases. Because the lumenal environment is more acidic in the small intestine than in the colon, the genomic potential of histamine-secretion may be expressed in the small intestine but not in the colon. As a result, the presence of putative histamine-secreting bacteria in prominent families found in the small intestines, such as Enterobacteriaceae and Lactobacillaceae, may have significant immunological implications (Table S3) [62]. The deeper relationships between histamine-secreting bacteria and host in terms of inflammatory and immunological diseases are yet to be elucidated.

To the best of our knowledge, the histidine/arginine decarboxylation system observed in a number of histamine-secreting bacteria from three phyla is the first three-component decarboxylation system identified in a cluster of colocalized genes. Our findings suggest that some amino acid antiporters are likely capable of transporting multiple amino acids with similar structural and chemical properties. The capability to transport multiple amino acids using a single antiporter may increase the robustness and efficiency of acid resistance in the bacteria. Further studies may include the exploration of other three-component or multi-component decarboxylation systems other than the histidine/arginine and lysine/ornithine systems.

Conclusions

In conclusion, we have systematically annotated bacterial histidine decarboxylase in two large databases. 117 putative histamine-secreting bacteria species were identified throughout GTDB and UHGG including those in Fusobacteriota and Verrucomicrobiota where no histamine-secreting bacteria had been previously described. We have identified a novel three-component decarboxylase system which contains an arginine decarboxylase, a histidine decarboxylase and an antiporter. Differential abundance analysis showed that the putative histamine-secreting bacteria were significantly enriched in IBD patients but not in CRC patients. Our findings would expand the knowledge and provide a comprehensive understanding of histamine-secreting bacteria in the human gut and facilitate advances in potential therapeutic targets toward histamine-related inflammatory and immunological diseases.

Methods

Histidine decarboxylase operon identification

The Unified Human Gastrointestinal Genome (UHGG) [26] and Genome Taxonomy Database (GTDB) (release 95) were used as the databases for histamine-secreting bacteria identification [24, 25]. To identify PLP-dependent histidine decarboxylases, we performed a multiple sequence alignment of known PLP-dependent histidine decarboxylase genes (WP_191935110.1,WP_068969528.1,WP_152135723.1 and WP_136342781.1) with MUSCLE 3.8.31 [63] and then built a customized profile HMM with hmmbuild from HMMER 3.3.1 package [64]. This customized profile was used to search query protein databases using hmmscan. Hits with e-value 1e-100 or less were selected and manually curated and annotated as PLP-dependent histidine decarboxylases. To identify pyruvoyl-dependent histidine decarboxylases, we searched the query database using hmmscan against Pfam profile PF02329 (Histidine carboxylase PI chain) with cutoff e-value 1e-40. We searched for amino acid antiporters using hmmscan against PANTHER profile PTHR42770 with cutoff e-value 1e-40. The genomic loci where histidine decarboxylase and antiporters are present in the same gene neighborhood (less than three genes away) were identified as histamine secretion gene clusters. The gene structure in putative histamine-secreting bacteria was then manually screened in Geneious Prime 2021.0.3 [65].

Histidine decarboxylase 3D modeling

For the putative pyruvoyl-dependent decarboxylase hdcA, the amino acid sequences of histamine-secreting bacteria from UHGG and GTDB were aligned with Clustal Omega [66, 67]. BLAST and HHpred servers were used to identify potential structural templates from the Protein Data Bank [68, 69]. We used the 2.5 Å X-ray crystal structure of a bacterial HhdA from Lactobacillus saerimneri 30a (UniProt ID P00862, PDB entry 1PYA) as the sole template [33]. Due to the inter-oligomeric contacts and active site, a trimer C3 symmetry was defined based on 1PYA and used in modeling. The 3mer and 9mer fragment files were obtained from the Robetta server (http://old.robetta.org/). RosettaCM was then employed to generate at least 1000 models of each protein. We selected the top ten models based on the Rosetta energy and, among these ten, selected the lowest backbone RMSD to the 1PYA trimer crystal structure as the final protein models. Unlike pyruvoyl-dependent hdcA with auto-cleavage and multi-chain modeling with only one good template, no crystal structure was reported for bacterial hdcP and the PLP-dependent hdcA therefore no model was built upon hdcP and the PLP-dependent hdcA.

Phylogenetic analysis

To construct the phylogeny of pyruvoyl-dependent hdcA genes, a multiple sequence alignment was performed by Clustal Omega v1.2.3 with fifty-eight pyruvoyl hdcA proteins [66] and trimmed with trimAl v1.2 on strictplus mode [70]. A maximum-likelihood phylogenetic tree was inferred by IQ-TREE v2.1.2 using its suggested LG + I + G4 model with 1000 ultrafast bootstrap replicates [71]. This tree was rooted using the Minimal Ancestor Deviation (MAD) method via mad v2.2 [72] and was visualized and annotated using iTOL v5 (https://itol.embl.de/) [73].

Metagenomic data processing

The raw sequencing reads of metagenomic samples used in this study were downloaded and extracted using National Center for Biotechnology Information (NCBI)‘s SRA Toolkit [74] v2.10.9 under the accession numbers of PRJEB1220 [59], HMP2 (PRJNA398089 [58], PRJNA389280 [57]), PRJNA385949 [55], and PRJNA400072 [56] for IBD studies and accession numbers for studies are PRJEB6070 [54], PRJEB7774 [51], PRJEB12449 [52], and PRJEB10878 [53] for colorectal cancer (CRC) studies. Quality control and adapter trimming of the fastq sequence files were done with Trim Galore wrapper v0.6.6 [75]. Quality trimmed sequences were screened against the human assembly NCBI build 37 (hg19) using Bowtie2 v2.4.2 alignment software [76] and Samtools v1.11 [77] was then used to remove human genome contamination unmapped sequence from SAM files. Taxonomic assignment of filtered reads was performed using Kraken2 v2.0.8-beta (with default settings) [78] against a pre-built database of the UHGG catalog (http://ftp.ebi.ac.uk/pub/databases/metagenomics/mgnify_genomes/human-gut/v1.0/uhgg_kraken2-db/) [26].

Differential taxonomic abundance analysis

The kraken species-level abundance outputs and metadata were imported into the phyloseq R package (v.1.34.0) for analyses [79]. Samples with sequencing depth less than 1 million counts were excluded. Rare taxa with a relative abundance of less than 0.01% across 10% of all samples were filtered. To minimize the bias of a single method and report a robust analysis, we employed three commonly used statistical methods including DESeq2 [80], Multivariate Association with Linear Models 2 (MaAsLin2) [81], and Linear discriminant analysis effect size (LEfSe) [82] for differential abundance analysis. For DEseq2, the “local” method was employed to estimate dispersion, and the “poscounts” size factor estimator was employed in the normalization step to exclude the zeros when calculating the geometric mean. The DEseq2 “poscounts” normalized counts were used as input for LEfSe. For MaAsLin2, the trimmed mean of M-values (TMM) method was applied as the normalization method based on its satisfactory performance in a recent benchmark [83]. We defined species with more than 50% of the strains encoding histidine decarboxylase gene clusters as histamine-secreting species. One-tailed, two-proportion Z-tests with continuity correction were performed to compare the difference between the proportion of disease-associated species in defined histamine secreting species and the proportion of disease-associated species in all species.

Histidine decarboxylase operon abundance estimation and analysis

We created a set of reference genomes mainly from the representative genomes from the UHGG. If the representative genome for a species in the UHGG does not contain the histidine decarboxylase operon and a different genome of the same species has, we replaced it with the genome with the histidine decarboxylase operon instead. We aligned the cleaned metagenomic reads to the reference genomes using bowtie2 [76]. Samples with total read counts less than 1 million were filtered out. The abundance of histidine decarboxylase operon was estimated by the number of per million reads mapped to the histidine decarboxylase operon divided by the total number reads mapped to the reference genomes (namely measured in counts per million reads mapped, CPM). One sided Wilcoxon’s rank-sum test with continuity correction was performed to test the difference of the histidine decarboxylase operon abundance in the patient’s sample and that in the healthy control.

Availability of data and materials

Genomes analysed during the current study are available in the GTDB, https://data.ace.uq.edu.au/public/gtdb/data/releases/release95/95.0/ and UHGG, http://ftp.ebi.ac.uk/pub/databases/metagenomics/mgnify_genomes/human-gut/v1.0/uhgg_catalogue/. The raw sequencing reads of metagenomic samples were downloaded from Sequence Read Archive (SRA) from the National Center for Biotechnology Information (NCBI) database (SRA accession numbers of PRJEB1220, PRJNA389280, PRJNA398089, PRJNA385949, and PRJNA400072 for IBD studies and PRJEB6070, PRJEB7774, PRJEB12449, and PRJEB10878 for CRC studies). All data generated or analyzed during this study are included in this published article (and its supplementary information files).

Abbreviations

CRC:

Colorectal cancer

GTDB:

Genome Taxonomy Database

HDC:

Histidine decarboxylase

HMP2:

Integrative Human Microbiome Project

HSB:

Histamine-secreting bacteria

IBD:

Inflammatory bowel disease

PDB:

Protein data bank

PLP:

Pyridoxal-5′-phosphate

RMSD:

Root mean square deviation

UHGG:

Unified Human Gastrointestinal Genome

References

  1. Parsons ME, Ganellin CR. Histamine and its receptors. Br J Pharmacol. 2006;147(Suppl 1):S127–35. https://doi.org/10.1038/sj.bjp.0706440.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  2. White MV. The role of histamine in allergic diseases. J Allergy Clin Immunol. 1990;86(4 Pt 2):599–605. https://doi.org/10.1016/S0091-6749(05)80223-4.

    Article  PubMed  CAS  Google Scholar 

  3. O’Mahony L, Akdis M, Akdis CA. Regulation of the immune response and inflammation by histamine and histamine receptors. J Allergy Clin Immunol. 2011;128(6):1153–62. https://doi.org/10.1016/j.jaci.2011.06.051.

    Article  PubMed  CAS  Google Scholar 

  4. Shahid M, Tripathi T, Sobia F, Moin S, Siddiqui M, Khan RA. Histamine, histamine receptors, and their role in immunomodulation: an updated systematic review. Open Immunol J. 2009;2(1):9–41. https://doi.org/10.2174/1874226200902010009.

    Article  CAS  Google Scholar 

  5. Gallagher T, Snell EE, Hackert ML. Pyruvoyl-dependent histidine decarboxylase: active site structure and mechanistic analysis. J Biol Chem. 1989;264(21):12737–43. https://doi.org/10.1016/S0021-9258(18)63917-1.

    Article  PubMed  CAS  Google Scholar 

  6. Pishko EJ, Robertus’ JD. Site-Directed Alteration of Three Active-Site Residues of a Pyruvoyl-Dependent Histidine Decarboxylase+. Biochemistry. 1993;32:4943–8.

  7. Pishko EJ, Potter KA, Robertus JD. Site-directed mutagenesis of intersubunit boundary residues in histidine decarboxylase, a pH-dependent allosteric enzyme. Biochemistry. 1995;34(18):6069–73. https://doi.org/10.1021/bi00018a009.

    Article  PubMed  CAS  Google Scholar 

  8. Moniente M, García-Gonzalo D, Ontañón I, Pagán R, Botello-Morte L. Histamine accumulation in dairy products: microbial causes, techniques for the detection of histamine-producing microbiota, and potential solutions. Compr Rev Food Sci Food Saf. 2021;20(2):1481–523. https://doi.org/10.1111/1541-4337.12704.

    Article  PubMed  CAS  Google Scholar 

  9. Visciano P, Schirone M, Paparella A. An overview of histamine and other biogenic amines in fish and fish products. Foods. 2020;9(12):1795. https://doi.org/10.3390/foods9121795.

    Article  PubMed Central  CAS  Google Scholar 

  10. Alvarez MA, Moreno-Arribas MV. The problem of biogenic amines in fermented foods and the use of potential biogenic amine-degrading microorganisms as a solution. Trends Food Sci Technol. 2014;39(2):146–55. https://doi.org/10.1016/j.tifs.2014.07.007.

    Article  CAS  Google Scholar 

  11. Ladero V, Calles-Enriquez M, Fernandez M. A. Alvarez M. toxicological effects of dietary biogenic amines. Curr Nutr Food Sci. 2010;6(2):145–56. https://doi.org/10.2174/157340110791233256.

    Article  CAS  Google Scholar 

  12. Krulwich TA, Sachs G, Padan E. Molecular aspects of bacterial pH sensing and homeostasis. Nat Rev Microbiol. 2011;9(5):330–43. https://doi.org/10.1038/nrmicro2549.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  13. Diaz M, Del Rio B, Ladero V, Redruello B, Fernández M, Martin MC, et al. Histamine production in Lactobacillus vaginalis improves cell survival at low pH by counteracting the acidification of the cytosol. Int J Food Microbiol. 2020;321:108548. https://doi.org/10.1016/j.ijfoodmicro.2020.108548.

    Article  PubMed  CAS  Google Scholar 

  14. Foster JW. Escherichia coli acid resistance: tales of an amateur acidophile. Nat Rev Microbiol. 2004;2(11):898–907. https://doi.org/10.1038/nrmicro1021.

    Article  PubMed  CAS  Google Scholar 

  15. Richard H, Foster JW. Escherichia coli glutamate- and arginine-dependent acid resistance systems increase internal pH and reverse transmembrane potential. J Bacteriol. 2004;186(18):6032–41. https://doi.org/10.1128/JB.186.18.6032-6041.2004.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  16. Kostic AD, Xavier RJ, Gevers D. The microbiome in inflammatory bowel disease: current status and the future ahead. Gastroenterology. 2014;146(6):1489–99. https://doi.org/10.1053/j.gastro.2014.02.009.

    Article  PubMed  CAS  Google Scholar 

  17. Huang YJ, Boushey HA. The microbiome in asthma. J Allergy Clin Immunol. 2015;135(1):25–30. https://doi.org/10.1016/j.jaci.2014.11.011.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Barcik W, Pugin B, Westermann P, Perez NR, Ferstl R, Wawrzyniak M, et al. Histamine-secreting microbes are increased in the gut of adult asthma patients. J Allergy Clin Immunol. 2016;138(5):1491–4. https://doi.org/10.1016/j.jaci.2016.05.049.

    Article  PubMed  Google Scholar 

  19. Ferstl R, Frei R, Schiavi E, Konieczna P, Barcik W, Ziegler M, et al. Histamine receptor 2 is a key influence in immune responses to intestinal histamine-secreting microbes. J Allergy Clin Immunol. 2014;134(3):744–6. https://doi.org/10.1016/j.jaci.2014.04.034.

    Article  PubMed  CAS  Google Scholar 

  20. Shi Z, Fultz RS, Engevik MA, Gao C, Hall A, Major A, et al. Distinct roles of histamine H1- and H2-receptor signaling pathways in inflammation-associated colonic tumorigenesis. Am J Physiol Gastrointest Liver Physiol. 2019;316(1):G205–16. https://doi.org/10.1152/ajpgi.00212.2018.

    Article  PubMed  CAS  Google Scholar 

  21. Gao C, Ganesh BP, Shi Z, Shah RR, Fultz R, Major A, et al. Gut microbe-mediated suppression of inflammation-associated Colon carcinogenesis by luminal histamine production. Am J Pathol. 2017;187(10):2323–36. https://doi.org/10.1016/j.ajpath.2017.06.011.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  22. Barcik W, Wawrzyniak M, Akdis CA, O’Mahony L. Immune regulation by histamine and histamine-secreting bacteria. Curr Opin Immunol. 2017;48:108–13. https://doi.org/10.1016/j.coi.2017.08.011.

    Article  PubMed  CAS  Google Scholar 

  23. Tian L, Wang X-W, Wu A-K, Fan Y, Friedman J, Dahlin A, et al. Deciphering functional redundancy in the human microbiome. Nat Commun. 2020;11(1):6217. https://doi.org/10.1038/s41467-020-19940-1.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  24. Parks DH, Chuvochina M, Waite DW, Rinke C, Skarshewski A, Chaumeil P-A, et al. A standardized bacterial taxonomy based on genome phylogeny substantially revises the tree of life. Nat Biotechnol. 2018;36(10):996–1004. https://doi.org/10.1038/nbt.4229.

    Article  PubMed  CAS  Google Scholar 

  25. Parks DH, Chuvochina M, Chaumeil P-A, Rinke C, Mussig AJ, Hugenholtz P. A complete domain-to-species taxonomy for Bacteria and Archaea. Nat Biotechnol. 2020;38(9):1079–86. https://doi.org/10.1038/s41587-020-0501-8.

    Article  PubMed  CAS  Google Scholar 

  26. Almeida A, Nayfach S, Boland M, Strozzi F, Beracochea M, Shi ZJ, et al. A unified catalog of 204,938 reference genomes from the human gut microbiome. Nat Biotechnol. 2021;39(1):105–14. https://doi.org/10.1038/s41587-020-0603-3.

    Article  PubMed  CAS  Google Scholar 

  27. López-Sabater EI, Rodríguez-Jerez JJ, Hernández-Herrero M, Mora-Ventura MT. Incidence of histamine-forming bacteria and histamine content in scombroid fish species from retail markets in the Barcelona area. Int J Food Microbiol. 1996;28(3):411–8. https://doi.org/10.1016/0168-1605(94)00007-7.

    Article  PubMed  Google Scholar 

  28. Mukhopadhya I, Hansen R, El-Omar EM, Hold GL. IBD-what role do Proteobacteria play? Nat Rev Gastroenterol Hepatol. 2012;9(4):219–30. https://doi.org/10.1038/nrgastro.2012.14.

    Article  PubMed  CAS  Google Scholar 

  29. Zhang J, Guo Z, Xue Z, Sun Z, Zhang M, Wang L, et al. A phylo-functional core of gut microbiota in healthy young Chinese cohorts across lifestyles, geography and ethnicities. ISME J. 2015;9(9):1979–90. https://doi.org/10.1038/ismej.2015.11.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Falony G, Joossens M, Vieira-Silva S, Wang J, Darzi Y, Faust K, et al. Population-level analysis of gut microbiome variation. Science. 2016;352(6285):560–4. https://doi.org/10.1126/science.aad3503.

    Article  PubMed  CAS  Google Scholar 

  31. Guarner F, Malagelada J-R. Gut flora in health and disease. Lancet. 2003;361(9356):512–9. https://doi.org/10.1016/S0140-6736(03)12489-0.

    Article  PubMed  Google Scholar 

  32. Forster SC, Kumar N, Anonye BO, Almeida A, Viciani E, Stares MD, et al. A human gut bacterial genome and culture collection for improved metagenomic analyses. Nat Biotechnol. 2019;37(2):186–92. https://doi.org/10.1038/s41587-018-0009-7.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  33. Gallagher T, Rozwarski DA, Ernst SR, Hackert ML. Refined structure of the pyruvoyl-dependent histidine decarboxylase from Lactobacillus 30a. J Mol Biol. 1993;230(2):516–28. https://doi.org/10.1006/jmbi.1993.1168.

    Article  PubMed  CAS  Google Scholar 

  34. Schelp E, Worley S, Monzingo AF, Ernst S, Robertus JD. pH-induced structural changes regulate histidine decarboxylase activity in Lactobacillus 30a. J Mol Biol. 2001;306(4):727–32. https://doi.org/10.1006/jmbi.2000.4430.

    Article  PubMed  CAS  Google Scholar 

  35. Calles-Enríquez M, Eriksen BH, Andersen PS, Rattray FP, Johansen AH, Fernández M, et al. Sequencing and transcriptional analysis of the Streptococcus thermophilus histamine biosynthesis gene cluster: factors that affect differential hdcA expression. Appl Environ Microbiol. 2010;76(18):6231–8. https://doi.org/10.1128/AEM.00827-10.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  36. Economou V, Gousia P, Kemenetzi D, Sakkas H, Papadopoulou C. Microbial quality and histamine producing microflora analysis of the ice used for fish preservation: microbial quality of ice used to preserve fish. J Food Saf. 2017;37(1):e12285. https://doi.org/10.1111/jfs.12285.

    Article  CAS  Google Scholar 

  37. Mu Q, Tavella VJ, Luo XM. Role of Lactobacillus reuteri in human health and diseases. Front Microbiol. 2018;9:757. https://doi.org/10.3389/fmicb.2018.00757.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Gao C, Major A, Rendon D, Lugo M, Jackson V, Shi Z, et al. Histamine H2 receptor-mediated suppression of intestinal inflammation by probiotic Lactobacillus reuteri. MBio. 2015;6(6):e01358–15. https://doi.org/10.1128/mBio.01358-15.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  39. Simmons SS, Isokpehi RD, Brown SD, McAllister DL, Hall CC, McDuffy WM, et al. Functional annotation analytics of Rhodopseudomonas palustris genomes. Bioinform Biol Insights. 2011;5:115–29. https://doi.org/10.4137/BBI.S7316.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  40. Oda Y, Larimer FW, Chain PSG, Malfatti S, Shin MV, Vergez LM, et al. Multiple genome sequences reveal adaptations of a phototrophic bacterium to sediment microenvironments. Proc Natl Acad Sci U S A. 2008;105(47):18543–8. https://doi.org/10.1073/pnas.0809160105.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Ceapa C, Davids M, Ritari J, Lambert J, Wels M, Douillard FP, et al. The variable regions of Lactobacillus rhamnosus genomes reveal the dynamic evolution of metabolic and host-adaptation repertoires. Genome Biol Evol. 2016;8(6):1889–905. https://doi.org/10.1093/gbe/evw123.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  42. Marri PR, Hao W, Golding GB. Gene gain and gene loss in streptococcus: is it driven by habitat? Mol Biol Evol. 2006;23(12):2379–91. https://doi.org/10.1093/molbev/msl115.

    Article  PubMed  CAS  Google Scholar 

  43. Yu G. Using ggtree to visualize data on tree-like structures. Curr Protoc Bioinformatics. 2020;69(1):e96. https://doi.org/10.1002/cpbi.96.

    Article  PubMed  Google Scholar 

  44. Wüthrich D, Berthoud H, Wechsler D, Eugster E, Irmler S, Bruggmann R. The histidine decarboxylase gene cluster of Lactobacillus parabuchneri was gained by horizontal gene transfer and is Mobile within the species. Front Microbiol. 2017;8:218. https://doi.org/10.3389/fmicb.2017.00218.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Trip H, Mulder NL, Rattray FP, Lolkema JS. HdcB, a novel enzyme catalysing maturation of pyruvoyl-dependent histidine decarboxylase. Mol Microbiol. 2011;79(4):861–71. https://doi.org/10.1111/j.1365-2958.2010.07492.x.

    Article  PubMed  CAS  Google Scholar 

  46. Lucas PM, Wolken WAM, Claisse O, Lolkema JS, Lonvaud-Funel A. Histamine-producing pathway encoded on an unstable plasmid in Lactobacillus hilgardii 0006. Appl Environ Microbiol. 2005;71(3):1417–24. https://doi.org/10.1128/AEM.71.3.1417-1424.2005.

  47. Martín MC, Fernández M, Linares DM, Alvarez MA. Sequencing, characterization and transcriptional analysis of the histidine decarboxylase operon of Lactobacillus buchneri. Microbiology. 2005;151(Pt 4):1219–28. https://doi.org/10.1099/mic.0.27459-0.

    Article  PubMed  CAS  Google Scholar 

  48. Satomi M, Furushita M, Oikawa H, Yoshikawa-Takahashi M, Yano Y. Analysis of a 30 kbp plasmid encoding histidine decarboxylase gene in Tetragenococcus halophilus isolated from fish sauce. Int J Food Microbiol. 2008;126(1-2):202–9. https://doi.org/10.1016/j.ijfoodmicro.2008.05.025.

    Article  PubMed  CAS  Google Scholar 

  49. Romano A, Trip H, Lolkema JS, Lucas PM. Three-component lysine/ornithine decarboxylation system in Lactobacillus saerimneri 30a. J Bacteriol. 2013;195(6):1249–54. https://doi.org/10.1128/JB.02070-12.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  50. Smith CB, Graham DE. Outer and inner membrane proteins compose an arginine-agmatine exchange system in Chlamydophila pneumoniae. J Bacteriol. 2008;190(22):7431–40. https://doi.org/10.1128/JB.00652-08.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  51. Feng Q, Liang S, Jia H, Stadlmayr A, Tang L, Lan Z, et al. Gut microbiome development along the colorectal adenoma-carcinoma sequence. Nat Commun. 2015;6(1):6528. https://doi.org/10.1038/ncomms7528.

    Article  PubMed  CAS  Google Scholar 

  52. Vogtmann E, Hua X, Zeller G, Sunagawa S, Voigt AY, Hercog R, et al. Colorectal Cancer and the human gut microbiome: reproducibility with whole-genome shotgun sequencing. PLoS One. 2016;11(5):e0155362. https://doi.org/10.1371/journal.pone.0155362.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  53. Yu J, Feng Q, Wong SH, Zhang D, Liang QY, Qin Y, et al. Metagenomic analysis of faecal microbiome as a tool towards targeted non-invasive biomarkers for colorectal cancer. Gut. 2017;66(1):70–8. https://doi.org/10.1136/gutjnl-2015-309800.

    Article  PubMed  CAS  Google Scholar 

  54. Zeller G, Tap J, Voigt AY, Sunagawa S, Kultima JR, Costea PI, et al. Potential of fecal microbiota for early-stage detection of colorectal cancer. Mol Syst Biol. 2014;10(11):766. https://doi.org/10.15252/msb.20145645.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  55. Hall AB, Yassour M, Sauk J, Garner A, Jiang X, Arthur T, et al. A novel Ruminococcus gnavus clade enriched in inflammatory bowel disease patients. Genome Med. 2017;9(1):103. https://doi.org/10.1186/s13073-017-0490-5.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  56. Franzosa EA, Sirota-Madi A, Avila-Pacheco J, Fornelos N, Haiser HJ, Reinker S, et al. Gut microbiome structure and metabolic activity in inflammatory bowel disease. Nat Microbiol. 2019;4(2):293–305. https://doi.org/10.1038/s41564-018-0306-4.

    Article  PubMed  CAS  Google Scholar 

  57. Schirmer M, Franzosa EA, Lloyd-Price J, McIver LJ, Schwager R, Poon TW, et al. Dynamics of metatranscription in the inflammatory bowel disease gut microbiome. Nat Microbiol. 2018;3(3):337–46. https://doi.org/10.1038/s41564-017-0089-z.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  58. Lloyd-Price J, Arze C, Ananthakrishnan AN, Schirmer M, Avila-Pacheco J, Poon TW, et al. Multi-omics of the gut microbial ecosystem in inflammatory bowel diseases. Nature. 2019;569(7758):655–62. https://doi.org/10.1038/s41586-019-1237-9.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  59. Nielsen HB, Almeida M, Juncker AS, Rasmussen S, Li J, Sunagawa S, et al. Identification and assembly of genomes and genetic elements in complex metagenomic samples without using reference genomes. Nat Biotechnol. 2014;32(8):822–8. https://doi.org/10.1038/nbt.2939.

    Article  PubMed  CAS  Google Scholar 

  60. Bjornsdottir-Butler K, May S, Hayes M, Abraham A, Benner RA Jr. Characterization of a novel enzyme from Photobacterium phosphoreum with histidine decarboxylase activity. Int J Food Microbiol. 2020;334:108815. https://doi.org/10.1016/j.ijfoodmicro.2020.108815.

    Article  PubMed  CAS  Google Scholar 

  61. Bjornsdottir K, Bolton GE, McClellan-Green PD, Jaykus L-A, Green DP. Detection of gram-negative histamine-producing bacteria in fish: a comparative study. J Food Prot. 2009;72(9):1987–91. https://doi.org/10.4315/0362-028X-72.9.1987.

    Article  PubMed  CAS  Google Scholar 

  62. Donaldson GP, Lee SM, Mazmanian SK. Gut biogeography of the bacterial microbiota. Nat Rev Microbiol. 2016;14(1):20–32. https://doi.org/10.1038/nrmicro3552.

    Article  PubMed  CAS  Google Scholar 

  63. Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32(5):1792–7. https://doi.org/10.1093/nar/gkh340.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  64. Mistry J, Finn RD, Eddy SR, Bateman A, Punta M. Challenges in homology search: HMMER3 and convergent evolution of coiled-coil regions. Nucleic Acids Res. 2013;41(12):e121–1. https://doi.org/10.1093/nar/gkt263.

  65. Geneious. 2019. https://www.geneious.com/. Accessed 7 Apr 2021.

  66. Sievers F, Wilm A, Dineen D, Gibson TJ, Karplus K, Li W, et al. Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal omega. Mol Syst Biol. 2011;7(1):539. https://doi.org/10.1038/msb.2011.75.

    Article  PubMed  PubMed Central  Google Scholar 

  67. Sievers F, Higgins DG. Clustal omega for making accurate alignments of many protein sequences. Protein Sci. 2018;27(1):135–45. https://doi.org/10.1002/pro.3290.

    Article  PubMed  CAS  Google Scholar 

  68. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215(3):403–10. https://doi.org/10.1016/S0022-2836(05)80360-2.

    Article  PubMed  CAS  Google Scholar 

  69. Söding J, Biegert A, Lupas AN. The HHpred interactive server for protein homology detection and structure prediction. Nucleic Acids Res. 2005;33(Suppl 2):W244–8. https://doi.org/10.1093/nar/gki408.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  70. Capella-Gutiérrez S, Silla-Martínez JM, Gabaldón T. trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics. 2009;25(15):1972–3. https://doi.org/10.1093/bioinformatics/btp348.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  71. Nguyen L-T, Schmidt HA, von Haeseler A, Minh BQ. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol. 2015;32(1):268–74. https://doi.org/10.1093/molbev/msu300.

    Article  CAS  PubMed  Google Scholar 

  72. Tria FDK, Landan G, Dagan T. Phylogenetic rooting using minimal ancestor deviation. Nat Ecol Evol. 2017;1(1):193. https://doi.org/10.1038/s41559-017-0193.

    Article  PubMed  Google Scholar 

  73. Letunic I, Bork P. Interactive tree of life (iTOL) v5: an online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 2021;49(W1):W293–6. https://doi.org/10.1093/nar/gkab301.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  74. sra-tools. Github. https://github.com/ncbi/sra-tools. Accessed 7 Apr 2021.

  75. Babraham Bioinformatics - Trim Galore! https://www.bioinformatics.babraham.ac.uk/projects/trim_galore/. Accessed 14 Apr 2021.

  76. Langmead B, Salzberg SL. Fast gapped-read alignment with bowtie 2. Nat Methods. 2012;9(4):357–9. https://doi.org/10.1038/nmeth.1923.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  77. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25(16):2078–9. https://doi.org/10.1093/bioinformatics/btp352.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  78. Wood DE, Lu J, Langmead B. Improved metagenomic analysis with kraken 2. Genome Biol. 2019;20(1):257. https://doi.org/10.1186/s13059-019-1891-0.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  79. McMurdie PJ, Holmes S. phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS ONE. 2013;8:e61217.

  80. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550. https://doi.org/10.1186/s13059-014-0550-8.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  81. Mallick H, Rahnavard A, McIver LJ, Ma S, Zhang Y, Nguyen LH, et al. Multivariable Association Discovery in Population-scale Meta-omics Studies. Cold Spring Harbor Laboratory. 2021;2021.01.20.427420. https://doi.org/10.1101/2021.01.20.427420.

  82. Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, et al. Metagenomic biomarker discovery and explanation. Genome Biol. 2011;12(6):R60. https://doi.org/10.1186/gb-2011-12-6-r60.

    Article  PubMed  PubMed Central  Google Scholar 

  83. Pereira MB, Wallroth M, Jonsson V, Kristiansson E. Comparison of normalization methods for the analysis of metagenomic gene abundance data. BMC Genomics. 2018;19(1):274. https://doi.org/10.1186/s12864-018-4637-6.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

Download references

Acknowledgments

We would like to thank Keith Dufault-Thompson for his helpful feedback on a draft of the paper.

Funding

Z.M., Y.Y, and X.J. are supported by the Intramural Research Program of the NIH, National Library of Medicine. B.H. is supported by startup funding from the University of Maryland. Open Access funding provided by the National Institutes of Health (NIH).

Author information

Authors and Affiliations

Authors

Contributions

X.J. conceived and designed the study. Z.M. performed genomic and metagenomic data analysis and prepared Fig. 3. Y.Y performed gene cluster abundance analysis and phylogenetic analysis, and prepared Figs. 1 and 2. Z.M., Y.Y., and X.J. wrote the draft manuscript. X.J. and B.H. reviewed and revised the manuscript. X.J. supervised the work. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Xiaofang Jiang.

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.

Additional information

Publisher’s Note

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

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Mou, Z., Yang, Y., Hall, A.B. et al. The taxonomic distribution of histamine-secreting bacteria in the human gut microbiome. BMC Genomics 22, 695 (2021). https://doi.org/10.1186/s12864-021-08004-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-021-08004-3

Keywords