key: cord-0947624-gbrx801u authors: Kurina, Irina; Popenko, Anna; Klimenko, Natalia; Koshechkin, Stanislav; Chuprikova, Liubov; Filipenko, Maxim; Tyakht, Alexander; Alexeev, Dmitry title: Development of qPCR platform with probes for quantifying prevalent and biomedically relevant human gut microbial taxa date: 2020-04-15 journal: Mol Cell Probes DOI: 10.1016/j.mcp.2020.101570 sha: b810e561131d85023f028d4c09da6209190f59ed doc_id: 947624 cord_uid: gbrx801u Abstract Nowadays the advent of innovative high-throughput sequencing allows obtaining high-quality microbiome profiling. However, PCR-based tests are still considered the “golden standard” for many clinical applications. Here, we designed a qPCR-based platform with fluorescent-labeled oligonucleotide probes for assessing human gut microbiome composition. The system allows conducting qualitative and semiquantitative analysis for 12 prokaryotic taxa that are prevalent in the human gut and associated with diseases, diet, age and other factors. The platform was validated by comparing microbiome profile data obtained with two different methods - the platform and high-throughput 16S rRNA sequencing - across 42 stool samples. The test can form the basis for precise and cost-efficient microbiome assay for large-scale surveys including clinical trials with interventions related to diet and disease risks. The gut microbiome is virtually an organ of the human body. The gut is colonized with 10 13 of microbes belonging to thousands of species. Alterations of microbiome composition are linked to lifestyle changes, short-and long-term diet [1] , probiotics intake [1, 2] , medications including antibiotics [3] and diseases including obesity, type 2 diabetes [4] , colorectal cancer [5] , Crohn's disease [6] and others. Based on the discovered associations, a number of microbial taxa have been proposed as potential biomarkers. For example, the abundance of Akkermansia muciniphila, a mucin-degrading bacterium is inversely associated with body fat mass and glucose intolerance [7] . Disrupted gut community structure is also one of the risk factors for cardiovascular diseases [8] . Many of the bacteria belonging to Clostridiales order are efficient producers of butyrate an other short-chain fatty acids known to possess anti-inflammatory, immunomodulatory and other important activities including gut epithelium integrity maintainence [9, 10] . The task of typing and enumerating gut microbial taxa is highly important. There are a plethora of methods for addressing it, including cultivation-based methods, PCR approaches and high-thoughput sequencing. Despite the recent progress with culturomics [11] , routine stool cultivation allows to enumerate a narrow range of gut microbes with low accuracy. The methods based on microbiome sequencing are becoming more popular -the "shotgun" metagenomics as well as amplicon-based surveys (using 16S rRNA gene for prokaryotes and ITS region for yeast, etc.). Unlike the culture-dependent approaches, the sequencing gives a total semiquantitative portrait of the microbiomecapturing the uncultivable diversity as well. An alternative method for taxonomic and functional profiling of microbial communities are hybridization microarrays that allow highthroughput detection of thousands of targets at a time including prokaryotes, viruses, yeast and protozoa [12] [13] [14] . Although microarrays capture low-abundant microbes and can be cost-efficient on a large scale, they only detect a fixed set of microbial targets and require a specialized microarray analysis instrument. Another category for microbiome profiling methods is based on PCR using oligonucleotide sequences specific to selected microbial clades. The real-time PCR (qPCR) is superior to the above-mentioned methods in terms of generally available equipment, cheaper sample preparation, flexibility in application and robustness for diagnostic applications. Human gut-specific PCR panels have been described before. One of them describes a qPCR test for detecting 4 taxa using SYBR Green: Bacteroidetes, Firmicutes, Bacteroides thetaiotaomicron and Enterococcus spp. [15] . In a more extensive manner, a real-time quantitative PCRbased analysis platform termed 'GUt Low-Density Array' (GULDA) with SYBR Green for 27 bacterial taxa was introduced [16] . Detailed cladefocused applications are available -for example, for detecting all known species of human gut-dwelling bifidobacteria [17] using genus-and species-specific primers and probes for qPCR and further proceeding with PCR product restriction for species identification. In this article we describe a qPCR-based platform for taxonomic profiling of gut microbiome targeting taxa that are prevalent in the gut and/or associated with important diseases, diet, populations, age and other factors. The platform provides a flexible tool for analyzing an individual microbiome composition (e.g., from a stool sample) to assess its general deviation from the reference population distributions, as well as quantify potential microbial biomarkers of disbyosis -universal as well as disease-specific. The target taxa were selected based on the principles of combining the taxa differentially abundant in subjects with major diseases (obesity, type 2 diabetes, inflammatory bowel diseases [IBD] and atherosclerosis), the taxa enriched in probiotic representatives and the "driver" taxa of the human gut microbiome. We collected a statistical summary of the targets' abundance in the 16S rRNA datasets of stool samples collected from the Russian general population during the OhMyGut crowdfunded project [1] and other studies (see Supplementary Table 2 ). The HITdb v1.00 database was used as a non-redundant curated set of 16S rRNA genes specific for human gut microbiota [18] . In order to identify primers and probes targeting specific groups of gut microbes, the gene sequences were analyzed to identify the regions that were both highly conservative within a given group and sufficiently variable to separate the group from other members of the gut microbiome. On the first step, nucleotide sequences belonging to a certain taxon were extracted from the database and aligned using the Aligner tool [19] . The obtained alignment was then used to define regions conservative within a given taxon. For this, we split the alignment into overlapping 20-mer blocks (set of columns) and only highly conservative blocks (showing 90% agreement in the consensus sequence with < 3 degenerate positions) were selected. On the next step, we checked the taxon-specificity of each conservative block. For this, a BLAST database was constructed using all the sequences from HITdb except those for the analyzed taxon. A conservative block was considered taxon-specific if, during a search against the BLAST database, it had aligned to < 1% of non-specific targets with 100% identity and coverage. In case when no such taxonspecific conservative blocks were identified, all the blocks were combined into pairs of blocks separated from each other by 50-300 bases (typical qPCR fragment size). A pair of blocks was considered taxonspecific if the blocks had not aligned to the same non-specific target simultaneously. Stool DNA was extracted as described previously [1] . For facilitating PCR amplification, additional purification from PCR inhibitors was performed using AmpliPrime DNK-Sorb-B (Federal State Institution Central Research Institute of Epidemiology of Rospotrebnadzor, Russia) according to the manufacturer's instructions. A 1200-1400 bp fragment of 16S rDNA was amplified as follows. The master mix in a reaction volume of 50 μL included: 1x Taq polymerase buffer (65 mM Tris-НCl [pH 8.9]); 16 mM (NH 4 ) 2 SO 4 ; 0.05% Tween 20; 3.5 mM MgCl 2 , 0.2 mM dNTPs, 50 ng genomic bacterial DNA or DNA isolated from human feces, 0.3 μM solutions of oligonucleotide primers, 1U Taq-polymerase (Biosan), 0.5U Pfu-polymerase (Biosan). Amplification was performed in Tercyc thermal cycler (DNK-Tehnologiya) under the following program: 15 min at 95°C for initial denaturation, 30 cycles of 10 s at 95°C for denaturation, 10 s at 60°C for annealing, 60 s at 72°C for elongation. The amplification products were hydrolyzed using the respective restriction endonucleases (SibEnzyme) and ligated using pBluscriptII SK(+) vector hydrolyzed with the same endonucleases within 3 h with 100U T4 DNA ligases (Biosan). Competent cells of E. coli strain XL1-Blue (Stratagene) were transformed with a ligase mixture. In the plasmid clones selected according to the results of specific PCR, the nucleotide sequence of the insert was determined by Sanger sequencing for confirmation. Plasmid DNA was extracted from 100 mL of overnight culture on LB medium using QIAGEN Plasmid Midi Kit (QIAGEN) according to the manufacturer's instructions. Concentrations of the obtained reference plasmid DNA were measured by spectrophotometry and fluorometry (Qubit BR kit, Invitrogen). Two μg DNA was hydrolyzed by restriction endonuclease SalGI for linearization. The obtained linear standards were diluted to a concentration of 10 7 to 10 0 copies of plasmid DNA per μL in a sterile buffer containing 10 mM TrisHCl pH 7.6 and poly-A RNA (5 ng per μL). DNA concentration in the obtained standards was measured using digital PCR on a QX100 Droplet Digital PCR System (Bio-Rad) platform according to the manufacturer's instructions. For this purpose, 20 μL PCR mixture containing the analyzed DNA (< 66 ng per 20 μL), 1X PCR mixture (Bio-Rad), 300 nM oligonucleotide primers and 100 nM TaqMan probe were prepared. To obtain microdroplets, 20 μl of the prepared PCR mixture and 70 μl of oil to generate drops were placed in open DG8 cartridges. 40 μL of the obtained microdroplets was transferred to a 96-well PCR plate, sealed with foil and placed in a thermal cycler. Amplification program: 96°C-10 min then 50 cycles were performed as follows: 96°C-30 s, 57°C-60 s with final warming for 10 min at 98°C. Then the microdroplets were read and analyzed using a Droplet Reader instrument; the obtained data were processed in the QuantaSoft program (Bio-Rad). Quantitative qPCR was carried out using a CFX96 thermal cycler (BioRad) in 96-well PCR plate (PCR-96M2-HS-C; Axygen). Sequences of primers and TaqMan probes are provided in Supplementary Tables 2 and 3 . Amplification for each sample was carried out in two replicates in a total volume of 20 μL containing: 65 mM Tris-HCl (pH 8.9), 16 mM (NH 4 ) 2 SO 4 , 3.0 mM MgCl 2 , 0.05% Tween 20, 0.2 mM dNTPs (BIOSAN LLC), 0.3 μM oligonucleotide primers with 0.1 μM TaqMan probe, 1 U hot start Taq polymerase (ICBFM SB RAS), 2 μL target DNA or plasmid standard. Amplification was carried out according to the following protocol: To quantify the DNA of the target microbial taxa in a sample, the method of calibration curves was used: successive 10-fold dilutions of plasmid standards were obtained and then amplified. The concentration range of calibration standards is shown in Supplementary Table 4 . Gene amplification levels obtained by PCR were recorded as Ct values that are inversely related to the initial DNA concentration. Quantity Mean was calculated automatically by the CFX96 (BioRad) software control program based on the set calibration curve values. After automatic data processing, the inaccuracies of automatic reading for some of the curves were corrected by manual alignment. The quantity values were exported to a spreadsheet for further analysis and comparison of data with the results of 16S rRNA gene sequencing. The library preparation for stool DNA samples and Illumina MiSeq high-throughput sequencing was performed as described previously [2] targeting V4 region of the gene with the following modification of 515F-806R primers: GTGBCAGCMGCCGCGGTAA and GACTACNVGG-GTMTCTAATCC. Data analysis has been done using Knomics-Biota system [20] . Reads were denoised using DADA2 algorithm [21] with the variable trimming length (from 251 to 253 bp). The taxonomic classification of the denoised reads has been done using QIIME2 naivebayes classifier [22] and GreenGenes database [18] with pre-trimmed and clustered with 97% identity sequences using TaxMan [23] and CD-HIT software (version 4.8.1) [24] . The relative abundance of each taxon was calculated by rarefying to 3000 reads per sample and normalizing to 100%. Proportions of genera and families were obtained by summing the levels of their member amplicon-sequencing variants (ASV). The qPCR and 16S rRNA values for target taxa abundance were compared using Spearman's correlation. Multiple comparison adjustment has been done using the Benjamini-Hochberg method. During PERMANOVA analysis, abundance values of PCR targets were used as predictors and Bray-Curtis dissimilarity matrix calculated from 16S rRNA composition -as output. The developed qPCR platform allows semiquantitative detection of major gut microbial taxa using qPCR with oligonucleotide probes. It provides an efficient way to characterize the major components of the human gut microbiome in relation to diseases, dietary interventions as well as for comparison with other factors. The platform includes specific oligonucleotide probes for detection and quantification of the microbes listed in Table 1 . The final list included 12 taxa shown in Table 1 . The number of taxa was selected to balance the scope of detection with the economical expenses due to the probe design and analysis cost. Following the in silico analysis described in the Materials and Methods, the consensus sequences were derived from taxon-specific conservative blocks (or pairs of blocks). Primers and probes were designed based on the consensus sequences using Vector NTI software [54] . The following requirements had to be fulfilled: melting temperature (T m ) -58-62°C for a primer and 65-75°C for a probe, GC content -40-65%, no hairpins with the ΔG < -3 kcal/mol and > 2 kcal/mol, no primer-dimers with the ΔG < −5 kcal/mol. The 12 pairs of primers and probes selected in this way are specific to the investigated taxa. The qPCR-based platform contains a pair of primers and an oligonucleotide probe for each taxon. The platform characteristics were calculated based on the amplification of consecutive 10-fold dilutions of plasmid standards. An example of amplification for one of the primer sets is shown in Fig. 1 . All primer sets worked in the efficiency range of 81-128.9%, the coefficient of determination was high (R 2 > 0.979) confirming high reproducibility of the results. Detailed analytical characteristics are listed in Supplementary Table 5 . Due to multitude of targeted taxa, the specificity of the qPCR platform was evaluated via comparison with the results of high-throughput 16S rRNA sequencing. For validation of the qPCR platform, we selected 42 stool samples from generally healthy Russian population that, according to previously performed 16S rRNA sequencing, were representative of global gut microbiome diversity and contained widely varying proportions of the target taxa of the qPCR platform. Each of the stool DNA samples were amplified using each of the primer pairs and probes (listed in Supplementary Tables 2 and 3 ). In each qPCR run, the analysis was performed for each taxon on all samples accompanied with negative control and a series of positive controls. The positive controls were plasmids in 4 dilutions each used to assess PCR efficacy (concentrations are shown in Supplementary Table 4 ). Two repeats were performed for each qPCR experiment to obtain more precise results. Standard curve was drawn by plotting the logarithmic input DNA concentration vs. mean cycle threshold (mCt). Overall, the differences between the replicates for stool samples were significantly lower than between samples (delta Ct 0.54 ± 0.29 vs. 5.31 ± 4.63, p = 0; see Supplementary Fig. S1 ). In case if the discrepancy between the two replicates is more than 1 cycle it is recommended that the experiment should be repeated and the outlying replicate be removed. The same stool DNA samples were also subject to 16S rRNA library preparation (V4 region) and sequencing. The taxonomic classification of each read was used to obtain the total composition of each sample and the relative abundance of each identified taxon. Overall, when the abundance values across all taxa and samples were pooled, the correlation between the qPCR and 16S rRNA sequencing was significant (Spearman's r = 0.66, p = 0). For 11 of the 12 target taxa, their levels were significantly correlated between the two methods (FDR-adjusted p < 0.05) (see Fig. 2 and Supplementary Tables 6 and 7; per-sample comparison is shown in Supplementary Fig. S2 ). The remaining taxon [Ruminococcus] produced insignificant correlation (r = 0.03, FDR-adjusted p = 0.8414). Fig. 3 shows the percentage of variance in community structure according to NGS explained by the level of each PCR target; totally all targets explain 36.5% of it. Considerable conformance of our qPCR platform with the results of high-throughput sequencing suggests that the panel is a prototype applicable for profiling major components of the human gut microbiome. The analyzed sample types can range from stool samples -as the noninvasive method of choice in such studies -to swabs and biopsies as well as contents of human intestinal tract simulators. Large-scale datasets obtained via the standard described protocol and using the same DNA extraction protocol can be subsequently subject to statistical analysis like the ones used for NGS microbiome data to conduct exploratory I. Kurina, et al. analysis and reveal associations with factors like diseases, diet and lifestyle [20] . The platform was designed in a way to reach equilibrium between the coverage of the human gut microbiome diversity and cost per analysis. Subsequently, as the gut microbiome is highly complex, many prevalent taxa are not currently included among the targets. According to our comparison with NGS results, only a third of the total composition variance is explained by the platform. Our results show that the other families and genera that are abundant and variable across the test samples according to 16S rRNA sequencing but not targeted by the PCR test include Lachnospiraceae spp. (9.2 ± 4.3%), Faecalibacterium (5.9 ± 5.0%) and Ruminococcus (3.4 ± 3.0%). These along with their members are potential candidates to be included in the prospective extended panel -along with various members of the Firmicutes phylum including Dorea, Roseburia and others. However, the taxonomy within the phylum especially its Clostridiales order is rather ambiguous hence further elaboration of the primers and probes will be required. One of the methods for microbiome profiling that can be viewed as an alternative to the proposed panel is Nanopore sequencing. Its advantages includes the possibility to obtain the total microbial composition in an untargeted way. However, although the analysis time (~12 h) takes shorter than for 16S rRNA on Illumina platform, compared to qPCR it takes longer (qPCR: 2-3 h), is less cost-efficient and scalable and requires more expensive equipment [55] . For some of the taxa, we observed a reduced correlation of qPCR levels with the 16S rRNA sequencing results. At least for some of the taxa they could be due to a lower dynamic range provided by NGS (at given coverage) -so that the values below 1 read could be detected only by qPCR. A low correlation for [Ruminococcus] could be due to abovedescribed complexity of Clostridiales phylogeny and need to improve the alignment of taxonomic annotations with the 16S rRNA database. While our method is cultivation-independent and allows to capture a variety of taxa undetected during a classic stool test, its taxonomic detection accuracy is limited by the level of genus. A single genus can include species and strains with a highly variable gene content resulting in dramatically different phenotypes -from probiotic to human pathogen like in the case of Escherichia coli [6] -or correlating to opposing dietary patterns as shown for Prevotella copri [56] . Despite this limitation, a general-level profile of the gut microbiome composition provided by the proposed platform can be useful for assessing the general degree of dysbiosis. It can help in more advanced tasks like predicting response to drug therapy (as suggested by the discovered link between cancer immunotherapy response and Bacteroides levels [57] ) or diet (for example, high-fiber diet [1] ), measuring contribution of gut microbes to microbiome of built environments [58] as well as rapidly assessing the gut microbiome dysbiosis in patients with viral infections like COVID- Table 1 List of the 12 target taxa in the platform. Associations References Inversely associated with metabolic signatures of obesity including insulin resistance and fasting glucose level. Suggested marker for calorie restriction diet efficiency prediction. [7, [25] [26] [27] [28] . One of the dominant genera in adult gut microbiome; linked to "Western diet". Can change in abundance after high-fiber dietary interventions. [1, 29, 30] Bifidobacterium Dominates in the gut community of healthy infants. Inversely associated with obesity and Crohn's disease. Abundant in certain Eurasian populations. Contains many probiotic strains. Can increase in abundance after intervention with probiotics. [31] [2,32-34] One of the dominant genera in adult gut microbiome. Linked to type 2 diabetes. [4, 33] Christensenellaceae Inversely associated with obesity and IBD and directly -to Parkinson's disease; members of the family can promote reduction of excessive weight. [32, 35, 36] Enterobacteriaceae Positively associated with metabolic disorders and IBD. [6, 32, 37, 38] Linked to inflammation, can carry a wide range of virulence factors and antibiotic resistance genes. Its intestinal domination increases the risk of bacteremia in immunocompromised patients. Prevalent in infant microbiome. [39, 40] ; [41] . Inversely associated with type 2 diabetes, ulcerative colitis. Contains many probiotic strains. [25, 42] Linked to eating behavior, anorexia nervosa, body mass index, obesity, IBD. [43] [44] [45] One of the dominant genera in adult gut microbiome. Increased in rural populations of the world. Linked to diet high in fiber and low levels in animal protein and fat. Maternal carriage during pregnancy is reversely linked to food allergy in the offspring. [ [46] [47] [48] [Ruminococcus] (from Lachnospiraceae family) High abundance is linked to IBD. [49, 50] Streptococcus High abundance is associated with liver cirrhosis, as well as to the increased risk of bacteremia in immunocompromised patients. [51, 52] . One of the problems of the PCR approach is related to normalization. To address this issue, the set could be extended by adding a universal pair of primers (and a probe) targeting total prokaryotic content that can be used for normalization purposes (for example, by dividing signals from other taxa by it). Although this would reflect the microbial concentration in the analyzed DNA sample, this concentration could not directly correspond to the concentration in the subject's stool, as it can change considerably during the extraction. One of the suggested ways to approach absolute normalization is to collect an exactly equal volume of stool each time. Still, our platform allows us to assess the pairwise ratios between the targeted taxa. Consideration of ratios in the context of microbiome data is especially relevant as for many techniques -including metagenomics -they are often inherently compositional implying use of statistical approaches based on ratios, for instance, the balances [60] . The ratios and their derivatives can further be used as microbial community features to explore the associations with clinical status, diet and other factors. Finally, the proposed platform can be used for prioritization purposes: basing on its results for a specific subject, one can opt for performing additional in-depth analysis using other qPCR platforms targeting microbial genes (like virulence factors and antibiotic resistance determinants) or pathogenic taxa. The developed qPCR platform is a perspective technique for rapid evaluation of human gut microbiome composition. Positive correlation with the data obtained using a different method demonstrates that the qPCR method can be applied for larger cohorts during the clinical and observational studies related to diet and diseases. Further larger clinical trials are needed to identify the method specificity and accuracy when This work was financially supported by the Skolkovo Foundation for Development of the Center for Elaboration and Commercialization of New Technologies (#G94/16 to Knomics LLC). Microbiome responses to an uncontrolled short-term diet intervention in the frame of the citizen science project Human gut microbiome response induced by fermented dairy product intake in healthy volunteers Shifts in the human gut microbiota structure caused by quadruple eradication therapy Gut microbiota and diet in patients with different glucose tolerance Meta-analysis of fecal metagenomes reveals global microbial signatures that are specific for colorectal cancer Genetic diversity of Escherichia coli in gut microbiota of patients with Crohn's disease discovered using metagenomic and genomic analyses Akkermansia muciniphila and improved metabolic health during a dietary intervention in obesity: relationship with gut microbiome richness and ecology Symptomatic atherosclerosis is associated with an altered gut metagenome A decrease of the butyrate-producing species Roseburia hominis and Faecalibacterium prausnitzii defines dysbiosis in patients with ulcerative colitis Causal relationships among the gut microbiome, short-chain fatty acids and metabolic diseases Culturing the human microbiota and culturomics Axiom microbiome Array, the next generation microarray for high-throughput pathogen and microbiome analysis Development and validation of a microarray for the investigation of the CAZymes encoded by the human gut microbiome Development and application of the human intestinal tract chip, a phylogenetic microarray: analysis of universally conserved phylotypes in the abundant microbiota of young and elderly adults Tine Rask Licht, Freezing fecal samples prior to DNA extraction affects the Firmicutes to Bacteroidetes ratio determined by downstream quantitative PCR analysis, FEMS (Fed Introducing GUt low-density Array (GULDA) -a validated approach for qPCR-based intestinal microbial community analysis, FEMS (Fed A PCR-based method for identification of bifidobacteria from the human alimentary tract at the species level, FEMS (Fed Improved taxonomic assignment of human intestinal 16S rRNA sequences by a dedicated reference database Ribosomal database project: data and tools for high throughput rRNA analysis Knomics-biota -a system for exploratory analysis of human gut microbiota data DADA2: high-resolution sample inference from Illumina amplicon data Author correction: reproducible, interactive, scalable and extensible microbiome data science using QIIME 2 TaxMan: a server to trim rRNA reference databases and inspect taxonomic coverage CD-HIT: accelerated for clustering the next-generation sequencing data Disentangling type 2 diabetes and metformin treatment signatures in the human gut microbiota Cross-talk between Akkermansia muciniphila and intestinal epithelium controls diet-induced obesity Akkermansia muciniphila: key player in metabolic and gastrointestinal disorders Molecular and Cellular Probes xxx (xxxx) xxxx Specific gut microbiota features and metabolic markers in postmenopausal women with obesity Diet rapidly and reproducibly alters the human gut microbiome Linking long-term dietary patterns with gut microbial enterotypes Dynamics and stabilization of the human gut microbiome during the first year of life Interplay of host genetics and gut microbiota underlying the onset and clinical presentation of inflammatory bowel disease Human gut microbiota community structures in urban and rural populations in Russia Metabolic activities and probiotic potential of bifidobacteria Human genetics shape the gut microbiome Analysis of gut microbiota in patients with Parkinson's disease An opportunistic pathogen isolated from the gut of an obese human causes obesity in germfree mice Host-mediated inflammation disrupts the intestinal microbiota and promotes the overgrowth of enterobacteriaceae Comparative genomics of enterococci: variation in Enterococcus faecalis, clade structure in E. Faecium, and defining characteristics of E. Gallinarum and E. Casseliflavus Diversification and evolution of vancomycin-resistant Enterococcus faecium during intestinal domination Intestinal microflora in early infancy: composition and development The use of Lactobacillus casei and Lactobacillus paracasei in clinical trials for the improvement of human health Correlation between body mass index and gut concentrations of Lactobacillus reuteri, Bifidobacterium animalis, methanobrevibacter smithii and Escherichia coli A systematic review of studies on the faecal microbiota in anorexia nervosa: future Research may need to include microbiota from the small intestine Gut colonization with methanobrevibacter smithii is associated with childhood weight development Impact of diet in shaping gut microbiota revealed by a comparative study in children from Europe and rural africa Gut microbiome of the hadza hunter-gatherers Maternal carriage of Prevotella during pregnancy associates with protection against food allergy in the offspring Discovery of intramolecular trans-sialidases in human gut microbiota suggests novel mechanisms of mucosal adaptation Mucolytic bacteria with increased prevalence in IBD mucosa augment in vitro utilization of mucin by other bacteria Alterations of the human gut microbiome in liver cirrhosis Intestinal domination and the risk of bacteremia in patients undergoing allogeneic hematopoietic stem cell transplantation Links of gut microbiota composition with alcohol dependence syndrome and alcoholic liver disease Vector NTI, a balanced all-in-one sequence analysis suite Intawat Nookaew, An assessment of oxford Nanopore sequencing for human gut metagenome profiling: a pilot study of head and neck cancer patients Distinct genetic and functional traits of human intestinal Prevotella copri strains are associated with different habitual diets Gut microbiome modulates response to anti-PD-1 immunotherapy in melanoma patients Co-occurrence patterns of bacteria within microbiome of moscow subway Management of corona virus disease-19 (COVID-19): the Zhejiang experience Balances: a New perspective for microbiome analysis, mSystems Molecular and Cellular Probes xxx (xxxx) xxxx The authors would like to thank Maria Kardakova for performing language editing. Supplementary data to this article can be found online at https:// doi.org/10.1016/j.mcp.2020.101570. The authors declare that they have no competing interests.