Resilience of human gut microbial communities for the long stay with multiple dietary shifts.

We read with great interest the article by Bennet et al 1 that investigated the impact of dietary shifts on gut microbiota. Together with several recent studies,2–5 they have shown that short-term dietary shifts can alter the composition of gut microbiota for both patients and healthy humans. However, these works drawn the conclusion with only a few time points and the subjects under investigation had undergone just one dietary shift. The dynamic patterns of microbial communities across longer time scales with multiple dietary shifts remain unclear.

In this study, we recruited a Chinese volunteer team (VT) composing of 10 people who departed from Beijing, conducted a long stay of 6 months in Trinidad and Tobago (TAT) and returned to Beijing. Their faecal samples (188 samples), together with detailed dietary information, were collected using a high-density longitudinal sampling strategy (19 time points on average for VT members). We partitioned the whole longitudinal study into six phases: T1 represents pre-travel phase (20 samples), T2 (28 samples), T3 (60 samples) and T4 (21 samples) represent three time-slots when VT stayed in TAT, and T5 (35 samples) and T6 (20 samples) represent two time-slots after VT returned to Beijing. Meanwhile, we also collected samples from healthy natives of Beijing (BJN, 57 samples), TAT healthy natives (TTN, 28 samples), TAT patients (TTP, six samples) and TAT Chinese (TTC, eight samples) as control datasets. Finally, we sequenced the V4 hypervariable region of microbial 16S rRNA genes for 287 faecal samples from 41 individuals and analysed the high-quality reads with QIIME.6

Interestingly, we found that gut microbial communities are bidirectionally plastic and resilient across the long stay with multiple dietary shifts. First of all, BJN and TTN exhibited different patterns of microbial communities (figure 1A). The microbial communities of VT members, however, shifted from those similar to BJN towards TTN pattern accompanying their stays in TAT, and reverted to their original patterns within 1 month after VT returned to Beijing (figure 1B–F). Such plastic pattern is also verified by DeBlur7 analysis. Besides, though we have found that locations and populations have strong effects on differentiation of samples (figure 1G,H), the dynamic changes for each member of VT showed specific trends (figure 1I,J), indicating the possible enterotypes-dependent plastic patterns among VT members. Additionally, the relative abundances of phylum Firmicutes and Bacteroidetes have shown a strong resilience along the timeline, and they are negatively correlated along the timeline (figure 1K,L). Similarly, the relative abundance of Proteobacteria and Actinobacteria also exhibited plastic patterns (figure 1M,N). By tracing and comparing the common operational taxonomic units (OTUs) shared by at least 10% of VT members, we found that individual OTUs within Firmicutes, Bacteroidetes, Proteobacteria and Actinobacteria have distinctive temporal dynamics during VT’s long stay (figure 1O–R).

Figure 1

Pattern of human gut microbial community during the long stay with multiple dietary shifts. The dynamics of human gut microbial communities are shown through comparing different groups/phases of samples, including (A) BJN and TTN groups, (B) BJN, TTN and T1 groups, (C) BJN, TTN, T1 and T2-T4 groups, (D) BJN, TTN, T1, T2-T4 and T5-T6 groups, (E) BJN, TTN, T1, T2-T4, T5-T6 and TTC groups, (F) BJN, TTN, T1, T2-T4, T5-T6, TTC and TTP groups. The large arrows indicate the progressive add-in of groups for comparison. The profile of all faecal samples labelled by (G) collected locations, and (H) populations (Chinese or TATs). In (H), the label ’Chinese' represents all samples from BJN, VT and TTC, while the label ’TATs' represents all samples from TTN and TTP. The dynamic changes of gut microbiota for two members of the VT, including (I) VT5, and (J) VT6. The curves with arrows in (I) and (J) represent the directions of changes of microbial communities. Panels (K), (L), (M) and (N) show the dynamics of the average relative abundance of Firmicutes, Bacteroidetes, Proteobacteria and Actinobacteria along the temporal axis. Results were visualised using PCoA of weighted UniFrac distances of 16S rRNA amplicon profiles. OTUs shared by at least 10% VT samples during same phase (green bars) are tracked using Sankey plots for (O) Firmicutes, (P) Bacteroidetes, (Q) Proteobacteria and (R) Actinobacteria. The heights of bars are proportional to the number of OTUs, and phases are arranged in chronological order. The lines represent the transfer of OTUs between phases and are coloured according to the first phase of their appearance. The abbreviations used in this figure are listed as follows: BJN, Beijing natives; OTUs, operational taxonomic units; PCoA, principal coordinates analysis; TTN, Trinidad and Tobago natives; T1-T6, six phases; TTC, Trinidad and Tobago Chinese; TTP, Trinidad and Tobago patient; VT, volunteer team.

The bi-directional plasticity of human microbial community was largely mediated by dietary shifts. During their stay in TAT, VT members consumed local foods (online supplementary table 1) mainly including fish and seafood, dairy products, and refined grains, which were considerably different from those in Beijing (online supplementary table 2). Significant correlation was observed between dietary shifts and the dynamics of gut microbial community structures (p<0.05, Mantel’s test, online supplementary table 3 and 4).

The bi-directional plasticity of gut microbiome shown in this study could elicit effects on clinical practice for gut microbiome-related diseases in at least two ways: First, it indicated that diagnostics should incorporate dietary habits and travel records. Second, it suggested that due to the strong resilience of gut microbial communities, clinical practices, such as faecal microbiota transplantation and antibiotic treatment, should monitor therapeutic effects across a longer time scale.8