Table 1 Differentially abundant KEGG level 2 pathways at each sampling time based on PICRUSt

10.03 ± 0.75A

9.99 ± 0.54A

9.80 ± 0.63AB

8.94 ± 1.22B

0.74 ± 0.13AB

0.79 ± 0.09A

0.78 ± 0.11AB

0.66 ± 0.13B

2.00 ± 0.62A

1.92 ± 0.37AB

1.74 ± 0.39AB

1.41 ± 0.56B

3.70 ± 0.54A

3.58 ± 0.37AB

3.62 ± 0.5A

3.01 ± 0.77B

5.95 ± 0.24AB

5.82 ± 0.2A

5.89 ± 0.2A

6.17 ± 0.34B

2.68 ± 0.13A

2.73 ± 0.11AB

2.79 ± 0.1AB

2.88 ± 0.2B

2.37 ± 0.12A

2.32 ± 0.13AB

2.29 ± 0.12AB

2.19 ± 0.21B

1.71 ± 0.12A

1.61 ± 0.07B

1.54 ± 0.09B

1.58 ± 0.09B

1.82 ± 0.32A

1.8 ± 0.24A

1.7 ± 0.25AB

1.41 ± 0.47B

4.15 ± 0.34A

4.31 ± 0.36AB

4.4 ± 0.33AB

4.67 ± 0.5B

9.24 ± 1.05A

9.24 ± 0.9A

9.6 ± 1.04A

11.01 ± 1.85B

1.61 ± 0.23A

1.54 ± 0.16A

1.44 ± 0.18AB

1.23 ± 0.36B

2.46 ± 0.17A

2.71 ± 0.14B

2.7 ± 0.23B

2.44 ± 0.17A

6.61 ± 1.3A

6.61 ± 1.37A

6.96 ± 1.39A

8.88 ± 2.54B

2.61 ± 0.32A

2.39 ± 0.26AB

2.12 ± 0.24B

2.07 ± 0.48B