Microbial community changes in the digestive tract of the clam Meretrix petechialis in response to Vibrio parahaemolyticus challenge*

Hongxia WANG , Xin YUE , Jiajia YU Rui WANG Shuangshuang TENG , Jun FANG , Baozhong LIU ,

1 CAS Key Laboratory of Experimental Marine Biology, Institute of Oceanology, Center for Ocean Mega-Science, Chinese Academy of Sciences, Qingdao 266071, China

2 Laboratory for Marine Biology and Biotechnology, Qingdao National Laboratory for Marine Science and Technology, Qingdao 266000, China

3 Zhejiang Mariculture Research Institute, Wenzhou 325005, China

4 Zhejiang Key Laboratory of Exploitation and Preservation of Coastal Bio-resource, Wenzhou 325005, China

Abstract Disease in clams frequently occurred over the last decade and has become a serious threat to the clam aquaculture industry and natural stocks. Mass clam mortality events were reported to be associated with the presence of opportunistic pathogen vibrio. However, the complexity of infection that occurs in the natural environment remains poorly understood. In this study, we smulated a natural disease outbreak by vibrio immersion infection to study the diversity and dynamics of microbiota in the digestive tract of clam Meretrix petechialis during the infection process. Dramatic changes in operational taxonomic unit richness and phylum composition of the bacterial communities were observed during pathogen invasion. In addition, we investigated the potential relationship between microbiota dynamics and host status during disease progression. Results reveal that, at the end stage of vibrio infection, interindividual variation in the digestive tract microbiota increased, as did the difference in individual health status. The moribund clams displayed signs of microbial community shifts to low diversity, and the microbial community was characterized by mass proliferation of a few operational taxonomic units.

Keyword: Meretrix petechialis; vibrio challenge; 16S ribosomal RNA; bacterial community; digestive tract

1 INTRODUCTION

According to the latest the State of World Fisheries and Aquaculture (SOFIA) report, bivalves are one of the most important food products for the aquaculture industry, and worldwide production (mainly oysters, mussels, clams, and scallops) reached a volume of 15 t with an economic value of more than $18 billion (FishStatJ, Food and Agriculture Organization of the United Nations (FAO) http://www.fao.org/fishery/statistics/global-aquaculture-production/en). China is one of the main bivalve-producing countries, and the clam culture industry has become very important for many regions in China. However, some infectious diseases have frequently occurred over the last decade and have become the primary limiting factors for the viability and development of the bivalve aquaculture industry. For example, Vibrio parahaemolyticus is an opportunistic pathogen that usually causes mass mortality of the cultured clam Meretrix petechialis in summer. This gram-negative bacterium isolated from diseased clams was demonstrated to be the main pathogenic by artificial infection, and a correlation exists between adult clam mortality and an increased concentration of V. parahaemolyticus (Yue et al., 2010).

As invertebrates, clams largely depend on innate immunity to defend against and resist infectious agents. The hepatopancreas is part of the digestive tract and an important member of the immune system in which epithelial cells produce the majority of immune molecules (Rőszer, 2014). In addition, epithelia lining the digestive tract of bivalves phagocytize biotic and abiotic particles and colloids, so it has dual function in nutrition and defense by promoting nutrient absorption and by keeping pathogens in check to improve disease resistance (Allam and Espinosa, 2016). Additionally, the digestive tract of bivalves has abundant microorganisms, including pathogens and commensals from the surrounding environment. For example, it was reported that oysters are capable of filtering over 10 L/(h∙g) dry tissue (Dupuy et al., 2000), which equated to over 25 000 microbes/s for a 1-g oyster in a modest 104microbes/mL seawater (Allam and Espinosa, 2016). Therefore, the entire digestive tract microflora may have crucial effects on the health status of the host during vibrio infection.

It is becoming clear that diseases could substantially change the intestinal microbiota composition in both vertebrates and invertebrates (Yang et al., 2017; Xiong et al., 2018; Yu et al., 2019b). Studies on the microbiota-disease relationship in shrimp intestines revealed that dysbiosis in unhealthy intestinal microbiota preceded the appearance of physical symptoms of the disease, which emphasized that microorganisms are the first to respond to changes and can augment disease progression (Xiong, 2018). In addition, studies identified that the degree of change in intestinal bacterial communities of shrimp is closely associated with the disease severity (Xiong et al., 2015; Chen et al., 2017). Currently, research on bacterial community alterations in bivalves with vibrio disease is still rare and has mostly focused on hemolymph (Wendling et al., 2014; Lokmer and Wegner, 2015; Petton et al., 2015). Lokmer and Wegner (2015) previously reported that the microbial community of oyster hemolymph could be significantly affected by temperature stress, and more stable communities had greater resistance to invasion by pathogens. However, little is known about the role of the pathogen vibrio and the dynamics of how the microbial communities of the digestive tract change during disease progression when clams are infected with a high concentration of vibrio.

Increasing evidence has revealed complex interactions among pathogens, hosts and commensal microbiota (Boutin et al., 2013). The rapid development of high-throughput sequencing technologies and bioinformatic analyses facilitates indepth characterization of complex and varied microbiota (Caporaso et al., 2010; Boutin et al., 2013) and further development a broader understanding of the complexity of host-bacteria interactions. In this study, a 16S rRNA sequencing approach was used to investigate the dynamic changes of digestive tract bacterial communities in the clam M. petechialis subjected to bacterial challenge with V. parahaemolyticus. The microbiological data obtained in this study help us to understand the role and composition of host-associated microbial communities that are affected by vibrio infection, which is important for diagnosing, predicting, and preventing disease outbreaks in clam aquaculture.

2 MATERIAL AND METHOD

2.1 Experimental clams

In total, 250 clams ( M. petechialis) with a mean body weight of 30.76±8.01 g were used for the experiment. The clams were caught from the cultured population and brought into the laboratory. Before vibrio challenge, they were acclimated for 2 weeks in a holding system at 26 °C with salinity at 22.0±1.0 and continuously aerated seawater, and fed Isochrysis galbana. This holding system maintained water quality with a filtration unit.

2.2 Challenge experiment

At the beginning of the challenge experiment, the digestive tracts of 30 acclimated clams were arbitrarily sampled and used to assay the initial vibrio load. By colony-forming unit (CFU) counts, no vibrio load was detected in 83.3% of individuals (25/30), and a vibrio load between 0.08-1.26 CFU/mg tissue was detected in the other five clams (the detailed sampling and vibrio load assay methods are described in Section 2.3). Therefore, it could be deduced that there was seldom vibrio infection in the clams before the artificial infection. Vibrio challenge was conducted according to Liang et al. (2017). Briefly, approximately 150 clams were arbitrarily collected and divided into five groups. Each group of 28-32 clams was put in a single cage, and all clams from the five cages were reared in a water tank with V. parahaemolyticus strain MM21 (Yue et al., 2010) added to the tank at a final concentration of 1×107CFU/mL at 30 °C. Seawater mixed with vibrio was renewed once per day. Three groups were used to monitor daily mortality and the other two groups were used for sample collection. After 4 days (i.e., on day 5), the artificial vibrio challenge was stopped, and the clams were reared in fresh seawater again. Meanwhile, the remaining clams remained in a fresh seawater tank and were classified as the control group, which was used to confirm that mortality was due to pathogen exposure. Clams were checked every 6 h per day, dead or moribund individuals were collected and recorded. The daily mortality was calculated based on three replicates.

2.3 Sampling and vibrio load assay

During the challenge experiment, three clams were collected from each of the two challenge groups and one control group at 1 day post infection (DPI), 3 DPI, 4 DPI, and 6 DPI respectively. All clams were checked for survival during sample collection. The digestive tract (whole hepatopancreas and stomach) of the clams was aseptically sampled, and each sample was dissected into two parts. One part was immediately frozen in liquid nitrogen and stored at -80 °C until it was used for 16S rRNA sequencing. The other part was immediately stored on ice, weighed, placed into 500 μL sterile seawater, and homogenized as soon as possible. The homogenates were diluted an appropriate number of times, 50 μL was then plated onto vibrio-selective thiosulfate-citrate-bile-saltssucrose (TCBS) agar, and the CFU were counted after 24 h of cultivation at 25 °C.

2.4 PCR amplification and 16S rRNA sequencing

Microbial DNA of 15 samples (three samples per time point) was extracted using the CTAB/SDS method. Polymerase chain reaction (PCR) was performed with Phusion®High-Fidelity PCR Master Mix (New England Biolabs). The variable V4-V5 region of the bacterial 16S rRNA gene sequences were amplified with the primer pair 515F (5′GTGCCAGCMGCCGCGG3′) and 907R (5′CCGTCAATTCMTTTRAGTTT3′) with a barcode. All PCRs, purification, and library constructions were performed as previously described (Li et al., 2018). Finally, 250-bp paired-end reads were generated on the Illumina HiSeq2500 platform.

2.5 Pyrosequencing and data analysis

The raw sequence reads were filtered using the Quantitative Insights Into Microbial Ecology (QIIME) quality controlled process (Caporaso et al., 2010). Potential chimeric sequences were identified and removed by the UCHIME algorithm (Edgar et al., 2011). After filtering, the high-quality sequences obtained were used in the following bioinformatics analyses. The operational taxonomic units (OTUs) were clustered with a cutoff of 97% similarity using Uparse (Edgar, 2013). The representative (most abundant) sequence for each OTU was screened and then taxonomically classified in the SILVA database (http://www.arb-silva.de), using the Ribosomal Database Project classifier (at 80% confidence threshold) (Wang et al., 2007; Quast et al., 2013). Sequences classified to the genus Vibrio were then compared with 16S rRNA sequences of the cultured strain MM21 to identify that the OTU corresponded to the MM21 strain. Singleton OTUs and sequences identified as chloroplasts or mitochondria were removed from the analysis. All sequences obtained from this study were deposited in GenBank Sequence Read Archive (SRA; accession Nos. SRR8741669-SRR8741683).

2.6 Statistical analysis

The rarefaction curves and rank-abundance curves, which were plotted to determine the abundance of communities and sequencing data of each sample, were constructed using R (version 3.2.3, http://www.r-project.org). Alpha diversity (OTU richness, phylogenetic diversity, and Simpson index) and beta diversity were analyzed by QIIME with default parameters. Beta diversity across the 15 samples was calculated based on Bray-Curtis and weighted UniFrac distances (Hamady et al., 2010) and analyzed by nonmetric multidimensional scaling (NMDS), which showed the overall dissimilarity of the bacterial community in different infection stages. An unweighted pair group method with arithmetic mean (UPGMA) tree was also constructed with QIIME. Discriminatory analysis of taxonomic groups among infection stages was performed with LEfSe ( P ＜0.01, log 10 linear discriminant analysis score＞4) (Segata et al., 2011).

3 RESULT

3.1 Mortality under vibrio challenge

Clams in the challenge group were infected with V. parahaemolyticus by immersion infection. Low mortality was observed (1.09%) during the first 2 days, but mortality then drastically increased and the daily mortality reached 16.3% on the fifth day. Although the artificial vibrio challenge was stopped on day 5, the daily mortality peaked at 34.78% and 29.35% on day 6 and 7, respectively (Fig.1). Additionally, nearly 70% of surviving individuals on day 6 died by the following time point (day 7). Experimental challenge with V. parahaemolyticus resulted in 86.96% accumulated mortality of clams within 7 days, as opposed to no death events in the control group ( P ＜0.001). The results indicated that the controlled experiments of clam immersion in vibrio-contaminated seawater was an effective way to mimic the acute disease outbreak.

Fig.1 Mortality curves for clams challenged with vibrio

3.2 General characteristics of the clam digestive tract microbiome

We characterized the microbial dynamics and composition of the communities in the healthy clams (control group) and challenged clams at different time points. The key time points were 1, 3, and 4 DPI under constant vibrio immersion infection, as clams were healthy and only sporadically died in these stages; at 6 DPI, 1 day after stopping the artificial infection, a death outbreak occurred and the daily mortality peaked at this time point. Briefly, these time points covered different infection phases and were in consistent with those identified by a previous infection experiment (Yu et al., 2019a), these time points could be divided into latency phase (1 DPI, no dead clams found), prodrome phase (3 and 4 DPI, dead clams were sporadically observed), and onset phase (6 DPI, death outbreak occurred).

In total, 913 225 valid sequence reads binned into 1 733 OTUs were retrieved from 15 individual samples of the control and challenge groups over time. The rarefaction curves of most samples were nearly asymptotic, which indicated that the sequencing depth per sample covered most of the microorganisms (Fig.2). A total of 120 OTUs were shared by all five groups, whereas there were 33 unique OTUs in the control group, 125 in the 1 DPI group, 485 in the 3 DPI group, 131 in the 4 DPI group, and 16 in the 6 DPI group (Fig.3). Among the groups, the highest number of unique OTUs was observed in the 3 DPI group, which indicated that more specific microbial species appeared in this phase.

Fig.2 Rarefaction curves of 15 samples

Table 1 Average relative abundance and taxonomy information of the top10 OTUs from digestive tract samples at different stages of vibrio challenge

Fig.3 Numbers of shared and specif ic OTUs in the f ive groups that represent different infection stages

Fig.4 Relative abundance of bacterial phyla in digestive tract samples at different stages of vibrio challenge

The relative OTU abundances fit the Fisher’s logseries distribution well, which indicated that few OTUs accounted for the majority of the reads. For example, the relative abundance of OTU_1 (Alphaproteobacteria) was 49.78%±33.15% in the samples (Table 1). The phyla Proteobacteria and TM6-(Dependentiae) encompassed the bulk of the OTU diversity and abundance across all 15 samples from the five groups (Fig.4). In detail, 67.5% of the phylotypes belonged to Proteobacteria, and 14.3% belonged to TM6-(Dependentiae). Proteobacteria and TM6-(Dependentiae), the two core phyla, were relatively similarly distributed in each sample in each group, except in one specimen, 4-1, which was vastly different from the other specimens in the 4 DPI group. The core phylum composition showed significant changes over the course of vibrio challenge. There was a significant decrease in Proteobacteria abundance, and TM6-(Dependentiae) abundance markedly increased and became the dominant bacteria at 3 DPI. In the subsequent stage, the Proteobacteria abundance increased again and accounted for 99% of the total community at 6 DPI. However, the composition of Proteobacteria at 6 DPI became far less diverse, with OTU_1 (Alphaproteobacteria) accounting for nearly 97.75%±0.25% of the bacterial diversity.

Fig.5 Microbial diversity comparisons of digestive tract samples at different infection stages

3.3 Variation in microbial community composition changes during disease progression

During the course of vibrio challenge, a significant increase in the number of OTUs was observed at 3 DPI ( P ＜0.05), followed by a distinct decline at 6 DPI ( P ＜0.05). The phylogenetic diversity also significantly increased at 3 DPI compared with the control group (Fig.5a), although the microbial communities of 1 and 4 DPI exhibited no significant variation compared with that of the healthy control group based on the alpha diversity analysis (Fig.5a).The UPGMA tree based on weighted UniFrac distances, which measures the phylogenetic similarity between bacterial communities, showed that the microbiomes of the control group and 1 DPI group (initial time of challenge) clustered first, and the microbial communities in the 3 DPI group differed more from those of other groups (Fig.5b). This finding was consistent with previous alpha diversity results (Fig.5a). These results both demonstrated that a dramatic change in microbial communities in the digestive tract occurred at 3 DPI during vibrio challenge. Notably, a significantly low OTU diversity was found in the moribund samples at 6 DPI (Fig.5a), which indicated that the communities of the moribund samples were characterized by the disappearance of some rare OTUs and proliferation of one OTU (OTU_1).

Fig.6 NMDS plots of Bray-Curtis distances between the microbial communities showing sample status stability

Furthermore, we connected the samples from one group and used the area of the triangles in the NMDS plot to represent the relative stability of community composition at different stages during vibrio challenge (Fig.6). The microbiomes in the early stage of vibrio infection (e.g., 1 and 3 DPI), were more stable compared with the microbiomes at the end stage of vibrio infection. Moreover, the health status of the clams appeared divergence with the progression of the vibrio challenge. The daily mortality was very low during the first 3 days, but 17.44% of surviving samples at 4 DPI were dead by the next day. Therefore, the unstable microbial dynamics at 4 DPI and 6 DPI might be partially caused by the individual variations in health status at the end stage of vibrio infection.

3.4 Effects of vibrio infection and indicator taxa

Fig.7 Vibrio CFU counts on TCBS agar at different infection stages

In the samples from natural environments, vibrio was detected at a relative abundance of approximately 0.2%-1.8%. The microbiomes of samples (1 DPI, 3 DPI, and 4 DPI) immersed in vibrio-contaminated seawater were unexpectedly dominated by Alphaproteobacteria (which could not be further classified to the genus level) and phylum TM6-(Dependentiae), and not by vibrio. The highest relative abundance of vibrio (14.9%±8.8%) was found at 1 DPI (Table 1). However, the relative abundance of vibrio subsequently decreased to 2.2%±0.7% and 0.1%±0.06%, respectively, at 3 DPI and 4 DPI (Table 1), even if a high concentration of V. parahaemolyticus was maintained in the surrounding environment. The relatively high persistence of bacterial residents (e.g., Alphaproteobacteria) before and after the vibrio challenge might indicate that the microbiome was buffering the colonization and amplification of the pathogen vibrio. Furthermore, the absolute abundance of vibrio in the clam digestive tract (half of which were used for 16S rRNA gene amplicon pyrosequencing) was calculated by CFU counts (Fig.7). The vibrio colony counts were closely correlated with the relative read abundances in the 16S rRNA gene amplicon pyrosequencing analysis (Spearman’s ρ=0.65, P ＜0.01).

We next asked whether any of these bacterial taxa were strongly associated with the disease progression. LEfSe ( P ＜0.01, log 10 linear discriminant analysis score＞4) analysis was used to determine the indicator species (Fig.8). These analyses identified 33 discriminatory taxa, of which the majority (21) distinguished the 4 DPI group from the other groups. The discriminating taxon for the control group was Actinobacteria, despite this phylum appearing in the entire dataset. The genus vibrio could be considered an indicator at 1 DPI, which indicated that vibrio abundances were associated with infection but not directly with mortality or disease. The 3 DPI group was largely distinguished by the presence of the phylum TM6-(Dependentiae) and the genus Mycoplasma, whereas the 4 DPI group was distinguished by high abundances of the families Comamonadaceae and Planococcaceae, order Sphingobacteriales, and phylum Planctomycetes. Together, these indicator taxa might have the potential to distinguish health status among samples during vibrio infection.

Fig.8 Identif ied bacterial taxonomic biomarkers by LEfSe ( P ＜0.01, log10 linear discriminant analysis score＞4) analysis

4 DISCUSSION

Vibrio is a gram-negative bacterium that is widely distributed in marine environments and is recognized as one of the main opportunistic pathogens that causes mass mortality in cultured clams in summer (Yue et al., 2010). The population structure of family Vibrionaceae (‘vibrios’) among healthy marine invertebrates in coastal habitats has been explored, and the low specificity of animal-associated vibrio populations was detected by comparison of the vibrio populations collected from habitat water samples and animal samples (Preheim et al., 2011). This indicates that vibrios do not have a high degree of host preference, which supports the hypothesis that environmental dynamics are important factors in vibrio colonization. Therefore, we used vibrio immersion infection to mimic the vibrio outbreak in wild environment.

In this process, clams could easily filter exogenous bacteria from seawater, and mucus- and gill-associated microbiota may constitute a first line of defense against pathogens compared with infection via injection (Espinosa et al., 2016; Lokmer et al., 2016). Immersion infection by vibrio-contaminated seawater clearly caused approximately 87% mortality within 7 days. Faster and higher mortality occurred here compared with previous experimental infection (Liang et al., 2017), mainly because of higher seawater temperature in this study. Thus, laboratory immersion infection is helpful for investigating microbial diversity and dynamics during disease progression in the clam Meretrix petechialis.

Given the functional importance of clam digestive tract in nutrition digestion and innate immunity, it is necessary to investigate how the invading pathogenic vibrio changes the microbial communities of the digestive tract during disease progression. Therefore, in this study, we investigated the bacteria changes in the clam digestive tract when challenged with vibrio using a 16S rRNA sequencing approach. The phylum analysis revealed that Proteobacteria, TM6-(Dependentiae), Firmicutes Bacteroidetes, Actinobacteria, Tenericutes and Deinococcus-Thermus were the most common phyla in the healthy clam samples. Proteobacteria was the most abundant phylum and accounted for 58.57%±9.16% of the average relative abundance in the control group. The result was generally consistent with those of recent studies in which the resident microbiota of oyster species was characterized by 16S rRNA pyrosequencing (Fernández et al., 2014; Bernal et al., 2017).

In addition, our results revealed dramatic changes in the microbial communities during vibrio disease progression that preceded the mortality dynamics. At the end stage of vibrio infection, high mortality was accompanied by an obvious decrease in the microbial diversity of the digestive tract, and the overproliferation of OTU_1 (Alphaproteobacteria) indicated the collapse of digestive tract microbe homeostasis at 6 DPI. The decrease in the diversity of the digestive tract microbial community coincided with impaired health, and therefore could be as an indicator of health decline. This phenomenon has been found in various animals (Garnier et al., 2007; Chang et al., 2008; Green and Barnes, 2010). Lokmer and Wegner (2015) also found that dead or moribund oysters displayed signs of hemolymph microbiome structure disruption that was characterized by very low diversity and proliferation of a few taxa. Overall, our results indicate that the digestive tract microbial community is not merely the result of a filter-feeding lifestyle, that is, the microbiome is not simply determined by its surrounding environment. The stability of these digestive tract microfloras and their possible role in stimulating immunity or in competitive exclusion of external pathogens may have crucially affected the health status of clams during vibrio infection.

Host resistance refers to the ability to reduce pathogen replication and limit parasite growth through behavioral, morphological, and/or immunological mechanisms (Råberg et al., 2007). Therefore, the absolute number of pathogens (pathogen burden) is an important indicator for evaluating host resistance. Our study showed an obvious decrease in vibrio CFU numbers and relative vibrio abundance in the digestive tract at 3 DPI, which demonstrated that the digestive tract might take certain defense strategies to eliminate infectious agents (vibrio). Previous studies showed that oxidative burst and antimicrobial peptides release would be potent reactions for intestine immunity to eliminate pathogens in invertebrates (Ha et al., 2005; Schmitt et al., 2012; Yang et al., 2016). This could partially explained why vibrio pathogens do not undergo mass proliferation in later infection stages. However, further studies are needed to determine how clams discriminate between commensal and pathogenic bacteria during immune elimination. Furthermore, we found that the digestive tract microflora at 3 DPI showed larger adjustments in microbial composition following vibrio challenge. This shift in the composition of endogenous OTUs might help repress the proliferation of external pathogens (e.g., vibrio), on the other hand this shift might indicate that the previous homeostasis of microflora was destroyed. So, although the absolute and relative vibrio abundance decreased at 3 DPI and the following days, both the damage of activated immune defense imposed on host tissues and the disruption of microbial homeostasis still possibly resulted in massive mortalities happened at 5 DPI.

Dramatic changes took place in the digestive tract microflora during the transition from healthy to diseased clams, and one of the most important issues is whether sensitive assemblages could serve as indicators of host health. In this study, we found that more than 30 bacterial taxa had significantly different abundances between the control and challenged sample clams. ε-Proteobacteria are usually rare in coastal seawater and sediments (Campbell et al., 2011; Gobet et al., 2012), and they are rare in oyster stomach, gut and gill microbiota (King et al., 2012; Wegner et al., 2013; Fernández et al., 2014). Arcobacter spp. was also rare in our control samples. Therefore, the sharp increase in Arcobacter spp. abundance observed after 3 and 4 days of vibrio infection indicates an indirect but pronounced role in mortality. High abundance of Arcobacter spp. strains have also been reported in moribund oysters (Lokmer and Wegner, 2015), necrotic sponges (Fan et al., 2013), and starved abalones (Tanaka et al., 2004), which indicates that Arcobacter spp. can act as opportunistic pathogens when they occur in high enough densities (Olson et al., 2014). Similarly, Mycoplasma is usually considered pathogens, and was common but in low abundance in our control samples; Mycoplasma significantly increased with vibrio infection at 3 DPI. Therefore, it is possible that these opportunistic pathogens are commonly present in the digestive tract of healthy clams. The immune system (or physical barriers in the body) is mostly capable of fending them off when the dose is small, but opportunistic pathogens may become infectious and fatal when the hosts are immunosuppressed because of (a) biotic stresses (de Lorgeril et al., 2018).

5 CONCLUSION

To date, little is known about the bacterial community in the clam. Our results provided insight into the bacterial community composition changes when the clams were infected with a high concentration of a vibrio strain and revealed potential taxa that indicate the disease progression. Our findings will provide useful information for diagnosing, predicting, and preventing disease outbreaks in clam aquaculture.

6 DATA AVAILABILITY STATEMENT

RNA-seq data and amplicon sequences for microbiota analysis have been made available through the SRA database (accession Nos. SRR8741669-SRR8741683). Other data analyzed during this study are included in the supplementary file. Complementary information is available from the corresponding authors on reasonable request.