Pd(II)-Catalyzed 1,2-Oxyarylation of Alkenes with O-Acylhydroxylamines as the Oxygen Source.

O-Acylhydroxylamine has been widely employed as an electrophilic amination reagent in transition-metal-catalyzed C-N coupling reactions, but its use as an electrophilic oxygen source has not been disclosed. Here, we report a Pd-catalyzed 1,2-oxyarylation of alkenes with O-acylhydroxylamines as an oxidant and an oxygen source for the first time. With simple amide as the monodentate directing group, this method features a broad substrate scope, good functional group tolerance, and mild conditions.

Host genetic regulation of human gut microbial structural variation.

Although the impact of host genetics on gut microbial diversity and the abundance of specific taxa is well established1-6, little is known about how host genetics regulates the genetic diversity of gut microorganisms. Here we conducted a meta-analysis of associations between human genetic variation and gut microbial structural variation in 9,015 individuals from four Dutch cohorts. Strikingly, the presence rate of a structural variation segment in Faecalibacterium prausnitzii that harbours an N-acetylgalactosamine (GalNAc) utilization gene cluster is higher in individuals who secrete the type A oligosaccharide antigen terminating in GalNAc, a feature that is jointly determined by human ABO and FUT2 genotypes, and we could replicate this association in a Tanzanian cohort. In vitro experiments demonstrated that GalNAc can be used as the sole carbohydrate source for F. prausnitzii strains that carry the GalNAc-metabolizing pathway. Further in silico and in vitro studies demonstrated that other ABO-associated species can also utilize GalNAc, particularly Collinsella aerofaciens. The GalNAc utilization genes are also associated with the host's cardiometabolic health, particularly in individuals with mucosal A-antigen. Together, the findings of our study demonstrate that genetic associations across the human genome and bacterial metagenome can provide functional insights into the reciprocal host-microbiome relationship.

Diphosphine Ligand-Enabled Nickel-Catalyzed Chelate-Assisted Inner-Selective Migratory Hydroarylation of Alkenes.

The precise control of the regioselectivity in the transition metal-catalyzed migratory hydrofunctionalization of alkenes remains a big challenge. With a transient ketimine directing group, the nickel-catalyzed migratory ß-selective hydroarylation and hydroalkenylation of alkenyl ketones has been realized with aryl boronic acids using alkyl halide as the mild hydride source for the first time. The key to this success is the use of a diphosphine ligand, which is capable of the generation of a Ni(II)-H species in the presence of alkyl bromide, and enabling the efficient migratory insertion of alkene into Ni(II)-H species and the sequent rapid chain walking process. The present approach diminishes organosilanes reductant, tolerates a wide array of complex functionalities with excellent regioselective control. Moreover, this catalytic system could also be applied to the migratory hydroarylation of alkenyl azahetereoarenes, thus providing a general approach for the preparation of 1,2-aryl heteroaryl motifs with wide potential applications in pharmaceutical discovery.

Gut microbial structural variation associates with immune checkpoint inhibitor response.

The gut microbiota may have an effect on the therapeutic resistance and toxicity of immune checkpoint inhibitors (ICIs). However, the associations between the highly variable genomes of gut bacteria and the effectiveness of ICIs remain unclear, despite the fact that merely a few gene mutations between similar bacterial strains may cause significant phenotypic variations. Here, using datasets from the gut microbiome of 996 patients from seven clinical trials, we systematically identify microbial genomic structural variants (SVs) using SGV-Finder. The associations between SVs and response, progression-free survival, overall survival, and immune-related adverse events are systematically explored by metagenome-wide association analysis and replicated in different cohorts. Associated SVs are located in multiple species, including Akkermansia muciniphila, Dorea formicigenerans, and Bacteroides caccae. We find genes that encode enzymes that participate in glucose metabolism be harbored in these associated regions. This work uncovers a nascent layer of gut microbiome heterogeneity that is correlated with hosts' prognosis following ICI treatment and represents an advance in our knowledge of the intricate relationships between microbiota and tumor immunotherapy.

Prevotella copri alleviates sarcopenia via attenuating muscle mass loss and function decline.

BACKGROUND: The gut microbiome and fecal metabolites have been found to influence sarcopenia, but whether there are potential bacteria that can alleviate sarcopenia has been under-investigated, and the molecular mechanism remains unclear. METHODS: To investigate the relationships between the gut microbiome, fecal metabolites and sarcopenia, subjects were selected from observational multi-ethnic study conducted in Western China. Sarcopenia was diagnosed according to the criteria of the Asian Working Group for Sarcopenia 2014. The gut microbiome was profiled by shotgun metagenomic sequencing. Untargeted metabolomic analysis was performed to analyse the differences in fecal metabolites. We investigated bacterium with the greatest relative abundance difference between healthy individuals and sarcopenia patients, and the differences in metabolites associated with the bacteria, to verify its effects on muscle mass and function in a mouse model. RESULTS: The study included 283 participants (68.90% females, mean age: 66.66 years old) with and without sarcopenia (141 and 142 participants, respectively) and from the Han (98 participants), Zang (88 participants) and Qiang (97 participants) ethnic groups. This showed an overall reduction (15.03% vs. 20.77%, P = 0.01) of Prevotella copri between the sarcopenia and non-sarcopenia subjects across the three ethnic groups. Functional characterization of the differential bacteria showed enrichment (odds ratio = 15.97, P = 0.0068) in branched chain amino acid (BCAA) metabolism in non-sarcopenia group. A total of 13 BCAA and their derivatives have relatively low levels in sarcopenia. In the in vivo experiment, we found that the blood BCAA level was higher in the mice gavaged with live P. copri (LPC) (P < 0.001). The LPC mice had significantly longer wire and grid hanging time (P < 0.02), longer time on rotor (P = 0.0001) and larger grip strength (P < 0.0001), indicating better muscle function. The weight of gastrocnemius mass and rectus femoris mass (P < 0.05) was higher in LPC mice. The micro-computed tomography showed a larger leg area (P = 0.0031), and a small animal analyser showed a higher lean mass ratio in LPC mice (P = 0.0157), indicating higher muscle mass. CONCLUSIONS: The results indicated that there were lower levels of both P. copri and BCAA in sarcopenia individuals. In vivo experiments, gavage with LPC could attenuate muscle mass and function decline, indicating alleviating sarcopenia. This suggested that P. copri may play a therapeutic potential role in the management of sarcopenia.

Gut dysbiosis associates with cytokine production capacity in viral-suppressed people living with HIV.

Background: People living with human immunodeficiency virus (PLHIV) are exposed to chronic immune dysregulation, even when virus replication is suppressed by antiretroviral therapy (ART). Given the emerging role of the gut microbiome in immunity, we hypothesized that the gut microbiome may be related to the cytokine production capacity of PLHIV. Methods: To test this hypothesis, we collected metagenomic data from 143 ART-treated PLHIV and assessed the ex vivo production capacity of eight different cytokines [interleukin-1ß (IL-1ß), IL-6, IL-1Ra, IL-10, IL-17, IL-22, tumor necrosis factor, and interferon-γ] in response to different stimuli. We also characterized CD4+ T-cell counts, HIV reservoir, and other clinical parameters. Results: Compared with 190 age- and sex-matched controls and a second independent control cohort, PLHIV showed microbial dysbiosis that was correlated with viral reservoir levels (CD4+ T-cell-associated HIV-1 DNA), cytokine production capacity, and sexual behavior. Notably, we identified two genetically different P. copri strains that were enriched in either PLHIV or healthy controls. The control-related strain showed a stronger negative association with cytokine production capacity than the PLHIV-related strain, particularly for Pam3Cys-incuded IL-6 and IL-10 production. The control-related strain is also positively associated with CD4+ T-cell level. Conclusions: Our findings suggest that modulating the gut microbiome may be a strategy to modulate immune response in PLHIV.

Ligand-Enabled NiII -Catalyzed Hydroxylarylation of Alkenes with Molecular Oxygen.

The use of molecular oxygen as the terminal oxidant in transition metal catalyzed oxidative process is an appealing and challenging task in organic synthetic chemistry. Here, we report a Ni-catalyzed hydroxylarylation of unactivated alkenes enabled by a ß-diketone ligand with high efficiency and excellent regioselectivity employing molecular oxygen as the oxidant and hydroxyl source. This reaction features mild conditions, broad substrate scope and incredible heterocycle compatibility, providing a variety of ß-hydroxylamides, γ-hydroxylamides, ß-aminoalcohols, γ-aminoalcohols, and 1,3-diols in high yields. The synthetic value of this methodology was demonstrated by the efficient synthesis of two bioactive compounds, (±)-3'-methoxyl citreochlorol and tea catechin metabolites M4.

CRISPR dynamics during the interaction between bacteria and phage in the first year of life.

Gut microbiomes in infancy have a profound impact on health in adulthood. CRISPRs play an essential role in the interaction between bacteria and phages. However, the dynamics of CRISPRs in gut microbiomes during early life are poorly understood. In this study, using shotgun metagenomic sequencing data from 82 Swedish infants' gut microbiomes, 1882 candidate CRISPRs were identified, and their dynamics were analysed. We found large-scale turnover of CRISPRs and their spacers during the first year of life. As well as changes in relative abundance of the bacteria containing CRISPR, acquisition, loss and mutation of spacers were observed within the same CRISPR array sampled over time. Accordingly, the inferred interaction network of bacteria and phage was distinct at different times. This research underpins CRISPR dynamics and their potential role in the interaction between bacteria and phage in early life.

Systematic analysis of relationships between plasma branched-chain amino acid concentrations and cardiometabolic parameters: an association and Mendelian randomization study.

BACKGROUND: Branched-chain amino acids (BCAAs; valine, leucine, and isoleucine) are essential amino acids that are associated with an increased risk of cardiometabolic diseases (CMD). However, there are still only limited insights into potential direct associations between BCAAs and a wide range of CMD parameters, especially those remaining after correcting for covariates and underlying causal relationships. METHODS: To shed light on these relationships, we systematically characterized the associations between plasma BCAA concentrations and a large panel of 537 CMD parameters (including atherosclerosis-related parameters, fat distribution, plasma cytokine concentrations and cell counts, circulating concentrations of cardiovascular-related proteins and plasma metabolites) in 1400 individuals from the Dutch population cohort LifeLines DEEP and 294 overweight individuals from the 300OB cohort. After correcting for age, sex, and BMI, we assessed associations between individual BCAAs and CMD parameters. We further assessed the underlying causality using Mendelian randomization. RESULTS: A total of 838 significant associations were detected for 409 CMD parameters. BCAAs showed both common and specific associations, with the most specific associations being detected for isoleucine. Further, we found that obesity status substantially affected the strength and direction of associations for valine, which cannot be corrected for using BMI as a covariate. Subsequent univariable Mendelian randomization (UVMR), after removing BMI-associated SNPs, identified seven significant causal relationships from four CMD traits to BCAA levels, mostly for diabetes-related parameters. However, no causal effects of BCAAs on CMD parameters were supported. CONCLUSIONS: Our cross-sectional association study reports a large number of associations between BCAAs and CMD parameters. Our results highlight some specific associations for isoleucine, as well as obesity-specific effects for valine. MR-based causality analysis suggests that altered BCAA levels can be a consequence of diabetes and alteration in lipid metabolism. We found no MR evidence to support a causal role for BCAAs in CMD. These findings provide evidence to (re)evaluate the clinical importance of individual BCAAs in CMD diagnosis, prevention, and treatment.

Microbial vitamin production mediates dietary effects on diabetic risk.

Adequate levels of essential vitamins are important for the prevention of diabetes. While the main efforts to address this are currently focused on the intake of vitamin supplements, improving and maintaining intrinsic vitamin production capacity, which is determined by gut microbes, has received insufficient attention. In this study, we systematically investigated the relationship between gut microbial vitamin production and factors related to diabetes and cardiometabolic health in a deeply phenotyped cohort, Lifelines-DEEP (N = 1,135). We found that blood glucose-related factors, lipids, circulating inflammation, and fecal short-chain fatty acids are associated with gut microbial vitamin production. Use of laxatives and metformin are associated with increased levels of vitamin B1/B6 biosynthesis pathways. We further reveal a mediatory role for microbial vitamin B1/B2 production on the influence of fruit intake on diabetes risk. This study provides preliminary evidence for microbiome-targeted vitamin metabolism interventions to promote health.

Ligand-enabled Ni-catalyzed hydroarylation and hydroalkenylation of internal alkenes with organoborons.

The transition metal-catalyzed hydrofunctionalization of alkenes offers an efficient solution for the rapid construction of complex functional molecules, and significant progress has been made during last decades. However, the hydrofunctionalization of internal alkenes remains a significant challenge due to low reactivity and the difficulties of controlling the regioselectivity. Here, we report the hydroarylation and hydroalkenylation of internal alkenes lacking a directing group with aryl and alkenyl boronic acids in the presence of a nickel catalyst, featuring a broad substrate scope and wide functional group tolerance under redox-neutral conditions. The key to achieving this reaction is the identification of a bulky 1-adamantyl ß-diketone ligand, which is capable of overcoming the low reactivity of internal 1,2-disubstituted alkenes. Preliminary mechanistic studies unveiled that this reaction undergoes an Ar-Ni(II)-H initiated hydroarylation process, which is generated by the oxidative addition of alcoholic solvent with Ni(0) species and sequential transmetalation. In addition, the oxidative addition of the alcoholic solvent proves to be the turnover-limiting step.

Plasticity of the adult human small intestinal stoma microbiota.

The human distal small intestine (ileum) has a distinct microbiota, but human studies investigating its composition and function have been limited by the inaccessibility of the ileum without purging and/or deep intubation. We investigated inherent instability, temporal dynamics, and the contribution of fed and fasted states using stoma samples from cured colorectal cancer patients as a non-invasive access route to the otherwise inaccessible small and large intestines. Sequential sampling of the ileum before and after stoma formation indicated that ileostoma microbiotas represented that of the intact small intestine. Ileal and colonic stoma microbiotas were confirmed as distinct, and two types of instability in ileal host-microbial relationships were observed: inter-digestive purging followed by the rapid postprandial blooming of bacterial biomass and sub-strain appearance and disappearance within individual taxa after feeding. In contrast to the relative stability of colonic microbiota, the human small intestinal microbiota biomass and its sub-strain composition can be highly dynamic.

Influence of the microbiome, diet and genetics on inter-individual variation in the human plasma metabolome.

The levels of the thousands of metabolites in the human plasma metabolome are strongly influenced by an individual's genetics and the composition of their diet and gut microbiome. Here, by assessing 1,183 plasma metabolites in 1,368 extensively phenotyped individuals from the Lifelines DEEP and Genome of the Netherlands cohorts, we quantified the proportion of inter-individual variation in the plasma metabolome explained by different factors, characterizing 610, 85 and 38 metabolites as dominantly associated with diet, the gut microbiome and genetics, respectively. Moreover, a diet quality score derived from metabolite levels was significantly associated with diet quality, as assessed by a detailed food frequency questionnaire. Through Mendelian randomization and mediation analyses, we revealed putative causal relationships between diet, the gut microbiome and metabolites. For example, Mendelian randomization analyses support a potential causal effect of Eubacterium rectale in decreasing plasma levels of hydrogen sulfite-a toxin that affects cardiovascular function. Lastly, based on analysis of the plasma metabolome of 311 individuals at two time points separated by 4 years, we observed a positive correlation between the stability of metabolite levels and the amount of variance in the levels of that metabolite that could be explained in our analysis. Altogether, characterization of factors that explain inter-individual variation in the plasma metabolome can help design approaches for modulating diet or the gut microbiome to shape a healthy metabolome.

Helicobacter pylori infection affects the human gastric microbiome, as revealed by metagenomic sequencing.

Helicobacter pylori infection is a prevalent infectious disease, associated with many gastric diseases, including gastritis, gastric ulcer, and gastric cancer. To reveal the characteristics of the gastric microbiome in patients infected with H. pylori, we performed metagenomic shotgun sequencing of stomach swab samples from 96 patients and then conducted metagenomic association analyses between alterations in the gastric microbiome and H. pylori infection status. The overall composition of the gastric microbiota in H. pylori-infected individuals was distinctly different from the negative controls; H. pylori became the dominant species after colonizing the human stomach and significantly decreased the α-diversity of the gastric community (P < 0.05, Wilcoxon rank-sum test). We also identified 6 HPI-associated microbial species (FDR < 0.05, Wilcoxon rank-sum test): Stenotrophomonas maltophilia, Stenotrophomonas unclassified, Chryseobacterium unclassified, Pedobacter unclassified, Variovorax unclassified, and Pseudomonas stutzeri. Furthermore, 55 gastric microbial pathways were enriched in the H. pylori-positive group, whereas only 2 pathways were more abundant in the H. pylori-negative group: dTDP-L-rhamnose biosynthesis and tetrapyrrole biosynthesis (FDR < 0.05, Wilcoxon rank-sum test). Gastritis was not associated with non-H. pylori species in the stomach (P > 0.05, Wilcoxon rank-sum test). This study revealed alterations in gastric microbial taxa and function associated with HPI in the Chinese population, which provides an insight into gastric microbial interactions and their potential role in the pathological process of gastric diseases.

Discovery, diversity, and functional associations of crAss-like phages in human gut metagenomes from four Dutch cohorts.

The crAss-like phages are a diverse group of related viruses that includes some of the most abundant viruses of the human gut. To explore their diversity and functional role in human population and clinical cohorts, we analyze gut metagenomic data collected from 1,950 individuals from the Netherlands. We identify 1,556 crAss-like phage genomes, including 125 species-level and 32 genus-level clusters absent from the reference databases used. Analysis of their genomic features shows that closely related crAss-like phages can possess strikingly divergent regions responsible for transcription, presumably acquired through recombination. Prediction of crAss-like phage hosts points primarily to bacteria of the phylum Bacteroidetes, consistent with previous reports. Finally, we explore the temporal stability of crAss-like phages over a 4-year period and identify associations between the abundance of crAss-like phages and several human phenotypes, including depletion of crAss-like phages in inflammatory bowel disease patients.

Chromosome-Level Genome Assembly and HazelOmics Database Construction Provides Insights Into Unsaturated Fatty Acid Synthesis and Cold Resistance in Hazelnut (Corylus heterophylla).

Corylus heterophylla (2n = 22) is the most widely distributed, unique, and economically important nut species in China. Chromosome-level genomes of C. avellana, C. heterophylla, and C. mandshurica have been published in 2021, but a satisfactory hazelnut genome database is absent. Northeast China is the main distribution and cultivation area of C. heterophylla, and the mechanism underlying the adaptation of C. heterophylla to extremely low temperature in this area remains unclear. Using single-molecule real-time sequencing and the chromosomal conformational capture (Hi-C) assisted genome assembly strategy, we obtained a high-quality chromosome-scale genome sequence of C. heterophylla, with a total length of 343 Mb and scaffold N50 of 32.88 Mb. A total of 94.72% of the test genes from the assembled genome could be aligned to the Embryophyta_odb9 database. In total, 22,319 protein-coding genes were predicted, and 21,056 (94.34%) were annotated in the assembled genome. A HazelOmics online database (HOD) containing the assembled genome, gene-coding sequences, protein sequences, and various types of annotation information was constructed. This database has a user-friendly and straightforward interface. In total, 439 contracted genes and 3,810 expanded genes were identified through genome evolution analysis, and 17 expanded genes were significantly enriched in the unsaturated fatty acid biosynthesis pathway (ko01040). Transcriptome analysis results showed that FAD (Cor0058010.1), SAD (Cor0141290.1), and KAT (Cor0122500.1) with high expression abundance were upregulated at the ovule maturity stage. We deduced that the expansion of these genes may promote high unsaturated fatty acid content in the kernels and improve the adaptability of C. heterophylla to the cold climate of Northeast China. The reference genome and database will be beneficial for future molecular breeding and gene function studies in this nut species, as well as for evolutionary research on species of the order Fagales.

Characterization of gut microbial structural variations as determinants of human bile acid metabolism.

Bile acids (BAs) facilitate intestinal fat absorption and act as important signaling molecules in host-gut microbiota crosstalk. BA-metabolizing pathways in the microbial community have been identified, but it remains largely unknown how the highly variable genomes of gut bacteria interact with host BA metabolism. We characterized 8,282 structural variants (SVs) of 55 bacterial species in the gut microbiomes of 1,437 individuals from two cohorts and performed a systematic association study with 39 plasma BA parameters. Both variations in SV-based continuous genetic makeup and discrete clusters showed correlations with BA metabolism. Metagenome-wide association analysis identified 809 replicable associations between bacterial SVs and BAs and SV regulators that mediate the effects of lifestyle factors on BA metabolism. This is the largest microbial genetic association analysis to demonstrate the impact of bacterial SVs on human BA composition, and it highlights the potential of targeting gut microbiota to regulate BA metabolism through lifestyle intervention.

Clinical and CT Radiomics Nomogram for Preoperative Differentiation of Pulmonary Adenocarcinoma From Tuberculoma in Solitary Solid Nodule.

AIM: To investigate clinical and computed tomography (CT) radiomics nomogram for preoperative differentiation of lung adenocarcinoma (LAC) from lung tuberculoma (LTB) in patients with pulmonary solitary solid nodule (PSSN). MATERIALS AND METHODS: A total of 313 patients were recruited in this retrospective study, including 96 pathologically confirmed LAC and 217 clinically confirmed LTB. Patients were assigned at random to training set (n = 220) and validation set (n = 93) according to 7:3 ratio. A total of 2,589 radiomics features were extracted from each three-dimensional (3D) lung nodule on thin-slice CT images and radiomics signatures were built using the least absolute shrinkage and selection operator (LASSO) logistic regression. The predictive nomogram was established based on radiomics and clinical features. Decision curve analysis was performed with training and validation sets to assess the clinical usefulness of the prediction model. RESULTS: A total of six clinical features were selected as independent predictors, including spiculated sign, vacuole, minimum diameter of nodule, mediastinal lymphadenectasis, sex, and age. The radiomics nomogram of lung nodules, consisting of 15 selected radiomics parameters and six clinical features showed good prediction in the training set [area under the curve (AUC), 1.00; 95% confidence interval (CI), 0.99-1.00] and validation set (AUC, 0.99; 95% CI, 0.98-1.00). The nomogram model that combined radiomics and clinical features was better than both single models (p < 0.05). Decision curve analysis showed that radiomics features were beneficial to clinical settings. CONCLUSION: The radiomics nomogram, derived from unenhanced thin-slice chest CT images, showed favorable prediction efficacy for differentiating LAC from LTB in patients with PSSN.

Glycosaminoglycan biosynthesis pathway in host genome is associated with Helicobacter pylori infection.

Helicobacter pylori is a causative pathogen of many gastric and extra-gastric diseases. It has infected about half of the global population. There were no genome-wide association studies (GWAS) for H. pylori infection conducted in Chinese population, who carried different and relatively homogenous strain of H. pylori. In this work, we performed SNP (single nucleotide polymorphism)-based, gene-based and pathway-based genome-wide association analyses to investigate the genetic basis of host susceptibility to H. pylori infection in 480 Chinese individuals. We also profiled the composition and function of the gut microbiota between H. pylori infection cases and controls. We found several genes and pathways associated with H. pylori infection (P < 0.05), replicated one previously reported SNP rs10004195 in TLR1 gene region (P = 0.02). We also found that glycosaminoglycan biosynthesis related pathway was associated with both onset and progression of H. pylori infection. In the gut microbiome association study, we identified 2 species, 3 genera and several pathways had differential abundance between H. pylori infected cases and controls. This paper is the first GWAS for H. pylori infection in Chinese population, and we combined the genetic and microbial data to comprehensively discuss the basis of host susceptibility to H. pylori infection.

A Benchmark of Genetic Variant Calling Pipelines Using Metagenomic Short-Read Sequencing.

Microbes live in complex communities that are of major importance for environmental ecology, public health, and animal physiology and pathology. Short-read metagenomic shotgun sequencing is currently the state-of-the-art technique for exploring these communities. With the aid of metagenomics, our understanding of the microbiome is moving from composition toward functionality, even down to the genetic variant level. While the exploration of single-nucleotide variation in a genome is a standard procedure in genomics, and many sophisticated tools exist to perform this task, identification of genetic variation in metagenomes remains challenging. Major factors that hamper the widespread application of variant-calling analysis include low-depth sequencing of individual genomes (which is especially significant for the microorganisms present in low abundance), the existence of large genomic variation even within the same species, the absence of comprehensive reference genomes, and the noise introduced by next-generation sequencing errors. Some bioinformatics tools, such as metaSNV or InStrain, have been created to identify genetic variants in metagenomes, but the performance of these tools has not been systematically assessed or compared with the variant callers commonly used on single or pooled genomes. In this study, we benchmark seven bioinformatic tools for genetic variant calling in metagenomics data and assess their performance. To do so, we simulated metagenomic reads to mimic human microbial composition, sequencing errors, and genetic variability. We also simulated different conditions, including low and high depth of coverage and unique or multiple strains per species. Our analysis of the simulated data shows that probabilistic method-based tools such as HaplotypeCaller and Mutect2 from the GATK toolset show the best performance. By applying these tools to longitudinal gut microbiome data from the Human Microbiome Project, we show that the genetic similarity between longitudinal samples from the same individuals is significantly greater than the similarity between samples from different individuals. Our benchmark shows that probabilistic tools can be used to call metagenomes, and we recommend the use of GATK's tools as reliable variant callers for metagenomic samples.