Strongyloides stercoralis infection induces gut dysbiosis in chronic kidney disease patients
============================================================================================

* Nguyen Thi Hai
* Nuttanan Hongsrichan
* Kitti Intuyod
* Porntip Pinlaor
* Manachai Yingklang
* Apisit Chaidee
* Thatsanapong Pongking
* Sirirat Anutrakulchai
* Ubon Cha’on
* Somchai Pinlaor

## Abstract

**Background** *Strongyloides stercoralis* infection typically causes severe symptoms in immunocompromised patients. However, the progression of infection-driven chronic kidney disease (CKD) is not understood fully. Recent studies have shown that gut dysbiosis plays an important role in the progression of CKD. Hence, this study aims to investigate the effect of *S. stercoralis* infection on the gut microbiome in CKD patients.

**Methodology/Principal Findings** Among 838 volunteers from Khon Kaen Province, northeastern Thailand, 40 subjects with CKD were enrolled and divided into two groups (*S. stercoralis*-infected and -uninfected) matched for age, sex and biochemical parameters. Next-generation technology was used to amplify and sequence the V3-V4 region of the 16S rRNA gene to provide a profile of the gut microbiota. Results revealed that members of the *S. stercoralis*-infected group had lower gut microbial diversity than was seen in the uninfected group. Interestingly, there was significantly greater representation of some pathogenic bacteria in the *S. stercoralis*-infected CKD group, including *Escherichia-Shigella* (*P* = 0.013), *Rothia* (*P* = 0.013) and *Aggregatibacter* (*P* = 0.03). There was also a trend towards increased *Actinomyces, Streptococcus* and *Haemophilus* (*P* > 0.05) in this group. On the other hand, the *S. stercoralis*-infected CKD group had significantly lower representation of SCFA-producing bacteria such as *Anaerostipes* (*P* = 0.01), *Coprococcus*_1 (0.043) and a non-significant decrease of *Akkermansia, Eubacterium rectale* and *Eubacterium hallii* (*P* > 0.05) relative to the uninfected group. Interesting, the genera *Escherichia-Shigella* and *Anaerostipes* exhibited opposing trends, which were significantly related to sex, age, infection status and CKD stages. The genus *Escherichia-Shigella* was significantly more abundant in CKD patients over the age of 65 years and infected with *S. stercoralis*. A correlation analysis showed inverse moderate correlation between the abundance of the genus of *Escherichia-Shigella* and the level of estimated glomerular filtration rate (eGFR).

**Conclusions/Significance** Conclusion, the results suggest that *S. stercoralis* infection induced gut dysbiosis in the CKD patients, which might be involved in CKD progression.

**Author summary** Human strongyloidiasis is caused by a soil-transmitted helminth, *Strongyloides stercoralis*, which typically causes severe symptoms in immunocompromised. However the relationship between *S. stercoralis* and chronic kidney disease patients (CKD) progression was not known. This is the first study to investigate the gut microbiota of CKD patients with and without *S. stercoralis* using high-throughput sequencing of the V3–V4 region of the 16S rRNA gene. Infection with *S. stercoralis* was associated with reduced gut microbiota diversity than in the uninfected group. In addition, infection with this nematode led to reduced abundance of SCFA-producing bacteria and enrichment of pathogenic bacteria. In particular, there were significant differences in abundance of the beneficial genus *Anaerostipes* (a decrease) and the pathogenic taxon *Escherichia-Shigella* (an increase) in CKD patients infected with *S. stercoralis* relative to controls. In the infected group, the representation of genus *Escherichia-Shigella* was significant higher in patients over the age of 65 years. There was a significant inverse moderate correlation of *Escherichia-Shigella* with the estimated glomerular filtration rate (eGFR).

Key words
*   *Strongyloides stercoralis*
*   chronic kidney disease
*   gut microbiota
*   16S rRNA gene
*   next-generation technology

## Introduction

An imbalance within the microbiota in the gastrointestinal tract, termed gut dysbiosis, contributes to the development and progression of many diseases including chronic kidney disease (CKD) [1]. Many studies have shown a significant difference in the abundance of bacterial populations in the gastrointestinal tract (GI) between CKD and control individuals. Substantially lower proportions of *Bifidobacterium, Lactobacillaceae, Bacteroidaceae* and *Prevotellaceae* were seen in CKD patients, including those undergoing hemodialysis, while the proportions of *Enterobacteriaceae*, especially *Enterobacter, Klebsiella* and *Escherichia*, were notably higher [2-5]. The production of uremic toxins (indoxyl sulphate (IS), trimethylamine-N-oxide (TMAO)), which results from nutrient processing by gut microbiota, and the reduction of fiber-derived short-chain fatty acids, are linked with CKD progression [6-8]. Recent studies have found various factors involved in microbial dysbiosis and CKD, such as the use of antibiotics [9], decreased consumption of dietary fiber [10], and oral iron intake [11]. However, many etiological factors associated with CKD remain obscure [12], particularly those due to infection with intestinal parasites.

The ability of GI parasitic infection to changes the gut microbiota and host-microbiota interactions has been clearly identified. Infection with intestinal parasites either induced gut dysbiosis or provided protection against dysbiosis and inflammatory disease [13]. *Strongyloides stercoralis* is one of the most medically important parasites in northeastern Thailand, where the prevalence of CKD is also high [14]. Typically, *S. stercoralis* infection causes only mild GI symptoms. However, when immunity is suppressed by, for example, CKD, diabetes mellitus or HIV infection, the parasite can rapidly multiply and eventually leads to hyperinfection and disseminated strongyloidiasis, which is a life-threatening condition [15-18].

Recent studies have successfully demonstrated that *S. stercoralis* induces an increase in bacterial diversity and changes faecal microbiota [19, 20]. By using metagenomic analysis, microbial alpha diversity was found to increase, and beta diversity decrease, in the faecal microbial profiles of *S. stercoralis*-infected individuals compared to uninfected. Faecal metabolite analysis detected marked increases in the abundance of selected amino acids and decrease in short-chain fatty acids in *S. stercoralis* infection, relative to uninfected controls [19]. Taken together, we therefore hypothesize that *S. stercoralis* infection changes gut microbiome, contributing to progression of CKD.

To test this hypothesis, metagenomic analysis was done in patients with CKD to investigate the changes in the gut microbiota that can be attributed to *S. stercoralis* infection. The result from this study might be useful for identifying strategies for prevention of chronic kidney disease development and progression.

## Methods

### Study population, area and ethics statement

The study was conducted between January 2017 and May 2018 at Donchang sub-district, Khon Kaen Province, northeastern Thailand as a part of the Chronic Kidney Disease Northeastern Thailand (CKDNET) project. Included in the study were individuals (>35 years of age) of chronic kidney disease. Their diagnosis, done by a nephrologist [14], included clinically proven impaired kidney structure or renal function, as detected using ultrasonography, and a finding of reduced eGFR. The staging (stages 1 to 5) based on the estimated glomerular filtration rate (eGFR) [21, 22] was estimated for each individual. In patients with eGFR > 60 ml/min/1.73 m2, kidney damage was confirmed based on urine albumin-to-creatinine ratio (UACR), hematuria and abnormal renal ultrasound. Stool examination was performed on CKD patients using the modified formalin ethyl acetate concentration technique (FECT) and modified agar plate culture (mAPC) as previously reported [23].

Exclusion criteria were as follows: use of antibiotics or probiotics, diabetes, autoimmune disease, urinary tract infection, renal infection with other intestinal parasites. Twenty CKD patients with *S. stercoralis* infection (Ss+) [24] met these criteria and were included. The control group consisted of 20 CKD patients free of *S. stercoralis* infection (Ss−) who otherwise matched the characteristics (including sex, age and biochemical factors) of the Ss+ group (Table 1). These datasets were obtained from the medical records of CKDNET and from a recent study [23]. Absence of *S. stercoralis* infection in members of the control group was also confirmed using PCR tests. The human ethical review committee of Khon Kaen University (HE631200) approved the protocol of study.

View this table:
[Table 1.](http://medrxiv.org/content/early/2022/03/09/2022.03.07.22271995/T1)

Table 1. Characteristics of chronic kidney disease patients.
Data are presented as mean ± standard deviation of the mean. Independent t-tests and Mann-Whitney U tests (a) were used to calculate *P* values.

View this table:
[Table 2.](http://medrxiv.org/content/early/2022/03/09/2022.03.07.22271995/T2)

Table 2. Alpha diversity of the gut microbiota in Ss− and Ss+ groups, calculated according to several indices.

### Sample collection, DNA extraction

Faecal and serum samples (n=40) were collected from the CKD Northeast Thailand project and kept at −20 °C until analysis. DNA was extracted from faecal samples using the QIAamp Kit (Qiagen, Germany). A Nanodrop 2000 spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA) was used to measure DNA concentration and 1.5% agarose gel electrophoresis was used to check the DNA quality.

### 16S rRNA gene sequencing and analysis

To confirm that V3-V4 regions (about 450-500 bp) of the prokaryotic 16S rRNA gene could be amplified, PCR was used to amplify DNA from each faecal sample. Specific V3-V4 primers were used, V3-forward: 5’-CCTACGGGNGGCWGCAG and V4-reverse: 5’-TACNVGGGTATCTAATCC [25]. Each reaction (20.1 μL) contained 2 μL of 10x buffer MgCl2, 0.4 μL of 50mM MgCl2, 0.6 μL of 10mM dNTP, 1 μL of 5 μM of each primer, 0.1μL Platinum *Taq* DNA polymerase and distilled water. The amplification profile was initial denaturation at 94 °C for 5 min, at 94 °C for 40 sec, then 35 cycles of 52.8 °C for 30 sec, 72 °C for 2 min, followed by a final extension at 72 °C for 10 min. PCR product was electrophoresed in a 1.5% agarose gel to confirm the presence of a band of the expected size.

A sequencing library was generated for each sample using NEBNext® UltraTM DNA Library Prep Kit for Illumina (Thermo Scientific) following the manufacturer’s recommendations. The library quality was assessed on the Qubit@ 2.0 Fluorometer (Thermo Scientific) and Agilent Bioanalyzer 2100 system. Finally, the library was sequenced on an Illumina platform and 250 bp paired-end reads were generated. Processing and quality control of these reads used the following steps: 1) Data split (based on their unique barcode) and truncation by cutting off the barcode and primer sequences; 2) Sequence assembly (Paired-end reads were merged using FLASH [26] and generated the raw tags; 3) Data filtration (quality filtering on the raw tags was performed to obtain high-quality clean tags [27] according to the QIIME(V1.7.0, [http://qiime.org/index.html](http://qiime.org/index.html)) quality controlled process [28]; 4) Chimera removal (the tags were compared with a reference database [29] using the UCHIME algorithm [30] to detect chimera sequences, which were then removed [31].

Sequences sharing ≥97% similarity were assigned to the same operational taxonomic unit (OTU) by using Uparse v7.0.1001) [32] and species annotated by reference to the GreenGene Database[33] based on RDP 3 classifier [34] algorithm to annotate taxonomic information. Phylogenetic relationship construction was conducted using the MUSCLE software (Version 3.8.31) [35]. Information on abundance of each OTU was normalized relative to the sample with the fewest sequences. Subsequent analyses of alpha diversity were performed based on this normalized data. Alpha diversity indicates species diversity for a sample using six indices (observed-species, Chao1, Shannon, Simpson, ACE and Good’s coverage). All these indices were calculated using QIIME (Version 1.7.0) and displayed with the help of R software (Version 2.15.3). To show beta diversity, weighted and unweighted UniFrac metrics were used to evaluate differences of samples in species complexity by using QIIME software (Version 1.7.0). Principal Coordinate Analysis (PCoA) was performed to visualize complex and multidimensional data. PCoA analysis was done using the WGCNA package, stat packages and ggplot2 package in R software (Version 2.15.3). Metastats was used to detect taxa with significant intra-group variation. The potential biomarkers were detected by using LEfSe (linear discriminant analysis (LDA) Effect Size) [36]. Raw data are available from Mendeley Data (DOI: 10.17632/hvbvrtc34x.1)

### Polymerase Chain Reaction (PCR)

Conventional PCR was used to confirm the absence of *S. stercoralis* infection in members of the control group (n=20). Primers were designed specifically to amplify a 125 bp fragment from a *S. stercoralis* dispersed repetitive sequence, GenBank: [AY08262](http://medrxiv.org/lookup/external-ref?link_type=GEN&access_num=AY08262&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F09%2F2022.03.07.22271995.atom) [37]. The forward primer was SSC-F 5′ CTCAGCTCCAGTAAAGCAACAG 3′ and reverse primer was SSC-R 5′ AGCTGAATCTGGAGAGTG AAGA 3′. PCR amplification was performed in a 12.5 μL volume with Dream Taq PCR Master-mix (Thermo Fisher Scientific, Vilnius, Lithuania), 1 μL of each primer, 5 μL of 9–155 ng/μL DNA, and PCR-grade water. The amplification profile was initial denaturation at 95 °C for 10 min, followed by 35 cycles of 95 °C for 1 min, 60 °C for 1 min 30 s and 72 °C for 1 min; then a final extension at 72 °C for 10 min. To confirm amplification and amplicon size, the PCR products were resolved on a 2% agarose gel stained with ethidium bromide.

### Statistical analysis

Statistical analyses including ANOVA, independent t-test, Welch’s t-test, Kruskal-Wallis test, Mann-Whitney U tests and Pearson’s correlation coefficient were conducted using IBM SPSS Statistics version 20 (IBM, Armonk, New York). Statistically significant features were further examined with post-hoc tests (Tukey-Kramer) to determine which groups of profiles differed from each other. Non-normal distribution data of two genera (*Anaerostipes* and *Escherichia-Shigella*) in association with CKD were log-transformed into normal distributions. Statistical significance and 95% confidence intervals (95%CI) were calculated and considered as *P* < 0.05.

## Results

### Study population characteristics

Demographic, socioeconomic and clinical characteristics of CKD patients with and without *S. stercoralis* infection were matched (Table 1). No significant differences in these characteristics were found between the Ss+ and Ss− groups.

### Characterization of bacterial diversity and community structure

In total, 1551 OTUs were identified based on the 97%-similarity rule, with an average of 477 OTUs per sample. Sequences were classified into 16 bacterial phyla, 26 classes, 45 orders, 72 families and 258 genera and 189 species (including unidentified species). The species accumulation curves showed a saturation phase (Fig. 1). This indicates that the sample size was sufficient to capture the overall microbiota structure.

![Fig 1.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/03/09/2022.03.07.22271995/F1.medium.gif)

[Fig 1.](http://medrxiv.org/content/early/2022/03/09/2022.03.07.22271995/F1)

Fig 1. Species accumulation curve.
X-axis: Number of samples, Y-axis: number of OTUs. Following an initial sharp rise in the number of OTUs as number of samples increases, there is a levelling of the plot. The narrow spread of the boxplots as the total number of samples is approached indicates that the number of samples was adequate to capture most of the microbial diversity present.

The unweighted UniFrac distances, reflecting beta diversity, showed significant differences between the Ss+ and Ss− groups (*P* = 0.00019). In terms of alpha diversity overall, there were no significant differences in estimated OTU richness, Chao1, the ACE metric, the Shannon diversity index and Good’s coverage, (*P* > 0.05) (Table 3). In contrast, the alpha diversity in males in the Ss− group was significantly higher than in males in the Ss+ group (Shannon diversity index, *P* = 0.015; Simpson diversity index, *P* = 0.058) (Fig 2B, C).

View this table:
[Table 3.](http://medrxiv.org/content/early/2022/03/09/2022.03.07.22271995/T3)

Table 3. Taxa in the gut microbiome differing significantly between CKD patients with and without *S. stercoralis*.

![Fig 2.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/03/09/2022.03.07.22271995/F2.medium.gif)

[Fig 2.](http://medrxiv.org/content/early/2022/03/09/2022.03.07.22271995/F2)

Fig 2. Comparison of alpha diversity indexes and beta diversity in CKD patients with and without *S. stercoralis* infection.
(A) Shannon index (B) Shannon index in males (C) Simpson index in males. (D) Boxplot based on unweighted UniFrac distance. (E) Principal coordinate analysis (PCoA).

The principal coordinate analysis (PCoA) was used to illustrate the beta diversity based on the unweighted UniFrac distances (Fig 2D). PCoA analysis revealed that the gut microbiota of Ss+ subjects showed deviation from the Ss− group (Fig 2E).

LDA score showed the significantly in abundance between the two groups. The candidate biomarker for the Ss− category was order Bradymonadales and for CKD with *S. stercoralis* infection (Ss+) was species *E. coli*, genera *Escherichia-Shigella, Dialister*, family Veillonellaceae and order Selenomonadales (Fig 3).

![Fig 3.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/03/09/2022.03.07.22271995/F3.medium.gif)

[Fig 3.](http://medrxiv.org/content/early/2022/03/09/2022.03.07.22271995/F3)

Fig 3. Linear discriminant analysis (LDA).
The histogram of the LDA scores presents taxa (potential biomarkers) whose abundance differed significantly among groups (Ss+ vs. Ss−) Oder Bradymonadales in Ss**−** (green colour). Species *E. coli* belong to genus *Escherichia-Shigella*; genus *Dialister* belong to order Selenomonadales, class Negativicutes and family Veillonellaceae in Ss+ (red colour).

### Differences in bacterial abundance between the Ss+ and Ss− groups

Proportions of sequence reads were compared between groups at the phylum and genus levels using Metastats. At the phylum level, there were no significant differences (Fig 4). For example, relative abundances of some principal taxa were: Firmicutes (Ss− 64.14% vs. Ss+ 59.33%; *P* = 0.39), Proteobacteria (Ss− 15.01% vs. Ss+ 18.83%; *P* = 0.50), Bacteroidetes (Ss− 12.61% vs. Ss+ 12.68%; *P* = 0.98), Actinobacteria (Ss− 3.46% vs. Ss+ 4.34%; *P* = 0.67) and Fusobacteria (Ss− 1.88% vs. Ss+ 3.81%; *P* = 0.52). However, at the genus level, 42 taxa were differentially relative abundance (Table 4).

![Fig 4.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/03/09/2022.03.07.22271995/F4.medium.gif)

[Fig 4.](http://medrxiv.org/content/early/2022/03/09/2022.03.07.22271995/F4)

Fig 4. The gut microbiota composition.
(A) and (B), Control group compared with *S. stercoralis-*infection group (Ss+) at the phylum and genus levels, respectively.

To study the similarity among different samples, clustering analysis was applied. The unweighted pair group method with arithmetic mean (UPGMA), a type of hierarchical clustering method widely used in ecology, Ss− vs. Ss+ samples tended to cluster separately (Fig 5).

![Fig 5.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/03/09/2022.03.07.22271995/F5.medium.gif)

[Fig 5.](http://medrxiv.org/content/early/2022/03/09/2022.03.07.22271995/F5)

Fig 5. Unweighted pair group method with arithmetic mean (UPGMA).
UPGMA cluster tree based on unweighted UniFrac distances between CKD patients with or without *S. stercoralis* infection. The red branches present *S. stercoralis* infection (Ss+) subject and the dark blue branches indicate un-infection (Ss−) ones.

### The trends of some bacteria in *S. stercoralis* infection

Fig 6 shows the comparisons of abundance of some bacteria between Ss− and Ss+ groups at the genus level, with pathogenic bacteria groups were increased in Ss+, such as: *Escherichia-Shigella* (3.36% vs. 13.33%; *P* < 0.05), *Streptococcus* (0.97% vs. 2.18%; *P >* 0.05), *Haemophilus* (0.46% vs 0.71%; *P >* 0.05), *Rothia* (0.024% vs. 0.11%; *P* < 0.05), *Actinomyces* (0.038% vs. 0.067%; *P >* 0.05), *Aggregatibacter* (0.0013% vs. 0.025%; *P* < 0.05) and reducing of some SCFA-producing bacteria were observed, such as: *Eubacterium rectale_group* (4.51% vs. 3.78%; *P >* 0.05), *Eubacterium hallii_group* (1.24% vs. 0.94%; *P >* 0.05), *Anaerostipes* (0.54% vs. 0.074%; *P* < 0.05), *Coprococcus*_1 (0.11% vs. 0.057%; *P* < 0.05) and *Akkermansia* (0.081% vs. 0.043%; *P >* 0.05).

![Fig 6.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/03/09/2022.03.07.22271995/F6.medium.gif)

[Fig 6.](http://medrxiv.org/content/early/2022/03/09/2022.03.07.22271995/F6)

Fig 6. Comparisons of abundance (numbers of sequence reads) of some bacteria between Ss− and Ss+ group.
Pathogenic bacteria: *Escherichia-Shigella, Streptococcus, Haemophilus, Rothia, Actinomyces, Aggregatibacter*. SCFA-producing bacteria: *Eubacterium rectale\_group, Eubacterium hallii\_group, Anaerostipes, Coprococcus_1, Akkermansia*.

### *Anaerostipes* and *Escherichia-Shigella* exhibit opposing trends in abundance and correlate with sex, age and CKD stage

Fig 7 shows the opposing trends in abundance in the genera *Anaerostipes* and *Escherichia-Shigella*. The proportion of reads of *Anaerostipes* was lower in those aged over 65 years (Fig. 7A), in females (7B), and with increasingly advanced CKD stage (7C) and in those infected with *S. stercoralis* (Fig. 7A, B, C). The contrary was observed in the case of *Escherichia-Shigella*: higher proportions of sequence reads of this genus were seen in elderly (>65 years) individuals (Fig. 7 D), in females (Fig. 7 E) and in those infected with *S. stercoralis* (Fig. 7D, E, F). More reads of *Escherichia-Shigella* were seen with increasingly advanced CKD stage (Fig. 7F).

![Fig 7.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/03/09/2022.03.07.22271995/F7.medium.gif)

[Fig 7.](http://medrxiv.org/content/early/2022/03/09/2022.03.07.22271995/F7)

Fig 7. Opposing trends in abundance of two genera, *Anaerostipes* and *Escherichia-Shigella*.
(A), (D) Different trend related to age (B), (E) sex and (C), (F) CKD stages. * *P* < 0.05, ***P*<0.01, \***|*P*<0.001. Analysis of the difference among groups of sex, age and CKD stages based on one-way ANOVA test.

## Discussion

In this study, we first characterized the gut microbiota of CKD with and without *S. stercoralis* using high-throughput sequencing of the V3–V4 region of the 16S rRNA gene. The results showed that *S. stercoralis* infection altered gut microbiota composition in CKD patients, leading to lower microbial diversity. This study also suggested that microbial candidate biomarkers for CKD concurrent with *S. stercoralis* infection include *Escherichia coli* (genus *Escherichia-Shigella*, phylum Proteobacteria) and the genus *Dialister* (family Veilslonellaceae, order Selenomonadales, class Negativicutes, phylum Firmicutes).

Our results revealed that gut microbiota was altered by *S. stercoralis* infection in CKD patients. The alpha-diversity indices (Chao1, the ACE metric, the Shannon diversity index, Good’s coverage) did not significantly differ between the two groups. However, the Shannon diversity index in males (n=12) in Ss+ group was significantly lower than in the Ss− group. In addition, the beta diversity, based on the unweighted UniFrac distances, was significantly lower in the Ss+ group, suggesting a decrease in ecological diversity in CKD concurrent with *S. stercoralis* infection. These slight differences spanned all taxonomic levels of the microbiota. At the phylum level, abundance of Firmicutes was reduced while abundance of Proteobacteria and Fusobacteria increased in Ss+ subjects. At the family level, there was an increase of Clostridiaceae, Streptococcaceae, Desulfovibrionaceae and Enterobacteriaceae in the Ss+ group. A previous study demonstrated that these families were associated with trimethylamine (TMA) production [38]. The high abundances of families Enterobacteriaceae, Clostridiaceae and Veillonellaceae in Ss+ subjects are in agreement with a previous study [39]. These bacteria are associated with increasing the fecal pH [39] to a level where most opportunistic bacterial pathogens prefer to grow [40]. The change of microbiota composition that we observed may influence the environment in the gut, suggesting that *S. stercoralis* infection may influence the microbiota and modulate the pH of the gut environment.

Forty-two genera showed contrasting abundances between the two groups. Interestingly, there were significant increases of pathogenic bacteria including *Escherichia-Shigella, Rothia* and *Aggregatibacter* and some increase of *Actinomyces, Streptococcus* and *Haemophilus* in CKD patients infected with *S. stercoralis* compared to uninfected controls. In contrast, significant reduction of some SCFA-producing bacteria, such as *Anaerostipes* and *Coprococcus*_1 and some decrease of *Akkermansia, Eubacterium rectale_group* and *Eubacterium hallii_group* were noted in the Ss+ group. Specifically, abundance of the genus *Escherichia-Shigella* is known to be positively correlated with uremic toxins such as trimethylamine-N-oxide and indoxyl sulfate [41]. Our results demonstrated a significant inverse correlation of *Escherichia-Shigella* with the estimated glomerular filtration rate (eGFR r = −0.37, *P* = 0.018). eGFR is one criterion for diagnosis and staging of CKD. Thus, high abundance of *Escherichia-Shigella* was correlated with low eGFR value and higher CKD stage (Figure 6). In addition, *Enterobacteriaceae* and *E. coli* are markedly more abundant in individuals with impaired kidney function as demonstrated previously [42], highlighting that there is an association between the genus *Escherichia-Shigella* and CKD with concurrent *S. stercoralis* infection.

Interestingly, the genus *Anaerostipes* was less abundant in CKD patients with concurrent *S. stercoralis* infection than in those without. Members of this genus are Gram-variable, obligate anaerobes which produce acetate, butyrate and lactate from glucose fermentation [43]. Our findings were partially consistent with those of a previous study, which found that *Anaerostipes* had low relative abundance in CKD in an animal model and noted that this genus was negatively correlated with amount of intestinal urea. Nephrectomized mice with low levels of *Anaerostipes* exhibited negative effects on kidney parameters (BUN and creatinine) [44]. However, we found no correlation here between *Anaerostipes* and kidney parameters. This may be due to the limited sample size. Recent research indicated that elevated levels of *Anaerostipes* led to increased eGFR and improvement in renal function [45]. Furthermore, the relative abundance of the genus *Anearostipes* was markedly reduced in nonsurvivors with end-stage kidney disease (ESKD) [46]. Specifically, we found that one species, *Anaerostipes hadrus*, an important microbe in maintaining intestinal metabolic balance [47], was significantly reduced in CKD with concurrent *S. stercoralis* infection.

A previous study using a rat model revealed that increased levels of *Rothia* were positively associated with creatinine levels in acute kidney injury and with severity of kidney damage [48]. *Rothia* spp. are Gram-positive cocco-bacilli that cause a wide range of serious infections, especially in immunocompromised hosts. *Rothia* is often identified in blood cultures from patients with bacteremia [49]. In this study, we found a positive correlation between *Rothia* and the genus *Streptococcus* (r = 0.47, P = 0.002). A similar relationship was observed in a recent study, which found that the log-ratio between the presence of the genera *Rothia* and *Streptococcus* genera was the best predictor of creatinine level [48].

The main SCFA-producing bacteria in humans [50] include *Faecalibacterium prausnitzii* (butyrate-producing bacteria in the phylum Firmicutes), *Eubacterium rectale* and *E. hallii* (family Lachnospiraceae) and *Anaerostipes* spp. (sugar/lactate-utilizing bacteria producing butyrate from lactate and acetate) were all less abundant in *S. stercoralis* infection. A previous study observed a decrease in the level of serum SCFAs in CKD patients and an inverse correlation between butyrate level and renal function [51]. Our study suggested that SCFA-producing bacteria are depleted in CKD with concurrent *S. stercoralis* infection, which may affect CKD progression. However, further studies are needed to confirm this association.

Our study has strengths and limitations. An important strength of this study is that we used groups that were pair-matched for sex, age and biochemical factors. However, we did not obtain data for concentrations of uremic toxins (TMAO and IS) or for amounts of SCFAs. Moreover, we did not record the clinical manifestation of *S. stercoralis* infection in CKD patient so we were not able to show the association between some pathogenic bacteria and *S. stercoralis* infection in CKD. In addition, the sample size in this study was small due to the fact that participant who infected with single *S. stercoralis* and got CKD which was limited. These affect the estimation of microbiome alpha diversity and other statistical power in analyses.

## Conclusions

This study suggests that *S. stercoralis* infection reduces the diversity of the gut microbiota in CKD patients. An increased abundance of harmful bacteria and reduction of some SCFA-producing bacteria in *S. stercoralis* infection were found. In addition, the abundance of members of the genus *Escherichia-Shigella* was significantly and inversely correlated with eGFR levels. Significant elevation of members of this genus in CKD patients with *S. stercoralis* infection may indicate potential diagnostic markers for CKD in *S. stercoralis*-endemic areas. Thus, we suggest that these changes in the composition of the gut microbiome in *S. stercoralis* infection may result in disruption of the gut barrier structure and absorption of harmful products that can contribute to toxicity, inflammation and malnutrition, contributing to CKD progression. Future metabolomics studies are required to unravel the relationship between CKD and *S. stercoralis* infection.

## Data Availability

The Illumina datasets for the profile of gut microbiota (Raw data) are available from Mendeley Data (DOI: 10.17632/hvbvrtc34x.1)

## Funding

This study was supported by Research Program (RP65-2-001), CKDNET (grant no. CKDNET2559007) and Invitation research grant, Faculty of Medicine Research Grant (IV63136), Khon Kaen University. S.A. acknowledged research funds from Thailand Science Research and Innovation (TSRI), through Program Management Unit for Competitiveness (PMUC), number C10F630030.

## Author contributions

**Conceptualization:** Nuttanan Hongsrichan, Porntip Pinlaor, Somchai Pinlaor.

**Formal analysis:** Nguyen Thi Hai, Nuttanan Hongsrichan, Kitti Intuyod, Manachai Yingklang, Apisit Chaidee, Thatsanapong Pongking, Porntip Pinlaor, Sirirat Anutrakulchai, Ubon Cha’on, Somchai Pinlaor.

**Funding acquisition:** Nguyen Thi Hai, Sirirat Anutrakulchai, Ubon Cha’on, Somchai Pinlaor.

**Investigation:** Nguyen Thi Hai, Nuttanan Hongsrichan, Kitti Intuyod, Porntip Pinlaor, Manachai Yingklang, Apisit Chaidee, Thatsanapong Pongking, Sirirat Anutrakulchai, Ubon Cha’on, Somchai Pinlaor.

**Project administration:** Somchai Pinlaor.

**Resources:** Sirirat Anutrakulchai, Ubon Cha’on, Somchai Pinlaor.

**Supervision:** Nuttanan Hongsrichan, Porntip Pinlaor, Somchai Pinlaor.

**Visualization:** Nguyen Thi Hai, Thatsanapong Pongking, Kitti Intuyod, Apisit Chaidee.

**Writing - original draft:** Nguyen Thi Hai.

**Writing - review & editing:** Nguyen Thi Hai, Nuttanan Hongsrichan, Somchai Pinlaor.

## Acknowledgments

N.T.H. thanks the scholarship under the Doctoral Training Program from Research Affairs, Medicine Faculty and Graduate School, Khon Kaen University, Thailand. We would like to thank to all of the people who voluntarily participated in this study. We also acknowledge Prof. David Blair from Publication Clinic KKU, Thailand, for his comments and editing the manuscript.

*   Received March 7, 2022.
*   Revision received March 7, 2022.
*   Accepted March 9, 2022.


*   © 2022, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/)

## References

1.  1.Nallu A, Sharma S, Ramezani A, Muralidharan J, Raj D. Gut microbiome in chronic kidney disease: challenges and opportunities. Transl Res. 2017;179:24–37. pmid: 27187743.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.trsl.2016.04.007&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27187743&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F09%2F2022.03.07.22271995.atom) 

2.  2.Hida M, Aiba Y, Sawamura S, Suzuki N, Satoh T, Koga Y. Inhibition of the accumulation of uremic toxins in the blood and their precursors in the feces after oral administration of Lebenin, a lactic acid bacteria preparation, to uremic patients undergoing hemodialysis. Nephron. 1996;74(2):349–55. pmid: 8893154.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1159/000189334&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8893154&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F09%2F2022.03.07.22271995.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1996VL83900017&link_type=ISI) 

3.  3.Wang F, Jiang H, Shi K, Ren Y, Zhang P, Cheng S. Gut bacterial translocation is associated with microinflammation in end-stage renal disease patients. Nephrology (Carlton). 2012;17(8):733–8. pmid: 22817644.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1440-1797.2012.01647.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22817644&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F09%2F2022.03.07.22271995.atom) 

4.  4.Vaziri ND, Wong J, Pahl M, Piceno YM, Yuan J, DeSantis TZ, et al. Chronic kidney disease alters intestinal microbial flora. Kidney Int. 2013;83(2):308–15. pmid: 22992469.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ki.2012.345&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22992469&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F09%2F2022.03.07.22271995.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000314621500020&link_type=ISI) 

5.  5.De Angelis M, Montemurno E, Piccolo M, Vannini L, Lauriero G, Maranzano V, et al. Microbiota and metabolome associated with immunoglobulin A nephropathy (IgAN). PLoS One. 2014;9(6):e99006. pmid: 24922509.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0099006&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24922509&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F09%2F2022.03.07.22271995.atom) 

6.  6.Coresh J, Selvin E, Stevens LA, Manzi J, Kusek JW, Eggers P, et al. Prevalence of chronic kidney disease in the United States. Jama. 2007;298(17):2038–47. pmid: 17986697.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.298.17.2038&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17986697&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F09%2F2022.03.07.22271995.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000250718200024&link_type=ISI) 

7.  7.Ramezani A, Raj DS. The gut microbiome, kidney disease, and targeted interventions. J Am Soc Nephrol. 2014;25(4):657–70. pmid: 24231662.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiam5lcGhyb2wiO3M6NToicmVzaWQiO3M6ODoiMjUvNC82NTciO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wMy8wOS8yMDIyLjAzLjA3LjIyMjcxOTk1LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

8.  8.Wong J, Piceno YM, DeSantis TZ, Pahl M, Andersen GL, Vaziri ND. Expansion of urease-and uricase-containing, indole-and p-cresol-forming and contraction of short-chain fatty acid-producing intestinal microbiota in ESRD. Am J Nephrol. 2014;39(3):230–7. pmid: 24643131.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1159/000360010&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24643131&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F09%2F2022.03.07.22271995.atom) 

9.  9.Jakobsson HE, Jernberg C, Andersson AF, Sjolund-Karlsson M, Jansson JK, Engstrand L. Short-term antibiotic treatment has differing long-term impacts on the human throat and gut microbiome. PLoS One. 2010;5(3):e9836. pmid: 20352091.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0009836&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20352091&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F09%2F2022.03.07.22271995.atom) 

10. 10.Kalantar-Zadeh K, Kopple JD, Deepak S, Block D, Block G. Food intake characteristics of hemodialysis patients as obtained by food frequency questionnaire. J Ren Nutr. 2002;12(1):17–31. pmid: 11823990.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1053/jren.2002.29598&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11823990&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F09%2F2022.03.07.22271995.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000173838600003&link_type=ISI) 

11. 11.Fouque D, Cruz Casal M, Lindley E, Rogers S, Pancirova J, Kernc J, et al. Dietary trends and management of hyperphosphatemia among patients with chronic kidney disease: an international survey of renal care professionals. J Ren Nutr. 2014;24(2):110–5. pmid: 24447438.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24447438&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F09%2F2022.03.07.22271995.atom) 

12. 12.Lunyera J, Mohottige D, Von Isenburg M, Jeuland M, Patel UD, Stanifer JW. CKD of uncertain etiology: a systematic review. Clin J Am Soc Nephrol. 2016;11(3):379–85. pmid: 26712810.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiY2xpbmphc24iO3M6NToicmVzaWQiO3M6ODoiMTEvMy8zNzkiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wMy8wOS8yMDIyLjAzLjA3LjIyMjcxOTk1LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

13. 13.Leung JM, Graham AL, Knowles SCL. Parasite-Microbiota interactions with the vertebrate gut: synthesis through an ecological lens. Front Microbiol. 2018;9:843. pmid: 29867790.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29867790&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F09%2F2022.03.07.22271995.atom) 

14. 14.Cha’on U, Wongtrangan K, Thinkhamrop B, Tatiyanupanwong S, Limwattananon C, Pongskul C, et al. CKDNET, a quality improvement project for prevention and reduction of chronic kidney disease in the Northeast Thailand. BMC Public Health. 2020;20(1):1299. pmid: 32854662.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32854662&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F09%2F2022.03.07.22271995.atom) 

15. 15.De Souza JN, Soares B, Goes LL, Lima CS, Barreto N, Jacobina BS, et al. Case report: Strongyloides stercoralis hyperinfection in a patient with HTLV-1: a case report of an infection with filariform and rhabditiform larvae, eggs, and free-living adult females. Am J Trop Med Hyg. 2018. pmid: 30277207.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30277207&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F09%2F2022.03.07.22271995.atom) 

16. 16.Winnicki W, Eder M, Mazal P, Mayer FJ, Sengölge G, Wagner L. Prevalence of Strongyloides stercoralis infection and hyperinfection syndrome among renal allograft recipients in Central Europe. Scientific Reports. 2018;8(1):15406.
    
    

17. 17.Qu TT, Yang Q, Yu MH, Wang J. A fatal Strongyloides stercoralis hyperinfection syndrome in a patient with chronic kidney disease: a case report and literature review. Medicine (Baltimore). 2016;95(19):e3638. pmid: 27175679.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27175679&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F09%2F2022.03.07.22271995.atom) 

18. 18.Krolewiecki A, Nutman TB. Strongyloidiasis: a neglected tropical disease. Infect Dis Clin North Am. 2019;33(1):135–51. pmid: 30712758.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30712758&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F09%2F2022.03.07.22271995.atom) 

19. 19.Jenkins TP, Peachey LE, Ajami NJ, MacDonald AS, Hsieh MH, Brindley PJ, et al. Schistosoma mansoni infection is associated with quantitative and qualitative modifications of the mammalian intestinal microbiota. Sci Rep. 2018;8(1):12072. pmid: 30104612.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41598-018-30412-x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30104612&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F09%2F2022.03.07.22271995.atom) 

20. 20.Pane S, Sacco A, Iorio A, Romani L, Putignani L. Strongyloides stercoralis infestation in a child: how a nematode can affect gut microbiota. International Journal of Molecular Sciences. 2021;22(4):2131. pmid: doi:10.3390/ijms22042131.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/ijms22042131&link_type=DOI) 

21. 21.Zoccali C, Vanholder R, Massy ZA, Ortiz A, Sarafidis P, Dekker FW, et al. The systemic nature of CKD. Nat Rev Nephrol. 2017;13(6):344–58. pmid: 28435157.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrneph.2017.52&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28435157&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F09%2F2022.03.07.22271995.atom) 

22. 22.Inker LA, Astor BC, Fox CH, Isakova T, Lash JP, Peralta CA, et al. KDOQI US commentary on the 2012 KDIGO clinical practice guideline for the evaluation and management of CKD. Am J Kidney Dis. 2014;63(5):713–35. pmid: 24647050.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1053/j.ajkd.2014.01.416&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24647050&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F09%2F2022.03.07.22271995.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000335456600412&link_type=ISI) 

23. 23.Kaewrat W, Sengthong C, Yingklang M, Intuyod K, Haonon O, Onsurathum S, et al. Improved agar plate culture conditions for diagnosis of Strongyloides stercoralis. Acta Trop. 2020;203:105291. pmid: 31816322.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31816322&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F09%2F2022.03.07.22271995.atom) 

24. 24.Fellows R, Denizot J, Stellato C, Cuomo A, Jain P, Stoyanova E, et al. Microbiota derived short chain fatty acids promote histone crotonylation in the colon through histone deacetylases. Nat Commun. 2018;9(1):105. pmid: 29317660.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-017-02651-5&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29317660&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F09%2F2022.03.07.22271995.atom) 

25. 25.Itthitaetrakool U, Pinlaor P, Pinlaor S, Chomvarin C, Dangtakot R, Chaidee A, et al. Chronic Opisthorchis viverrini infection changes the liver microbiome and promotes Helicobacter growth. PLOS ONE. 2016;11(11):e0165798.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0165798&link_type=DOI) 

26. 26.Magoč T, Salzberg SL. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics. 2011;27(21):2957–63. pmid: 21903629.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btr507&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21903629&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F09%2F2022.03.07.22271995.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000296099300005&link_type=ISI) 

27. 27.Bokulich NA, Subramanian S, Faith JJ, Gevers D, Gordon JI, Knight R, et al. Quality-filtering vastly improves diversity estimates from Illumina amplicon sequencing. Nat Methods. 2013;10(1):57–9. pmid: 23202435.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nmeth.2276&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23202435&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F09%2F2022.03.07.22271995.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000312810100036&link_type=ISI) 

28. 28.Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nature methods. 2010;7(5):335–6. pmid: 20383131.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nmeth.f.303&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20383131&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F09%2F2022.03.07.22271995.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000277175100003&link_type=ISI) 

29. 29.Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013;41(Database issue):D590–6. pmid: 23193283.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gks1219&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23193283&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F09%2F2022.03.07.22271995.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000312893300084&link_type=ISI) 

30. 30.Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011;27(16):2194–200. pmid: 21700674.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btr381&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21700674&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F09%2F2022.03.07.22271995.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000293620800004&link_type=ISI) 

31. 31.Haas BJ, Gevers D, Earl AM, Feldgarden M, Ward DV, Giannoukos G, et al. Chimeric 16S rRNA sequence formation and detection in Sanger and 454-pyrosequenced PCR amplicons. Genome Res. 2011;21(3):494–504. pmid: 21212162.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiZ2Vub21lIjtzOjU6InJlc2lkIjtzOjg6IjIxLzMvNDk0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDMvMDkvMjAyMi4wMy4wNy4yMjI3MTk5NS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

32. 32.Edgar RC. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nature Methods. 2013;10(10):996–8.
    
    

33. 33.DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, et al. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Applied and environmental microbiology. 2006;72(7):5069–72. pmid: 16820507.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYWVtIjtzOjU6InJlc2lkIjtzOjk6IjcyLzcvNTA2OSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzAzLzA5LzIwMjIuMDMuMDcuMjIyNzE5OTUuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

34. 34.Wang Q, Garrity GM, Tiedje JM, Cole JR. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Applied and environmental microbiology. 2007;73(16):5261–7. pmid: 17586664.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYWVtIjtzOjU6InJlc2lkIjtzOjEwOiI3My8xNi81MjYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDMvMDkvMjAyMi4wMy4wNy4yMjI3MTk5NS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

35. 35.Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Research. 2004;32(5):1792–7.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkh340&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15034147&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F09%2F2022.03.07.22271995.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000220487200025&link_type=ISI) 

36. 36.Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, et al. Metagenomic biomarker discovery and explanation. Genome Biology. 2011;12(6):R60.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/gb-2011-12-6-r60&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21702898&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F09%2F2022.03.07.22271995.atom) 

37. 37.Formenti F, La Marca G, Perandin F, Pajola B, Romano M, Santucci B, et al. A diagnostic study comparing conventional and real-time PCR for Strongyloides stercoralis on urine and on faecal samples. Acta Tropica. 2019;190:284–7.
    
    

38. 38.Al-Obaide MAI, Singh R, Datta P, Rewers-Felkins KA, Salguero MV, Al-Obaidi I, et al. Gut microbiota-dependent trimethylamine-N-oxide and serum biomarkers inpatients with T2DM and advanced CKD. J Clin Med. 2017;6(9). pmid: 28925931.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28925931&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F09%2F2022.03.07.22271995.atom) 

39. 39.Henrick BM, Hutton AA, Palumbo MC, Casaburi G, Mitchell RD, Underwood MA, et al. Elevated fecal pH indicates a profound change in the breastfed infant gut microbiome due to reduction of <em>Bifidobacterium</em> over the past century. mSphere. 2018;3(2):e00041–18.
    
    

40. 40.Duar RM, Kyle D, Casaburi G. Colonization Resistance in the Infant Gut: the role of B. infantis in reducing pH and preventing pathogen growth. High Throughput. 2020;9(2). pmid: 32230716.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32230716&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F09%2F2022.03.07.22271995.atom) 

41. 41.Tomlinson JAP, Wheeler DC. The role of trimethylamine N-oxide as a mediator of cardiovascular complications in chronic kidney disease. Kidney International. 2017;92(4):809–15.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28807612&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F09%2F2022.03.07.22271995.atom) 

42. 42.Gryp T, Huys GRB, Joossens M, Van Biesen W, Glorieux G, Vaneechoutte M. Isolation and quantification of uremic toxin precursor-generating gut bacteria in chronic kidney disease patients. Int J Mol Sci. 2020;21(6). pmid: 32183306.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32183306&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F09%2F2022.03.07.22271995.atom) 

43. 43.Schwiertz A, Hold GL, Duncan SH, Gruhl B, Collins MD, Lawson PA, et al. Anaerostipes caccae gen. nov., sp. nov., a new saccharolytic, acetate-utilising, butyrate-producing bacterium from human faeces. Syst Appl Microbiol. 2002;25(1):46–51. pmid: 12086188.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1078/0723-2020-00096&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12086188&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F09%2F2022.03.07.22271995.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000176285200006&link_type=ISI) 

44. 44.Chaves LD, McSkimming DI, Bryniarski MA, Honan AM, Abyad S, Thomas SA, et al. Chronic kidney disease, uremic milieu, and its effects on gut bacterial microbiota dysbiosis. Am J Physiol Renal Physiol. 2018;315(3):F487–f502. pmid: 29693447.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1152/ajprenal.00092.2018&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29693447&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F09%2F2022.03.07.22271995.atom) 

45. 45.Mazidi M, Shekoohi N, Covic A, Mikhailidis DP, Banach M. Adverse impact of Desulfovibrio spp. and beneficial role of Anaerostipes spp. on renal function: insights from a mendelian randomization analysis. Nutrients. 2020;12(8):2216. pmid: doi:10.3390/nu12082216.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/nu12082216&link_type=DOI) 

46. 46.Lin TY, Wu PH, Lin YT, Hung SC. Gut dysbiosis and mortality in hemodialysis patients. NPJ Biofilms Microbiomes. 2021;7(1):20. pmid: 33658514.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33658514&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F09%2F2022.03.07.22271995.atom) 

47. 47.Kant R, Rasinkangas P, Satokari R, Pietila TE, Palva A. Genome sequence of the butyrate-producing anaerobic bacterium Anaerostipes hadrus PEL 85. Genome Announc. 2015;3(2). pmid: 25838483.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MjoiZ2EiO3M6NToicmVzaWQiO3M6MTM6IjMvMi9lMDAyMjQtMTUiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wMy8wOS8yMDIyLjAzLjA3LjIyMjcxOTk1LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

48. 48.Andrianova NV, Popkov VA, Klimenko NS, Tyakht AV, Baydakova GV, Frolova OY, et al. Microbiome-metabolome signature of acute kidney injury. Metabolites. 2020;10(4). pmid: 32260384.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32260384&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F09%2F2022.03.07.22271995.atom) 

49. 49.Ramanan P, Barreto JN, Osmon DR, Tosh PK. Rothia bacteremia: a 10-year experience at Mayo Clinic, Rochester, Minnesota. J Clin Microbiol. 2014;52(9):3184–9. pmid: 24951810.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamNtIjtzOjU6InJlc2lkIjtzOjk6IjUyLzkvMzE4NCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzAzLzA5LzIwMjIuMDMuMDcuMjIyNzE5OTUuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

50. 50.Parada Venegas D, De la Fuente MK, Landskron G, González MJ, Quera R, Dijkstra G, et al. Short chain fatty acids (SCFAs)-mediated gut epithelial and immune regulation and its relevance for inflammatory bowel diseases. Frontiers in Immunology. 2019;10(277).
    
    

51. 51.Wang S, Lv D, Jiang S, Jiang J, Liang M, Hou F, et al. Quantitative reduction in short-chain fatty acids, especially butyrate, contributes to the progression of chronic kidney disease. Clin Sci (Lond). 2019;133(17):1857–70. pmid: 31467135.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6OToicHBjbGluc2NpIjtzOjU6InJlc2lkIjtzOjExOiIxMzMvMTcvMTg1NyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzAzLzA5LzIwMjIuMDMuMDcuMjIyNzE5OTUuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

52. 52.Morimoto J, Kaneoka H, Sasatomi Y, Sato YN, Murata T, Ogahara S, et al. Disseminated strongyloidiasis in nephrotic syndrome. Clin Nephrol. 2002;57(5):398–401. pmid: 12036202.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12036202&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F09%2F2022.03.07.22271995.atom)