The taxonomic distribution of histamine-secreting bacteria in the human gut microbiome
BMC Genomics volume 22, Article number: 695 (2021)
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.
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.
The findings of this study expand our knowledge of the taxonomic distribution of putative histamine-secreting bacteria in the human gut.
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 . 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 . 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 . For example, Lactobacillus vaginalis was found to produce histamine to maintain appropriate cytosolic pH in acidic conditions . 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 . 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 . 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) . 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.
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 .
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 . 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.
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) . 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 , Ile-59 as the “lid” on the substrate-binding pocket , the Asp-198 and Asp-53 pair as the pH-regulating bridge , 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 , supporting the hypothesis that mobile genetic elements could play an important role in the strain-specificity of histamine-secreting traits.
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 . 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 , 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) . 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.
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).
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 . 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 . 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 . 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 , 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) . 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.
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.
Histidine decarboxylase operon identification
The Unified Human Gastrointestinal Genome (UHGG)  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  and then built a customized profile HMM with hmmbuild from HMMER 3.3.1 package . 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 .
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 . 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.
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  and trimmed with trimAl v1.2 on strictplus mode . 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 . This tree was rooted using the Minimal Ancestor Deviation (MAD) method via mad v2.2  and was visualized and annotated using iTOL v5 (https://itol.embl.de/) .
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  v2.10.9 under the accession numbers of PRJEB1220 , HMP2 (PRJNA398089 , PRJNA389280 ), PRJNA385949 , and PRJNA400072  for IBD studies and accession numbers for studies are PRJEB6070 , PRJEB7774 , PRJEB12449 , and PRJEB10878  for colorectal cancer (CRC) studies. Quality control and adapter trimming of the fastq sequence files were done with Trim Galore wrapper v0.6.6 . Quality trimmed sequences were screened against the human assembly NCBI build 37 (hg19) using Bowtie2 v2.4.2 alignment software  and Samtools v1.11  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)  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/) .
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 . 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 , Multivariate Association with Linear Models 2 (MaAsLin2) , and Linear discriminant analysis effect size (LEfSe)  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 . 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 . 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).
Genome Taxonomy Database
Integrative Human Microbiome Project
Inflammatory bowel disease
Protein data bank
Root mean square deviation
Unified Human Gastrointestinal Genome
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.
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.
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.
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.
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.
Pishko EJ, Robertus’ JD. Site-Directed Alteration of Three Active-Site Residues of a Pyruvoyl-Dependent Histidine Decarboxylase+. Biochemistry. 1993;32:4943–8.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Geneious. 2019. https://www.geneious.com/. Accessed 7 Apr 2021.
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.
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.
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.
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.
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.
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.
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.
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.
sra-tools. Github. https://github.com/ncbi/sra-tools. Accessed 7 Apr 2021.
Babraham Bioinformatics - Trim Galore! https://www.bioinformatics.babraham.ac.uk/projects/trim_galore/. Accessed 14 Apr 2021.
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.
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.
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.
McMurdie PJ, Holmes S. phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS ONE. 2013;8:e61217.
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.
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.
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.
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.
We would like to thank Keith Dufault-Thompson for his helpful feedback on a draft of the paper.
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).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 4 Table S4. Enriched HSB of UHGG in IBD studies. This table corresponds to the species in
Fig. 3of the main text. (XLSX 13 kb)
Additional file 5 Table S5. One-tailed two-proportion Z tests for enriched HSB of UHGG in IBD and CRC studies. This table corresponds to the enriched proportion in
Fig. 3of the main text. (XLSX 11 kb)
About this article
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
- Histidine decarboxylase
- Histamine-secreting bacteria
- Human gut microbiota
- Inflammatory bowel disease